Technical  Report  GL-99-8 
August  1999 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


Improving  the  Traction  Prediction 
Capabilities  in  the  NATO  Reference 
Mobility  Model  (NRMM) 

by  Jody  D.  Priddy 


Approved  For  Public  Release;  Distribution  Is  Unlimited 

19991005  188 


DTIC  QUALITV  INSPECTiSD  4 

Prepared  for  Headquarters,  U.S.  Army  Corps  of  Engineers 


The  contents  of  this  report  are  not  to  be  used  for  advertising, 
publication,  or  promotional  purposes.  Citation  of  trade  names 
does  not  constitute  an  official  endorsement  or  approval  of  the  use 
of  such  commercial  products. 

The  findings  of  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position,  unless  so  desig¬ 
nated  by  other  authorized  documents. 


^  PRINTED  ON  RECYCLED  PAPER 


US  Army  Corps 

of  Engineers 

Waterways  Experiment 
Station 


HEMXXMRTERS 
BULOING 


FOR  INFORMATION  CONTACT: 

PUBLIC  AFFAIRS  OFFICE 
U.S.  ARMY  ENGINEER 
WATERWAYS  EXPERIMENT  STATION 
3909  HALLS  FERRY  ROAD 
VICKSBURG,  MISSISSIPPI  39180-6199 
PHONE:  (601)  634-2502 


STRUCTunES 

LABORATORY 


iWSAOF  BESOWATIOt)  -.r/KUT, 


Waterways  Experiment  Station  Cataloging-in-Publication  Data 

Priddy,  Jody  D. 

Improving  the  traction  prediction  capabilities  in  the  NATO  Reference  Mobility  Model 
(NRMM)  /  by  Jody  D.  Priddy ;  prepared  for  U.S.  Army  Corps  of  Engineers. 

215  p. :  ill. ;  28  cm.  —  (Technical  report ;  GL-99-8) 

Includes  bibliographic  references. 

1 .  Vehicles,  Military  — ^Traction  —  Computer  programs.  2.  NATO  Reference  Mobility 
Model  (NRMM)  3.  Tires,  Rubber  —  Traction  —  Computer  programs.  I.  United  States. 
Army.  Corps  of  Engineers.  II.  U.S.  Army  Engineer  Waterways  Experiment  Station. 

III.  Geotechnical  Laboratory  (U.S.  Army  Engineer  Waterways  Experiment  Station)  IV.  Title. 
V.  Series;  Technical  report  (U.S.  Army  Engineer  Waterways  Experiment  Station) ;  GL-99-8. 
TA7  W34  no.GL-99-8 


TABLE  OF  CONTENTS 


Page 

PREFACE . « 

ACRONYMS  . X 

NOMENCLATURE .  » 

CHAPTER 

1.  INTRODUCTION..... .  1 

Background .  ^ 

Objective  and  Scope .  ^ 

Definitions . ^ 

n.  DATABASES .  16 

General  Background  Information .  16 

Descriptive  Statistics .  1^ 

Trends  in  Tire  Shape  .  26 

m.  BASIC  MODELING  CONCEPTS. .  29 

Simplified  Tire/Soil  Interaction .  29 

Tire  Deflection  Effect .  21 

Model  Development  Tools .  24 

IV.  veil  MODELING .  40 

Overview .  40 

Existing  Models .  41 

New  Development  Models .  64 

Model  Comparisons .  67 


CHAPTER  Page 

V.  DRAWBAR  MODELING .  71 

Overview . 71 

Existing  Models .  72 

New  Development  Models .  83 

Model  Comparisons .  99 

VI.  SUMMARY  AND  CONCLUSIONS . 102 

Summary .  102 

Conclusions .  105 

Vn.  DIRECTIONS  FOR  FURTHER  RESEARCH .  107 

BIBLIOGRAPHY .  Ill 

APPENDIX 

A  TABLES  PRESENTING  DETAILED  DESCRIPTIONS 
OF  THE  OBSERVATIONS  USED  DURING 

THE  TRACTION  MODELING .  115 

B  SUMMARY  OF  ALL  MODEL  VARIANTS  EXPLORED 

DURING  IHE  veil  MODELING .  143 

C  SUMMARY  OF  ALL  MODEL  VARIANTS  EXPLORED 

DURING  THE  DRAWBAR  MODELING . 181 

D  TYPICAL  PROCEDURES  FOR  MEASURING 

veil  AND  DRAWBAR  PERFORMANCE . 


iv- 


LIST  OF  TABLES 


TABLE  Page 

1  General  statistics  for  the  VCIj  database .  19 

2  General  statistics  for  the  field  drawbar  database .  20 

3  General  statistics  for  the  laboratory  drawbar  database  .  21 

4  Quality  comparison  for  the  VCIi  modeling  results  .  68 

5  Comparison  of  exponents  between  best  CP2  and 

Numeric  relationships . 69 

6  Quality  comparison  for  the  drawbar  modeling  results  (field  data  only)  ...  100 

A1  Greneral  description  of  the  observations  used  for  the  VCIj  modeling .  116 

A2  Principal  vehicle  characteristics  associated  with  the  VCIi  observations  ....  121 

A3  General  description  of  the  field  observations  used  for 

the  drawbar  modeling .  126 

A4  Principal  vehicle  characteristics  associated  with  the 

field  drawbar  observations .  131 

AS  General  description  of  the  laboratory  observations  used  for  the  drawbar 

modeling  and  the  associated  principal  tire  characteristics  .  136 

B1  Quality  comparison  for  all  veil  model  variants  .  144 

Cl  Quality  comparison  for  all  drawbar  model  variants  (field  data  only)  .  182 


v- 


LISTOFFIGXJRES 


FIGURE  Page 

1  Prindpal  shape  variables  for  pneumatic  tires .  15 

2  WES  trafficability  equipment .  15 

3  Distribution  of  selected  system  variables  for  the  VCIi  database .  23 

4  Distribution  of  selected  system  variables  for  the  field  drawbar  database  ...  24 

5  Distribution  of  selected  system  variables  for  the  laboratory 

drawbar  database .  25 

6  Tire  shape  trends . 28 

7  Simple  stress  analysis  for  tire/soil  interaction  .  30 

8  Steady-state  force  analysis  for  tire/soil  interaction  .  31 

9  Long  chord  analogy  for  tire  contact  length  estimates  .  33 

10  Long  chord  versus  middle  ordinate .  34 

1 1  Standard  NRMM  relationship  for  VCIj  .  43 

12  Apparent  origin  of  the  Weight  Factor .  44 

13  Optimized  Weight  Factor  . .  45 

14  Optimized  NRMM  relationship  for  VCIj  .  47 

15  Relationship  between  Rowland’s  MMP  and  the  VCIj  data .  49 

16  Relationship  between  Larminie’s  MMP  and  the  VCIi  data  . .  51 

-vi- 


FIGURE  Page 

17  Relationship  between  Maclaurin’s  MMP  and  the  VCIj  data .  52 

18  Relationship  between  Freitag’s  Numeric  and  the  VCIj  data  .  54 

19  Basic  logic  behind  the  average  contact  pressure  parameter  CPI  .  57 

20  Best  CPI  relationship  for  VCIi .  .  60 

21  Basic  logic  behind  the  average  contact  pressure  parameter  CP2 .  62 

22  Best  CP2  relationship  for  VCIj . 64 

23  Best  numeric  relationship  for  VCIj .  66 

24  Similarity  between  best  CP2  and  Numeric  parameters  for  veil  .  70 

25  Standard  NRMM  relationship  for  drawbar .  75 

26  Optimized  NRMM  relationship  for  drawbar .  76 

27  Relationship  between  Freitag’s  Numeric  and  the  field  drawbar  data .  78 

28  Relationship  between  Tumage’s  Numeric  and  the  field  drawbar  data .  79 

29  Relationship  between  Smith’s  Numeric  and  the  field  drawbar  data .  81 

30  Relationship  between  Maclaurin’s  Numeric  and  the  field  drawbar  data  ....  83 

3 1  Relationship  between  new  numeric  variant  1  and  the 

field  drawbar  data . 

32  Relationship  between  new  numeric  variant  4  and  the 

field  drawbar  data .  90 

33  Relationship  between  new  numeric  variant  8  and  the 

field  drawbar  data .  92 

34  Comparison  between  laboratory  and  field  drawbar  measurements .  94 

3  5  New  numeric  variant  8  versus  Freitag’ s  Numeric  in 

laboratory  drawbar  data . 95 

-vii- 


HGURE 


Page 


36  New  numeric  variant  1 6  versus  Freitag’s  Numeric  in 

laboratory  drawbar  data .  97 

37  Relationship  between  new  numeric  variant  16  and  the 

field  drawbar  data .  99 

38  Sample  of  problems  due  to  modeling  drawbar  rather 

than  tractive  force  .  110 

D1  Sample  multi-pass  observations  used  for  determining  VCIj .  198 

D2  Sample  drawbar  data .  200 


-viii- 


PREFACE 


The  research  reported  herein  is  being  published  under  RDTE  Work  Unit  Number 
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Priddy.  Sincere  thanks  are  given  to  the  thesis  committee  members,  namely.  Dr.  Oswald 
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N  normal  force  at  tire/soil  interface 
■^3i  new  numeric  variant  3b  for  VCIj 
^c(F)  Freitag’s  numeric  for  wheels  in  clay  soils 

Freitag’s  numeric  for  wheels  in  clay  soils  cast  in  a  VCI  form 
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Q  torque  applied  to  the  wheel 

R  motion  resistance  force  at  tire/soil  interface 
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Sg  standard  error 

T  tractive  force  at  tire/soil  interface 

TCP  tire  construction  factor 

7EF  traction  element  factor  (part  of  MI) 

TF  transmission  factor  (part  of  MI) 
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one-pass  vehicle  cone  index 
fifty-pass  vehicle  cone  index 
average  axle  loading 
width  of  the  effective  contact  area 
wheel  load 

weight  factor  (part  of  MI) 
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hard-surface  tire  deflection 
normal  stress 
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shear  stress  component  in  the  X  direction 
shear  stress  component  in  the  Y  direction 
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CHAPTER  I 
INTRODUCTION 

Background 

A  ground-based  vehicle’s  ability  to  get  from  one  place  to  another  and  perform  work 
is  controlled  on  the  most  fundamental  level  by  its  traction  performance  characteristics.  A 
common  method  of  assessing  longitudinal  traction  performance  is  by  conducting  an  experi¬ 
ment  to  measure  the  drawbar  performance.  Drawbar  represents  a  vehicle’s  towing  ability  or 
its  acceleration  potential.  The  drawbar  force  (or  towing  force)  can  be  directly  measured 
during  a  controlled  experiment  whereas  the  tractive  force  cannot.  However,  the  drawbar 
force  is  equal  to  the  tractive  force  reduced  by  the  motion  resistance  forced  Therefore, 
drawbar  and  motion  resistance  forces  are  typically  the  parameters  of  interest  during  traction 
experiments. 

The  NATO  Reference  Mobility  Model  (NRMM)  was  developed  jointly  by  the  United 
States  Army  Tank-automotive  and  Armaments  Command  (TACOM)  and  the  United  States 
Army  En^eer  Waterways  Experiment  Station  (WES).  NRMM  was  developed  primarily  to 
predict  the  maximum  speed  a  vehicle  can  attain  over  various  terrain  conditions  and  environ¬ 
mental  scenarios  [Ahlvin  and  Haley  1992].  The  baseline  speed  prediction  generated  by  the 


*  -  A  section  containing  definitions  of  specialized  terms  is  provided  beginning  on  page  7. 
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model  is  based  on  traction  potential.  In  order  to  predict  traction  potential  in  off-road 
environments,  the  model  uses  relationships  for  predicting  drawbar  and  motion  resistance  as 
functions  of  vehicle  and  soil  characteristics. 

One  of  the  traction  relationships  in  NRMM  predicts  the  drawbar  potential  of  wheeled 
vehicles  operating  on  highly  plastic  clay  soil  (CH  classification  by  USCS^)  in  a  non-slippery 
state.  The  relationship  is  spedfically  applicable  to  all-drive  wheeled  vehicles  with  no  traction 
assist  devices  (e.g.,  tire  chains)  operating  on  level  terrain.  The  relationship  was  developed 
fi-om  the  results  of  several  traction  studies  which  involved  experiments  with  full-scale  vehicles 
in  a  field  environment  (i.e.,  naturally  occurring  soil  terrain).  These  experiments  were 
conducted  by  WES  during  the  1960s  through  the  1980s. 

The  NRMM  relationship  uses  the  One-Pass  Vehicle  Cone  Index  (VCIj)  to  predict 
drawbar  potential.  VCIj  is  the  minimum  soil  strength  necessary  to  support  a  vehicle  for  one 
pass  with  all  of  the  vehicle’s  traction  potential  being  expended  just  to  maintain  a  slow,  steady 
speed.  Therefore,  VCIj  represents  a  soil  strength  at  which  a  vehicle’s  drawbar  potential  is 
zero.  The  NRMM  relationship  makes  use  of  another  traction  relationship  that  was  developed 
specifically  for  predicting  VCIi.  WES  developed  the  VCIj  relationship  based  on  the  results 
of  soil  trafBcability  studies  with  full-scale  vehicles  conducted  during  the  1940s  through  the 
1960s,  and  WES  has  continuously  refined  the  VCIj  relationship  over  the  years  based  on  new 
information. 


^  -  Sections  describing  specialized  acronyms  and  nomenclature  used  throughout  this  document  are 
provided  on  pages  x  through  xiv. 
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TheNRMM  relationship  predicts  drawbar  potential  as  a  function  of  the  soil  strength 
above  veil  termed  the  excess  soil  strength  (i.e.,  excess  soil  strength  equals  the  soil  strength 
minus  VCIi).  Therefore,  the  NRMM  relationship  essentially  predicts  drawbar  potential  as  a 
function  of  the  soil  strength  only.  This  modeling  approach,  referred  to  as  the  VCI  Approach, 
works  fairly  well  because  many  wheeled  vehicles  have  similar  drawbar  performance  character¬ 
istics  above  their  VCIj  values  even  though  their  VCIi  values  can  differ  greatly.  Although  the 
approach  works  fairly  well,  it  gives  no  consideration  to  the  effects  that  important  vehicle 
characteristics  (e.g.,  tire  diameter,  tire  deflection,  etc.)  can  have  on  traction  performance 
above  VCIj.  The  relationship  would  be  much  more  beneficial  in  many  NRMM  applications 
if  it  included  these  effects. 

As  a  direct  consequence  of  this  shortcoming  in  the  VCI  Approach,  the  drawbar  rela¬ 
tionship  is  separated  into  two  prediction  categories.  One  of  the  categories  is  for  wheeled 
velucles  with  contact  pressure  factors  (CPF)  greater  than  4  psi,  and  the  other  is  for  wheeled 
vehicles  with  CPF  less  than  4  psi.  The  CPF  is  a  parameter  that  is  loosely  related  to  the 
average  contact  pressure  under  a  vehicle  on  a  non-yielding  surface.  There  are  two  separate 
categories  because  a  few  unusual  vehicles  tested  by  WES  [Schreiner  1971]  demonstrated 
significantly  better  drawbar  performance  than  other  vehicles  WES  had  tested.  The  vehicles 
demonstrating  better  performance  had  CPF  values  less  than  4  psi  whereas  all  the  other  vehi¬ 
cles  had  values  of  5  psi  or  higher.  The  form  of  the  NRMM  relationship  could  not  account  for 
the  observed  difference  in  performance,  and  hence  two  separate  equations  (or  curve-fits)  were 
developed  based  on  the  two  CPF  categories.  However,  it  has  never  been  substantially  verified 
that  differences  in  contact  pressure  caused  the  differences  in  drawbar  performance. 
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WES  has  attempted  to  develop  other  methods  of  predicting  drawbar  potential  that 
consider  the  effects  of  vehicle  characteristics  on  traction  performance.  A  large  research  effort 
was  conducted  during  the  1960s  and  1970s  in  which  laboratory  experiments  were  conducted 
with  a  single  tire  dynamometer  in  highly  processed  and  controlled  clay  and  sand  soils.  The 
primary  goal  of  this  research  effort  was  to  develop  a  single  dimensionless  parameter  (termed 
a  numeric)  that  could  be  used  to  predict  the  total  traction  performance  of  wheeled  vehicles 
[Tumage  1995],  This  approach,  referred  to  as  the  Numeric  Approach,  was  founded  on  a 
dimensional  analysis  of  tire/soil  interaction  [Freitag  1965],  and  the  VCIj  relationship  was  not 
considered  or  used  in  the  approach. 

The  Numeric  Approach  research  generated  a  lot  of  information  demonstrating  effects 
that  vehicle  characteristics  can  have  on  total  traction  performance,  but  the  primary  goal  was 
never  achieved.  One  significant  reason  the  laboratory  testing  was  abandoned  was  because 
relationships  developed  fi'om  the  single  tire,  laboratory  data  could  not  predict  real  vehicle 
performance  in  a  field  environment.  One  possible  reason  for  this  inability  to  predict  real 
vehicle  performance  was  that  the  processing  and  control  of  the  laboratory  soils  probably 
altered  them  to  such  a  high  degree  that  they  no  longer  had  the  trafficability  characteristics  of 
naturally  occurring  soils.  Another  possibility  arises  from  the  unnatural  process  in  which  the 
dynamometer  carriage  controlled  the  longitudinal  displacement  of  the  tire  rather  than  the 
more  realistic  process  of  the  tire/soil  interaction  controlling  the  longitudinal  displacement  of 
the  carriage.  There  is  no  proof  to  substantiate  either  of  these  explanations. 

In  a  recent  study  [Priddy  1995],  digital  databases  were  developed  for  use  in  a 
statistical  analysis  to  determine  the  inherent  variability  in  the  NRMM  traction  relationships. 
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The  development  of  the  digital  databases  involved  the  research  of  about  seventy  studies  in 
which  traction  experiments  were  conducted.  It  also  involved  a  detailed  exploration  of  the 
implemented  use  of  the  traction  relationships  in  NRMM.  The  databases  include  most  of  the 
ori^nal  data  used  in  the  development  of  the  current  NRMM  relationships  for  vehicle  traction 
on  fine-grained  soils  plus  some  data  fi'om  more  recent  (post  1980s)  traction  experiments. 
These  digital  databases  provide  a  readily  accessible  source  of  traction  performance  data  for 
real  vehicles  in  a  field  environment.  The  databases  also  provide  a  new  opportunity  for 
investigating  potential  improvements  to  the  traction  algorithms  in  NRMM. 

Traction  forms  the  baseline  limitation  on  speed,  but  it  is  not  the  only  thing  considered 
in  NRMM.  Other  speed-limiting  mobility  issues  are  considered  such  as  overriding  or  avoiding 
discrete  obstacles  and  vegetation,  climbing  or  descending  slopes,  limited  visibility  concerns, 
surfoce  roughness  effects  on  ride  quality,  and  tire  durability  concerns  due  to  large  operating 
deflections,  and  these  issues  don’t  include  the  on-road  issues  that  are  addressed  in  the  model. 
The  strength  of  NRMM  is  that  it  can  be  used  to  evaluate  the  mobility  of  one  or  more  vehicles 
based  on  all  these  mobility  concerns  using  real  terrain  data  fi'om  regions  of  consideration  that 
may  encompass  a  vast  amount  of  area.  Terrain  databases  are  mapped  for  use  with  NRMM, 
and  standards  have  been  established  for  the  types  of  terrain  data  that  must  be  used.  These 
standards  were  partly  the  result  of  practical  considerations  for  the  types  of  information  that 
could  effectively  be  obtained  for  vast  areas  of  terrain. 
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Objective  and  Scope 

The  goal  of  this  thesis  was  to  demonstrate  the  feasibility  of  utilizing  existing  perform¬ 
ance  data  to  develop  better  traction  relationships  for  NRMM.  The  approach  used  to  explore 
the  feasibility  was  to  conduct  a  statistical  modeling  effort  for  the  development  of  some  new 
relationships.  The  modeling  effort  encompassed  a  systematic  and  extensive  exploration  of 
new  and  old  ideas  from  both  the  VCI  and  Numeric  Approaches  along  with  some  new 
modeling  concepts. 

In  order  to  judge  whether  or  not  the  feasibility  existed,  a  set  of  decision  criteria  related 
to  the  quality  of  the  new  relationships  was  established.  First,  the  relationships  should  be 
lo^cal,  have  a  theoretical  foundation  where  possible,  and  be  dimensionally  correct.  Second, 
they  should  provide  a  more  sound  understanding  of  the  role  that  vehicle  characteristics  have 
in  traction  performance  throughout  the  full  range  of  soil  strengths  from  VCIj  up  to  maximum 
strengths.  Last,  they  should  describe  (i.e.,  collapse  into  a  predictable  trend)  the  existing  data 
as  well  or  better  than  the  current  NRMM  relationships. 

The  statistical  modeling  effort  focused  on  two  major  types  of  traction  relationships: 
(1)  veil  and  (2)  drawbar.  Three  types  of  traction  data  were  explored:  (1)  VCIj  data  from 
experiments  with  real  vehicles  in  a  field  environment,  (2)  drawbar  data  from  experiments  with 
real  vehicles  in  a  field  environment,  and  (3)  drawbar  data  from  experiments  with  a  single  tire 
dynamometer  in  a  laboratory  environment.  The  data  were  limited  to  experiments  with 
pneumatic  tires  in  non-slippery,  highly  plastic  clay  soils  (CH  classification  by  USCS). 

There  are  two  sides  of  the  vehicle/terrain  interaction  ^stem  that  must  be  characterized 
in  some  fashion  for  any  relationship  designed  to  predict  vehicle  traction:  (1)  the  terrain  and 
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(2)  the  vehicle.  This  research  concentrated  on  more  definitively  characterizing  the  vehicle 
(i.e.,  loading)  side  of  the  interaction  system  for  tires  in  soil.  A  simple,  yet  extremely  usable, 
method  of  characterizing  the  soil  (i.e.,  terrain)  side  of  the  interaction  system  has  been  devel¬ 
oped  and  standardized  by  WES.  WES’  method  involves  soil  characterization  by  USCS  soil 
type  and  rating  cone  index  (RCI)  (see  definitions  for  details).  This  method  produces 
prediction  accuracy  that  is  adequate  for  NRMM  requirements  and  has  thus  been  standardized 
as  the  soil  characterization  method  for  use  with  the  model.  Therefore,  WES  standard 
method  of  terrain  characterization  was  employed  in  this  research. 

Definitions 

The  following  are  definitions  of  specialized  terms  used  throughout  this  document. 
Pneumatic  Tire  Terms  fsee  Figure  H 

r.arr.agg  The  portion  of  a  tire  outside  the  lip  of  the  rim  flange  that  is  free  to  deflect. 

Diameter  rdV  The  outside  diameter,  not  counting  the  tread  height,  of  an  inflated  but 

unloaded  tire. 

Sertinn  height  ThV  The  distance  from  the  lip  of  the  rim  flange  to  the  outer  tire 

surface,  not  counting  the  tread  height,  of  an  inflated  but  unloaded  tire. 

.Serrinn  width  rhV  The  distance  between  the  extreme  points  on  the  tire  carcass  area 


of  an  inflated  but  unloaded  tire. 
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Deflection  (b\  The  distance  between  the  outer  surface  of  a  loaded  tire  and  the  outer 
sur&ce  of  an  inflated  but  unloaded  tire  taken  at  the  center  of  the  carcass  cross  section 
on  a  non-yielding  surface. 

Tread  hdght.  The  thickness  of  the  aggressive  features  on  the  outer  surface  of  a  tire. 
Aspect  ratio.  The  ratio  of  section  height  over  section  width. 

Vehicle  Terms 

All-drive.  Indicates  that  all  of  the  axles  (i.e.,  tires)  are  powered. 

Contact  pressure  factor  rCPFV  A  factor  that  is  loosely  related  to  the  average  hard- 
surface  contact  pressure  under  a  vehicle.  It  is  one  of  eight  factors  used  to  calculate 
the  Mobility  Index. 

Mobility  Index  fMTl.  A  parameter  that  is  related  to  the  VCI  performance  of  vehicles 
on  fine-grained  soils.  It  was  developed  in  the  United  States. 

Mean  Maximum  Pres.sure  fMMPV  A  parameter  that  is  related  to  the  average  ground 
contact  pressure  under  a  vehicle  and  is  used  to  represent  soft-soil  performance 
potential.  It  was  developed  in  the  United  Kingdom. 

Gross  vehicle  weight  rGVWl.  The  total  weight  of  a  vehicle. 
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Soil  Terms 

Trafficahilitv.  The  ability  of  terrain  to  support  the  passage  of  vehicles  [Meyer  et  al. 
1977]. 

r.nnft  Index  rcn  An  index  of  soil  shear  strength  obtained  using  a  trafficability  cone 
penetrometer  standardized  by  WES. 

T^emnlH  Index  fRIV  An  index  of  the  sensitivity  of  soil  to  strength  losses  under 
vehicular  traffic  obtained  using  remolding  equipment  standardized  by  WES. 

I?  atinp  r.nne  Index  (RCX)  An  index  of  soil  shear  strength  that  includes  consideration 
of  the  sensitivity  of  soil  to  strength  losses  under  vehicular  traffic.  It  is  defined  as  the 
product  of  cone  index  and  remold  index  for  the  particular  layer  of  soil.  For  example, 
the  6-12"  RCI  equals  the  6-12"  RI  times  the  6-12"  Cl. 

Cnne  penetrometer  Tsee  Figure  21.  An  instrument  consisting  of  a  circular  cylindrical 
shaft  (usually  18  to  36  inches  in  length)  with  a  30-degree  right  circular  cone  mounted 
on  one  end  and  a  calibrated  load  measuring  device  on  the  other  end.  For  fine-grained 
soils,  the  shaft  is  3/8  inch  diameter  steel  (previously  5/8  inch  aluminum)  and  the  cone 
has  a  0.5-square-inch  base  area.  The  output  measurement  (Cl)  is  the  average  of 
pressure  readings  (typically  in  pounds  per  square  inch)  taken  at  specified  depths  of 
penetration  of  the  base  of  the  cone  into  the  soil.  The  depths  of  penetration  used  in  the 
measurement  are  usually  those  taken  at  the  top,  mid-height,  and  bottom  of  the  critical 
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layer.  The  pressure  readings  are  the  result  of  the  penetration  force  divided  by  the  base 
area  of  the  cone  with  a  standard  penetration  rate  of  72  inches  per  minute. 

Remolding  equipment  (see  Figure  21.  Equipment  consisting  of  a  circular  cylindrical 
tube  mounted  on  a  steel  base  and  a  drop  hammer.  The  tube  has  1-7/8  inch  inside 
diameter.  The  drop  hammer  weighs  2-1/2  Ibf  and  has  12  inches  of  drop  travel.  In 
use,  soil  samples  approximately  6  inches  in  height  are  inserted  into  the  tube  using  a 
trafiBcability  (or  Hvorslev)  sampler.  For  fine-grained  soils,  cone  index  measurements 
are  taken  in  the  center  of  the  sample  before  and  after  1 00  blows  of  the  drop  hammer. 
The  cone  index  measurements  are  based  on  readings  taken  in  the  sample  at  depths  of 
0, 1, 2, 3,  and  4  inches.  The  output  measurement  (RI)  is  the  ratio  of  the  cone  index 
measurement  after  100  blows  over  the  cone  index  measurement  before  100  blows. 

Trafficability  for  Hvorslev^  sampler  (see  Figure  2).  A  piston-type  sampling  device 
that  is  used  to  obtain  an  undisturbed  sample  in  soft  soils.  It  has  a  circular  cylindrical 
tube  with  1-7/8  inch  inside  diameter  that  is  sharpened  on  the  open  end.  The  piston 
within  the  tube  retracts  during  penetration  into  the  soil  such  that  a  partial  vacuum  is 
maintained  above  the  sample  preserving  the  soil’s  in  situ  structure. 

Critical  layer.  A  layer  of  soil  lying  below  the  natural  terrain  surface  that  exerts  the 
greatest  influence  on  trafficability.  The  depth  of  the  critical  layer  is  dependent  upon 
vehicle  characteristics  and  the  nature  of  the  cone  index  (Cl)  profile  with  depth.  A  6 
inch  layer  of  soil  is  typically  used,  but  sometimes  a  12  inch  layer  of  soil  is  used.  It  is 
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typically  the  3-to-9  inch  layer  for  light  wheeled  vehicles  (wheel  loads  of 2,000  Ibf  and 
less)  and  the  6-to-12  inch  layer  for  normal  wheeled  vehicles  (wheel  loads  up  to  about 
10,000  Ibf). 

TTnified  Soil  Hassification  System  fUSCSV  A  system  which  identifies  (classifies)  soils 
according  to  their  textural  and  plasticity  qualities  and  to  their  grouping  with  respect 
to  their  performances  as  engineering  construction  materials  [Meyer  et  al.  1977]. 

Vehicle/Terrain  Interaction  Terms 

Traction.  The  process  by  which  a  ground-based  vehicle  develops  tractive  force  and 
overcomes  motion  resistance  to  produce  desired  motion  relative  to  the  terrain. 

Tractive  force  m  The  force  developed  at  the  vehicle/terrain  interface  by  the  traction 
elements  as  a  result  of  applied  torque  from  the  power  plant. 

Motion  resistance  fRV  Any  force  imposing  resistance  against  desired  motion.  For 
element-level  traction  considerations,  it  is  composed  of  rolling  resistance  forces  only. 

Rolling  resistance.  Motion  resistance  that  arises  fi’om  deformations  in  the  terrain 
(external)  and  the  traction  elements  (internal). 

Traction  element  Any  element  of  a  vehicle  that  is  designed  to  provide  support  and/or 
traction  for  a  vehicle  traveling  on  a  surface  (e.g.  tires,  tracks,  feet,  screws,  etc.) 
[Meyer  et  al.  1977]. 
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Drawfaar  (D).  The  amount  of  sustained  towing  force  a  self-propelled  vehicle  can 
produce.  It  is  the  resultant  of  tractive  force  reduced  by  motion  resistance. 

Drawbar  coefficient  CD/W).  Drawbar  divided  by  gross  vehicle  weight. 

Vehicle  Cone  Index  fVCIV  Minimum  soil  strength  necessary  for  a  vehicle  to  make 
a  specified  number  of  passes.  Consideration  is  most  often  given  to  one  pass  (VCIi) 
and  fifty  passes  (VCI50). 


Sinkage  (z).  The  depth  to  which  the  traction  elements  penetrate  the  terrain  measured 
normal  to  the  original,  undisturbed  surface. 


♦ 

Slip.  An  indication  of  how  the  speed  of  the  traction  elements  differs  fi'om  the  forward 


speed  of  the  vehicle.  It  is  defined  by  the  equation  [Meyer  et  al.  1977]: 


Slip  = 


''  r^co  - 
r^(o 


where: 


=  rolling  radius 

a>  =  angular  velocity  of  the  wheel  or  number  of  revolutions 
per  unit  time  divided  by  Itz  for  a  track 
V  =  forward  velocity  of  vehicle  or  wheel  axle 


Optimum  slip 


Slip  at  which  maximum  work  index  occurs. 


Work  Index  AVD.  A  dimensionless  number  that  represents  the  relative  efficiency  for 
a  particular  measure  of  drawbar.  It  is  defined  by  the  equation; 

Work  Index  =  — f  1  - 

W\  100  j 


-13- 


Maximum  pull  slip.  Slip  at  which  maximum  drawbar  occurs. 


Statistical  Modeling  Terms 


rnftffident  of  determination  A  measure  of  quality  that  indicates  the  amount  of 

variation  in  the  measurements  (Y)  that  is  accounted  for  by  the  relationship  predictions 
(7).  It  is  defined  as  the  ratio  of  the  explained  variation  in  Y  (i.e.,  the  variation 
explained  by  Y)  over  the  total  variation  in  Y,  and  it  is  usually  expressed  in  percent. 
It  can  be  calculated  using  the  equation: 

E  {r,  -  r,Y 

=  l  -  - 

E  (r,  -  yy 

/=l' 

where:  i  =  increment  counter  for  the  observations 
n  =  total  number  of  observations 
Y.  =  measured  value  for  i-th  observation 

Y.  =  predicted  value  for  i-th  observation 

Y  =  mean  measurement  for  all  observations 


Degrees  of  freedom.  A  quantity  related  to  quality  that  is  equal  to  the  total  number  of 
observations  less  the  number  of  empirical  constants. 


Residual.  The  difference  between  measured  (Y)  and  predicted  ( Y)  values  for  an 
individual  observation. 


Standard  ftrrnr  A  measure  of  quality  that  estimates  the  standard  deviation  of  the 
measurements  relative  to  the  equation  describing  the  relationship.  It  represents  the 
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amount  of  data  scatter  around  the  prediction  equation.  It  can  be  calculated  using  the 
equation; 


>=1 ' _ ' 

n  -  k 

where:  i  =  increment  counter  for  the  observations 
n  =  total  number  of  observations 

k  =  number  of  empirical  constants  in  the  prediction  equations 
Yj  =  measured  value  for  i-th  observation 

A 

7.  =  predicted  value  for  i-th  observation 


Adjusted  coefficient  of  determination.  A  measure  of  quality  that  is  identical  to 
with  the  exception  that  it  takes  the  degrees  of  freedom  associated  with  the  rela¬ 


tionship  into  account.  This  measure  of  quality  is  more  appropriate  when  comparing 
various  relationships  that  have  different  degrees  of  freedom.  It  can  be  calculated 


using  the  equation; 


where:  i  =  increment  counter  for  the  observations 
n  =  total  number  of  observations 

k  =  number  of  empirical  constants  in  the  prediction  equations 
Y.  =  measured  value  for  i-th  observation 

Y.  =  predicted  value  for  i-th  observation 

Y  =  mean  measurement  for  all  observations 


Figure  1.  Principal  shape  variables  for  pneumatic  tires 


Figure  2.  WES  trafScability  equipment 


CHAPTER  n 


DATABASES 

General  Background  Information 

As  described  in  the  introduction,  a  recent  study  [Priddy  1995]  produced  a  readily 
accessible  source  of  traction  data  from  experiments  with  real  vehicles  in  a  field  environment. 
The  digital  databases  generated  during  that  study  were  used  as  the  sole  source  of  field 
performance  data  in  this  research.  Only  two  of  the  major  databases  were  required:  the  Fine- 
Grained  One-Pass  Vehicle  Cone  Index  database,  and  the  Fine-Grained  Drawbar  at  Nominal 
Slip  database.  The  vehicle  characteristics  associated  with  the  performance  data  were  not 
originally  included  in  these  databases.  Therefore,  the  vehicle  characteristics  (primarily  tire 
dimenaons)  w^e  determined  and  added  as  part  of  this  research.  The  vehicle  characteristics 
were  determined  from  the  associated  reports  and  from  vehicle  files  maintained  by  WES 

The  wheeled  portion  of  the  Fine-Grained  One-Pass  Vehicle  Cone  Index  database 
includes  all  of  the  wheeled  vehicles  tested  by  WES  for  which  an  estimate  of  the  VCIj  could 
be  determined.  Some  of  the  data  were  from  fifty-pass  vehicle  cone  index  (VCIjo)  experiments 
and  some  had  vehicles  with  non-powered  traction  elements  (e.g.  towed  trailers).  VCI50  is  the 
nunimum  soil  strength  necessary  to  support  a  vehicle  for  fifty  passes,  and  VCI50  experiments 
did  not  usually  provide  sufficient  root  multi-pass  test  data  to  determine  VCIj  with  adequate 
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confidence  for  this  study.  Therefore,  the  observations  originating  from  VCI50  experiments 
were  reviewed,  and  those  having  low  confidence  estimates  were  not  considered.  The  obser¬ 
vations  having  non-powered  traction  elements  were  not  considered  either  since  the  VCIj 
relationships  are  for  all-drive  vehicles  only.  NRMM  has  algorithms  to  handle  vehicles  with 
non-powered  traction  elements,  but  these  algorithms  employ  other  traction  relationships  (i.e., 
drawbar  and  motion  resistance)  in  addition  to  the  VCIj  relationships.  The  final  database  used 
in  the  VCIj  modeling  included  experiments  that  spanned  a  period  of  thirty  years  fi’om  the 
1960s  to  the  1990s. 

The  wheeled  portion  of  the  Fine-Grained  Drawbar  at  Nominal  Slip  database  includes 
most  wheeled  vehicles  tested  by  WES  through  the  late  1980s  and  some  vehicles  fi'om  more 
recent  experiments.  The  database  includes  drawbar  performance  data  that  are  based  on  two 
critical  slip  zones:  (1)  20  percent  slip  or  optimum  slip  (typically  15  to  30  percent)  and  (2)  100 
percent  slip  or  maximum  pull  slip  (t5T)ically  20  to  100  percent).  Although  the  current  NRMM 
drawbar  relationships  are  based  on  100  percent  slip,  the  data  from  the  20  percent  slip  zone 
were  used  during  the  drawbar  modeling.  The  data  fi’om  the  100  percent  slip  zone  were  not 
used  because  the  performance  data  were  relatively  limited  and,  more  importantly,  because 
there  is  always  less  quality  control  in  drawbar  testing  at  the  100  percent  slip  zone  than  at  the 
20  percent  slip  zone.  The  final  field  database  used  in  the  drawbar  modeling  included 
experiments  that  spanned  a  period  of  thirty  years  fi’om  the  1960s  to  the  1990s. 

In  order  to  explore  the  difference  in  drawbar  performance  measurements  obtained  by 
WES  in  field  versus  laboratory  environments,  some  laboratory  drawbar  data  were  also 
explored.  As  will  be  shown  in  the  discussion  of  the  drawbar  modeling  (Chapter  V),  these  data 
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also  pro^dded  an  enhancement  to  the  drawbar  modeling  effort.  The  data  selected  for  use  were 
those  reported  by  Freitag  [1965]  in  his  dimensional  analysis  of  tire/soil  interaction.  This  was 
the  dimensional  analysis  research  that  laid  the  foundation  for  the  Numeric  Approach  research 
effort  described  in  the  background  portion  of  the  introduction  (Chapter  I).  The  data  contain¬ 
ed  experiments  with  four  different  treadless  tires  in  soil  bins  containing  clay  soils  that  were 
essentially  fnctionless  and  saturated. 

Descriptive  Statistics 

For  relationships  that  are  established  using  statistical  modeling  techniques,  it  is  imper¬ 
ative  that  the  limits  of  the  experimental  data  be  observed.  Use  of  these  relationships  outside 
the  experimental  limits  is  always  risky,  especially  if  there  is  no  underlying  theoretical  foun¬ 
dation.  For  this  reason,  tables  were  generated  that  show  the  total  number  and  types  of  obser¬ 
vations  and  the  minimum,  mean,  and  maximum  values  of  the  principal  modeling  variables. 
The  relative  magnitudes  of  variables  are  just  as  important  in  many  cases  as  the  absolute 
magnitudes.  Therefore,  the  minimum,  mean,  and  maximum  values  for  several  ratios  of  the 
principal  variables  were  also  included  in  the  tables.  This  was  especially  important  for  tire 
deflection,  a  control  variable  usually  spedfied  (through  recommended  tire  inflation  pressures) 
as  a  percent  of  the  tire  section  height  for  particular  operating  scenarios  such  as  highway, 
cross-country,  mud/sand/snow,  or  emergency. 

Table  1  shows  the  descriptive  statistics  for  the  VCI,  database.  As  shown,  there  were 
79  total  observations  used  in  the  VCIj  modeling.  About  40  percent  of  the  observations  had 
radial  tires  with  the  majority  having  bias  ply  tires,  and  there  were  two  observations  that  had 
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dual  tires  on  the  rear.  The  gross  vehicle  weights  ranged  from  less  than  1,200  to  over  86,000 
pounds,  and  these  loads  were  carried  on  anywhere  from  2  to  5  axles  translating  into  axle 
loadings  of  about  600  to  20,000  pounds.  Tire  diameters  ranged  from  17  to  63  inches  with 
section  widths  ranging  between  7  and  40  inches,  and  the  tire  shape  ratio  of  width-over¬ 
diameter  ranged  from  0.25  to  0.95.  The  deflections  expressed  in  percent  of  tire  section  height 
ranged  from  less  than  10  percent  to  nearly  50  percent  and  were  25  percent  on  average.  The 
measured  VCIj  values  ranged  from  extremely  low  at  5  to  fairly  high  at  40  and  were  22  on 
average. 


Table  1 .  General  statistics  for  the  VCIj  database 


Total  Observations:  79 


Radials:  3 1 


Bias  Ply:  48 


Duals:  2 


Singles:  77 


Vehicle 

Dimensions 

Average  Axle  Dimensions 

Meas. 

VCI, 

psi 

GVW 

Ibf 

D 

Q 

w 

Ibf 

n 

n 

b 

in. 

h 

in. 

5 

in. 

Min: 

1190 

2 

9.0 

595 

2 

17.4 

7.1 

4.5 

0.37 

5 

Mean: 

26038 

3 

13.2 

7969 

2 

41.7 

14.5 

10.8 

2.79 

22 

Max: 

86255 

5 

22.0 

20199 

3.333 

62.5 

40.0 

16.7 

5.93 

40 

5 

h 

5 

d 

H 

D 

B 

B 

K 

h 

Min: 

0.060 

0.016 

0.298 

0.248 

0.223 

0.216 

0.857 

Mean: 

0.250 

0.064 

0.765 

0.365 

0.257 

0.334 

1.304 

Max: 

0.467 

0.120 

0.944 

0.952 

0.317 

0.603 

2.340 
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Table  2.  General  statistics  for  the  field  drawbar  database 


Total  Observations:  84 


Radials:  24 


Bias  Ply:  60 


Duals:  0 


Singles:  84 


Vehicle 

Dimensions 

Average  Axle  Dimensions 

Meas. 

RCI 

psi 

Meas. 

D 

w,. 

GVW 

Ibf 

D 

w 

Ibf 

n 

d 

in. 

b 

in. 

h 

in. 

6 

in. 

Min: 

1900 

2 

950 

2 

22.5 

8.5 

6.2 

0.93 

9 

0.04 

Mean: 

18785 

3 

5910 

2 

10.7 

2.41 

199 

0.47 

Max: 

66000 

5 

16500 

2 

62.5 

750 

0.72 

5 

h 

6 

d 

B 

B 

h 

d 

0.130 

0.033 

0.300 

0.257 

0.230 

Mean: 

0.229 

0.059 

0.718 

0.416 

0.262 

Max: 

0.385 

0.097 

0.927 

0.952 

0.292 

Table  2  shows  the  descriptive  statistics  for  the  field  drawbar  database.  As  shown, 
there  were  84  total  field  observations  used  in  the  drawbar  modeling.  About  30  percent  of  the 
observations  had  radial  tires  with  the  vast  majority  having  bias  ply  tires,  and  none  of  the 
observations  had  dual  tires.  The  gross  vehicle  weights  ranged  from  1,900  to  66,000  pounds, 
and  these  loads  were  carried  on  anywhere  from  2  to  5  axles  translating  into  axle  loadings  of 
about  900  to  17,000  pounds.  Tire  diameters  ranged  jfrom  22  to  63  inches  with  section  widths 
ranging  between  8  and  40  inches,  and  the  tire  shape  ratio  of  width-over-diameter  ranged  fi'om 
0.26  to  0.95.  The  deflections  expressed  in  percent  of  tire  section  height  ranged  from  13 
percent  to  nearly  40  percent  and  were  23  percent  on  average.  The  measured  drawbar 
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performance  values  expressed  as  a  percent  of  gross  vehicle  weight  at  20  percent  slip  ranged 
from  near  zero  to  0.72.  These  drawbar  values  were  obtained  on  CH  soils  that  essentially 
spanned  the  entire  spectrum  of  soil  strengths  varying  from  a  very  soft  9  to  an  extremely  hard 
750  RCI. 


Table  3.  General  statistics  for  the  laboratory  drawbar  database 


Total  Observations  119 


Radials  0 


Bias  Ply  1 19 


Average  Axle  Dimensions 

Meas. 

Cl 

psi 

Meas. 

D 

w 

Ibf 

n 

d 

in. 

b 

in. 

h 

in. 

5 

in. 

Min: 

100 

1 

14.1 

4.2 

3.1 

0.26 

14 

-0.11 

Mean: 

472 

1 

25.4 

5.9 

4.6 

1.16 

39 

0.43 

Max: 

892 

1 

28.3 

8.3 

6.4 

2.24 

67 

1.51 

6 

h 

d 

B 

B 

h 

d 

Min: 

0.080 

0.009 

0.738 

0.150 

0.114 

Mean: 

0.254 

0.046 

0.771 

0.239 

0.184 

Max: 

0.450 

0.079 

0.803 

0.298 

0.226 

Table  3  shows  the  descriptive  statistics  for  the  laboratory  drawbar  database.  As 
shown,  a  total  of  119  laboratory  observations  were  explored  in  support  of  the  drawbar 
modeling.  There  were  no  laboratory  observations  with  radial  tires.  The  tire  loadings  ranged 
from  100  to  about  900  pounds.  These  tire  loadings  would  hypothetically  translate  to  realistic 
axle  loadings  of 200  to  1800  pounds  which  are  much  smaller  than  those  associated  with  the 
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field  (l.e.,  full-scale  vehicle)  observations.  Tire  diameters  ranged  from  14  to  28  inches  with 
section  widths  ranging  between  4  and  8  inches,  and  the  tire  shape  ratio  of  \wdth-over- 
diameter  ranged  Jfrom  0. 15  to  0.30.  The  deflections  expressed  in  percent  of  tire  section  height 
ranged  fi’om  8  percent  to  45  percent  and  were  25  percent  on  average.  The  measured  drawbar 
performance  values  expressed  as  a  percent  of  tire  loading  at  20  percent  slip  ranged  fi’om 
below  zero  to  an  astounding  1.51  (this  magnitude  of  drawbar  has  never  been  demonstrated 
with  real  vehicles).  The  soil  strengths  in  the  laboratory  prepared  CH  soils  ranged  fi’om  14  to 
67  Cl.  These  Cl  strengths  were  probably  equivalent  to  RCI  strengths  for  the  laboratory  soHs 
since  they  were  likely  fully  remolded  during  the  high  level  of  processing  necessary  to  prepare 
the  soil  bins  for  testing.  Therefore,  these  strengths  can  be  directly  compared  to  the  RCI 
strengths  fi’om  the  field  observations  for  comparing  the  strength  ranges. 

Another  important  concern  when  developing  relationships  with  statistical  modeling 
techraques  is  the  distribution  of  values  within  the  experimental  range  for  the  system  variables. 
It  is  desirable  to  have  values  distributed  throughout  the  experimental  range  rather  than  having 
a  small  number  of  clusters.  In  order  to  evaluate  the  distribution  of  values  for  the  system 
variables  in  the  three  databases,  dot  plots  were  generated  for  selected  variables. 

Figure  3  shows  the  dot  plots  for  the  VCIj  database,  and  it  demonstrates  that  there 
were  good  distributions  for  all  of  the  variables  except  for  section  width  which  had  nearly  all 
of  the  values  in  the  lower  half  of  the  experimental  range.  Figure  4  shows  the  dot  plots  for  the 
field  drawbar  database,  and  it  demonstrates  that  although  there  were  fair  distributions  for 
most  of  the  variables,  the  values  were  more  sparse  than  those  in  the  VCIj  database.  Figure 
5  shows  the  dot  plots  for  the  laboratory  drawbar  database,  and  it  demonstrates  that  although 
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Figure  4.  Distribution  of  selected  system  variables  for  the  field  drawbar  database 
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the  tire  diameter,  section  width,  and  section  height  variables  were  very  limited  relative  to  the 
field  databases,  the  load,  deflection,  soil  strength,  and  performance  metric  variables  had  good 
distributions. 

Trends  in  Tire  Shape 

When  statistical  modeling  techniques  are  used  to  develop  relationships,  it  is  imperative 
that  interrelationships  between  the  system  variables  be  understood.  In  this  research,  the  key 
task  in  the  modeling  eflfort  was  to  find  the  best  way  of  organizing  the  principal  vehicle 
variables  (i.e.,  w,  d,  b,  h,  and  6)  into  a  loading  parameter  that  had  a  strong  correlation  with 
the  performance  metrics.  Therefore,  interrelations  between  these  variables  were  investigated. 
Tire  load  (i.e.,  Wv^r  =  w/n)  and  deflection  are  control  variables  for  wheeled  vehicles,  but 
diameter,  section  width  and  section  height  are  “fixed”  variables  for  tires  that  could  possibly 
follow  trends  based  on  standards  in  tire  design,  so  the  interrelations  between  these  three  were 
explored. 

Figure  6  shows  three  x-y  plots  that  illustrate  the  trends  in  tire  shape  for  all  three  of  the 
databases  combined.  Figure  6a  shows  that  there  was  a  strong  linear  relationship  between 
diameter  and  section  hei^t  for  all  of  the  tires  in  the  three  databases  with  the  exception  of  one 
laboratory  tire  (4.00-20)  whose  section  hdght  was  very  small  relative  to  its  diameter.  Figure 
6b  shows  that  there  was  a  linear  relationship  between  diameter  and  section  width  for  most  of 
the  tires  in  the  three  databases,  but  several  tires  did  fall  far  away  from  the  normal  trend. 
Figure  6c  shows  that  there  was  also  a  linear  relationship  between  section  height  and  section 
vtidth  for  most  of  the  tires  in  the  three  databases.  Several  tires  fi’om  field  experiments  fell  far 
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away  from  the  normal  section  height-to-section  width  trend,  but  all  the  tires  from  the 
laboratory  experiments  fell  within  the  normal  trend. 


Diameter,  d  (in.)  P  Diameter,  d  (in.) 


Figure  6.  Tire  shape  trends 


CHAPTER m 


BASIC  MODELING  CONCEPTS 

Simplifled  Tire/Soil  Interaction 

For  a  locomotive,  there  is  the  rail;  for  an  airplane,  there  is  the  air;  and  for  a  wheeled 
vehicle,  there  is  the  road,  unless  of  course,  the  road  is  not  an  option.  Such  is  the  case  for 
wheeled  vehicles  designed  for  industries  like  military,  mining,  construction,  agriculture,  and 
so  forth.  For  this  breed  of  wheeled  vehicle,  which  is  designed  primarily  for  ofif-road  opera¬ 
tion,  there  is  the  terrain  as  it  is  found  in  a  natural  state.  The  variability  in  composition  and 
consistency  of  natural  terrain  is  vast,  and  the  mechamcs  that  control  the  propulsion  system 
of  tires  in  this  medium  are  immense  in  complexity.  However,  as  with  many  other  complex 
engineering  problems,  useful  knowledge  and  understanding  can  be  gained  from  a  relatively 
simple  evaluation  of  the  system. 

This  research  focuses  on  the  interaction  of  tires  in  soil.  A  simple  steady-state  analysis 
of  a  tire  interacting  with  soil  can  be  used  to  evaluate  the  traction  performance  of  wheeled 
vehicles.  Figure  7  shows  a  powered  tire  with  the  loading  forces,  which  are  transferred  from 
the  axle  to  center  of  the  wheel,  and  the  reaction  stresses,  which  are  transferred  from  the  soil 
to  the  tire  along  the  interface.  As  shown  in  the  figure,  the  reaction  stresses  can  be  aggregated 
into  meaningful  horizontal  and  vertical  forces.  All  of  the  vertical  stresses,  normal  (o)  and 
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shearing  (t),  are  converted  into  a  normal  reaction  force  (N),  all  of  the  horizontal  normal 
stresses  are  converted  into  a  motion  resistance  force  (R),  and  all  of  the  horizontal  shearing 
stresses  are  converted  into  a  tractive  force  (T). 


Figure  7.  Simple  stress  analysis  for  tire/soil  interaction 


Figure  8  shows  the  same  tire  but  with  the  reaction  forces  applied  at  the  interface 
instead  of  the  stresses.  The  reaction  forces  are  located  at  unique  points  along  the  interface, 
and  these  points  are  designated  by  the  force  eccentricities.  A  steady-state  analysis  of  the 
horizontal  and  vertical  forces  acting  on  the  tire  demonstrates  that  the  normal  force  is  equal 
to  the  weight  on  the  wheel  (W^)  and  that  the  tractive  force  is  equal  to  the  sum  of  the  motion 
resistance  force  and  the  drawbar  force  (D).  This  simple  force  analysis  forms  the  fundamental 
logic  used  in  NRMM  by  which  drawbar  and  motion  resistance  relationships  are  used  to 
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determine  traction  potential.  The  simple  force  analysis  was  presented  herein  to  clarify  why 
drawbar  and  VCIj  (i.e.,  zero  drawbar)  performance  are  major  considerations  in  traction 
modeling. 


Tire  Deflection  Effect 

Tire  deflection  is  a  control  variable  for  wheeled  vehicles.  The  magmtude  of  tire 
deflection  is  directly  related  to  wheel  loading  and  inflation  pressure.  Reduced  inflation  pres¬ 
sure  translates  into  increased  deflection  which  results  in  an  increased  area  of  contact  between 
the  tire  and  the  ground.  This  increased  contact  area  has  a  positive  impact  on  traction  per¬ 
formance,  especially  in  soft  soils,  but  increased  tire  deflection  also  has  a  negative  impact  on 
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vehicle  performance  issues  like  rollover  stability  and  tire  durability.  Therefore,  operating 
deflections  (i.e.,  tire  inflation  pressures)  are  recommended  based  on  a  consideration  of  the 
intended  operating  environment,  design  loading,  and  design  speed  such  that  an  acceptable 
compromise  is  achieved  between  the  positive  and  negative  impacts  on  vehicle  performance. 

The  principal  manner  in  which  deflection  affects  contact  area  is  through  an  increase 
in  the  length  of  the  contact  area.  The  width  is  affected  too,  but  the  increase  in  width  is  small 
compared  to  the  increase  in  length.  Also,  the  increase  in  width  occurs  not  at  the  tire  tread  but 
up  along  the  side  wall  of  the  tire  carcass.  Thus,  the  effect  of  tire  deflection  on  the  area  of 
contact  at  the  tire/soil  interface  is  almost  entirely  the  result  of  increased  contact  length. 

A  mathematical  expression  was  derived  for  estimating  tire  contact  lengths  and  the 
effect  that  deflection  has  on  those  lengths.  The  assumption  was  made  that  tire  contact  length 
is  analogous  to  the  length  of  a  long  chord  to  a  circle  as  illustrated  in  Figure  9.  The  figure 
shows  that  the  length  of  a  long  chord  to  a  circle  (LC)  can  be  calculated  based  on  the  diameter 
of  the  drcle  (0)  and  the  middle  ordinate  height  (M)  using  the  equation  LC  =  . 

The  long  chord  analogy  can  be  used  to  estimate  the  hard-surface  contact  length  (Ls)  by  using 
the  diameter  of  the  tire  for  0  and  the  hard-surface  deflection  for  M  as  shown  below. 

Zj  =  2^db  -  6^ 

The  long  chord  analogy  was  used  throughout  this  research  for  estimating  tire  contact  lengths 
which  included  projected  bearing  lengths  resulting  from  sinkage  in  soft  soils  in  addition  to 
hard-surface  contact  lengths. 
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Figure  9.  Long  chord  analogy  for  tire  contact  length  estimates 


The  basic  equation  can  be  rewritten  with  LC  and  M  expressed  as  fractions  of  the  circle 
diameter  as  shown  below. 


With  the  expression  in  this  form,  the  nature  of  variations  in  LC  due  to  variations  in  M  for  a 
given  diameter  can  be  observed  as  shown  in  Figure  10.  The  figure  shows  that  LC  increases 
rapidly  being  50  percent  of  the  diameter  when  M  is  only  0.067  times  the  diameter,  and  LC 
reaches  a  maximum  value  equal  to  the  diameter  when  M  is  0.5  times  the  diameter.  Figure  10 
also  illustrates  that  care  must  be  taken  when  applying  the  long  chord  analogy  if  M  has  the 
potential  to  increase  beyond  0.5  times  the  diameter.  This  realization  becomes  important  when 
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uang  the  long  chord  analogy  to  estimate  projected  bearing  lengths  for  tires  in  soft  soils  where 
large  sinkage  occurs. 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

M  (fraction  of  diameter) 


Figure  10.  Long  chord  versus  middle  ordinate 

Model  Development  Tools 

Over  one  hundred  statistical  models  were  developed  during  this  research.  Some  of 
the  models  involved  establishing  a  mathematical  equation  to  describe  the  relationship  between 
a  performance  metric  (VCIi  or  DAVjo)  and  one  or  more  clearly  defined  system  parameters. 
Most  of  the  models,  however,  involved  establishing  the  form  of  the  principal  system 
parameter  in  conjunction  with  establishing  the  mathematical  equation  that  would  describe  the 
relationship.  The  principal  system  parameters  were  mathematical  terms  that  included  some 
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combination  of  vehicle  characteristics  and/or  soil  strength.  Depending  on  their  origin,  the 
principal  system  parameters  may  or  may  not  have  been  dimensionally  correct.  Some  of  the 
principal  system  parameters  were  theoretically  founded,  some  were  theoretically  founded  with 
empirical  improvements,  and  others  were  purely  empirical  being  completely  defined  by  the 
experimental  data. 

Statistical  modeling  techniques  were  used  to  develop  the  mathematical  equations  and 
the  principal  system  parameters  necessary  to  describe  the  traction  relationships.  Least  squares 
linear  regression  was  typically  used  to  develop  the  mathematical  equations,  and  curvilinear 
transformations  were  used  when  necessary,  which  was  always  the  case  in  the  drawbar 
modeling.  Numerical  optimization  techniques  were  used  to  develop  the  principal  system 
parameters.  In  Some  cases,  numerical  optimization  techniques  were  also  used  to  define  the 
mathematical  equations  rather  than  using  linear  regression. 

The  numerical  optimization  techniques  involved  iterative  trial  and  error  routines. 
These  routines  systematically  changed  the  value  of  constants  (i.e.,  coefficients,  additives,  or 
exponents)  within  the  principal  system  parameter  in  directions  that  minimized  the  standard 
error  (SJ,  or  maximized  the  coefficient  of  determination  (R^),  for  the  mathematical  equation 
describing  the  relationship.  It  was  necessary  for  the  constants  to  be  initialized  with  reasonable 
estimates  first.  Then  the  iterative  trial  and  error  routines  would  converge  on  better  estimates 
for  the  constants  by  using  intelligent  mathematical  algorithms.  The  algorithms  that  were  used 
involved  forward  or  central  differencing  for  determining  estimates  of  partial  derivatives,  a 
quasi-Newton  or  conjugate  method  for  deciding  whether  to  increase  or  decrease  the  con¬ 
stants,  and  a  tangent  or  quadratic  approach  for  obtaining  the  next  estimates  for  the  constants. 


-36- 


When  linear  regression  was  used  to  develop  a  mathematical  equation  and  numerical 
optimization  was  used  to  define  the  principal  system  parameter,  the  two  techniques  were 
combined  in  a  larger  iterative  process.  Since  numerical  optimization  is  an  iterative  approach, 
the  linear  regression  algorithms  were  included  within  the  numerical  optimization  process. 
Each  time  the  numerical  optimization  algorithms  would  increment  the  values  of  the  constants 
within  the  principal  system  parameter,  the  linear  regression  algorithms  would  update  the 
constants  within  the  mathematical  equation.  This  combined  process  resulted  in  quicker 
convergence  to  the  optimum  values  of  the  empirical  constants  in  the  principal  system 
parameters.  Whether  this  combined  iterative  process  or  one  of  the  individual  statistical 
modeling  tools  (i.e.,  linear  regression  or  numerical  optimization)  was  used,  computerized 
spreadsheet  software  was  utilized  to  accomplish  the  complex  and  laborious  task. 

For  clarification,  a  fictitious  example  is  presented.  Let  there  be  a  performance  metric 
Y  for  which  a  means  of  prediction  is  desired.  If  Y  is  linearly  related  to  system  parameters 
and  X2,  a  least  squares  linear  regression  model,  as  shown  below,  could  be  used  to  develop 
a  mathematical  equation  that  would  predict  Y  as  a  function  of  X^  and  X^. 

r  =  p„  ^  +  p,  jr, 

The  task  of  the  linear  regression  would  be  to  determine  values  for  coefficients  ,  p  j ,  and  pj 
that  minimize  the  sum  of  the  squared  errors  between  the  measured  data  and  the  predictions 
from  the  mathematical  equation.  An  additional  level  of  complexity  is  added  if  X^  is  an 
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undefined  parameter  such  as  that  shown  below  where  A,  B,  and  C  are  system  variables  (such 
as  tire  diameter,  soil  strength,  etc.)  and  p^,  Cj ,  and  Cj  are  undefined  constants. 

+  c.B 

X.  =  - i— 

*  C  +  Cj 

For  this  second  case,  numerical  optimization  could  be  used  in  combination  with  linear 
regression  as  described  in  the  preceding  paragraph.  The  task  of  the  numerical  optimization 
would  be  to  determine  values  of  Py,  Cj ,  and  Cj  that  produce  a  stronger  linear  relationship 
between  Xy  and  Y,  and  the  linear  regression  would  continuously  update  the  coeflScients  p^, 
Pj ,  and  Pj  during  the  iterative  process.  The  numerical  optimization  process  could  also  be 
used  alone  to  determine  all  sax  of  the  empirical  constants  Pq  ,  p  j ,  ^2^  P\,  ^y,  Cj ,  but  this 
approach  was  only  used  when  linear  regression  was  not  applicable. 

As  stated,  curvdlinear  transformations  were  always  necessary  in  the  drawbar  modeling. 
They  were  necessary  because  trends  observed  in  the  drawbar  relationships  were  not  linear. 
Instead  the  trends  were  curvilinear  with  forms  that  were  described  well  by  a  rectangular 


hyperbola.  The  general  form  of  the  rectangular  hyperbola  equation  is  shown  below. 


Pi 


In  order  to  determine  the  empirical  constants  P^,  pj,  and  Pj,  the  curvilinear 


transformation  shown  below  was  applied. 


X=  by, 
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Least  squares  linear  regresaon  was  used  to  determine  b^,  and  b^ ,  and  then  estimates  for 
pQ,  Pj,  and  P2  were  determined  using  the  equations  below  to  translate  the  regression 
constants  from  the  transformed  space  into  regression  constants  for  the  true  space. 

Po  =  *1  .  Pi=*2+*o^  ’  P2  = 

This  transformation  process  typically  produces  good  estimates  for  pj, ,  P  j ,  and  pj .  However, 
it  mininuzes  the  sum  of  the  squared  residuals  in  the  transformed  space,  not  in  the  true  space. 
Therefore  these  estimates  were  used  as  initial  estimates  for  numerical  optimization,  and 
numerical  optimization  was  used  to  determine  new  estimates  for  Pq,  Pj  ,  and  Pj  that  mini¬ 
mized  the  sum  of  the  squared  residuals  in  the  true  space  by  minimizing  or  maximizing  R^. 

One  other  noteworthy  item  specific  to  the  drawbar  modeling  is  the  manner  in  which 
iS'gWas  calculated.  is  largely  determined  by  the  magnitude  of  the  residuals.  Since  the 
residuals  are  equal  to  the  difference  between  the  measured  and  predicted  values,  the  predic¬ 
tions  will  obviously  influence  S^.  The  rectangular  hyperbola  equations  were  developed  to 
describe  the  drawbar  relationships,  but  they  were  not  used  to  predict  the  drawbar  perform¬ 
ance  metric  (Y)  for  all  values  of  the  system  parameter  (X).  The  equations  were  only  used  to 
predict  Y  for  observations  that  had  X  values  greater  than  or  equal  to  the  X-intercept  of  the 
equation.  For  observations  that  had  X  values  less  than  the  X-intercept  of  the  equation,  a 
prediction  of  zero  was  used. 

This  approach  of  determining  the  predictions  to  be  used  in  calculating  was  used 
for  the  following  reasons.  First,  there  physically  is  no  such  thing  as  negative  drawbar  for  a 
self-propelled,  independent  vehicle  that  is  attempting  to  produce  forward  motion  on  level 
terrain.  Passive  motion-resisting  influences  such  as  the  soil  do  not  push  any  harder  on  the 
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vehicle  than  the  vehicle  pushes  on  them.  Therefore,  if  the  maximum  tractive  force  that  the 
vehicle  can  generate  is  less  than  the  maximum  motion  resistance  force  that  the  soil,  etc.  can 
generate,  then  the  drawbar  force  (i.e.,  tractive  force  less  motion  resistance  force)  would  be 
zero.  Another  reason  arises  from  the  nature  of  the  rectangular  hyperbola  equation.  The 
rectangular  hyperbola  equation  has  a  vertical  asymptote  that  is  less  than  and  usually  very  close 
to  the  X-intercept  of  the  equation.  The  potential  exists  for  some  observations  to  have  X 
values  less  than  the  value  of  the  asymptote,  and  this  situatiori  usually  produces  extemely 
unrealistic  predictions  and  hence  greatly  misleading  estimates  for  S^. 


CHAPTER  IV 


veil  MODELING 

Overview 

The  veil  modeling  consisted  of  two  major  areas  of  focus.  The  first  focus  was  on 
evaluating  existing  relationships  or  parameters  that  use  cone  penetrometer-based  soil 
measurements  and  relatively  simple  vehicle  characteristics  to  determine  traction  potential  in 
soft  soils.  Four  existing  models  were  investigated  which  included  the  current  NRMM  rela¬ 
tionship  and  three  versions  of  Mean  Maximum  Pressure  (MMP).  The  second  focus  was  on 
the  development  of  new  relationships.  Seven  major  types  of  new  models  were  developed 
which  included  six  different  types  of  contact  pressure  models  and  a  numeric-type  model. 

For  most  of  the  major  types  of  new  models  that  were  developed,  more  than  one 
variant  was  investigated.  The  term  variant  is  used  here  to  describe  a  slight  variation  on  the 
same  basic  modeling  theme.  Slight  variations  include  such  things  as:  (1)  using  an  empirical 
constant  times  hard-surface  deflection  to  account  for  the  possibility  that  tires  deflect  differ¬ 
ently  in  soft  soils  than  on  a  hard  surface  and  (2)  adding  a  tire  construction  variable  as  a  secon¬ 
dary  parameter  to  account  for  the  possibility  that  radial  and  bias  ply  tires  perform  differently, 
even  with  similar  size,  loading,  and  deflection.  Fifty-six  variants  of  the  contact  pressure 
models  and  nine  variants  of  the  numeric-type  model  were  developed.  Additionally,  five  model 
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variants  based  on  the  standard  NRMM  relationship  were  investigated,  the  primary  purpose 
being  to  optimize  the  current  approach  to  the  VCIj  database. 

For  the  sake  of  clarity  and  concision,  all  of  the  individual  model  variants  are  not 
addressed  in  the  main  text  of  this  document.  Only  the  four  existing  models  and  the  best 
variants  of  the  new  development  models  are  described.  However,  a  listing  of  all  the  model 
variants  explored  and  a  table  summarizing  their  quality  are  provided  in  Appendix  B. 


Existing  Models 

The  standard  NRMM  relationship  is  based  primarily  on  a  parameter  called  the  mobility 
index  (MI).  MI  was  initially  developed  based  on  the  results  of  trafBcability  studies  conducted 
by  WES  from  1945  through  1951  [1951],  and  was  then  upgraded  to  its  current  form  based 
on  additional  traflficability  studies  conducted  during  the  1960s  [Kennedy  and  Rush  1968].  MI 
was  developed  to  predict  VCI50,  but  it  was  also  adopted  later  as  a  parameter  for  predicting 
veil.  The  MI  parameter  is  actually  an  aggregation  of  eight  smaller  factors,  each  of  which 
is  based  on  some  vehicle  characteristic  considered  to  influence  vehicle  mobility  in  soft  soils. 
MI  for  wheeled  vehicles  is  calculated  as  shown  below. 


MI  =  [  (^)  +  WLF  -  Cf]  (EF)  (7F) 

\{TEF)(.GF) 

The  eight  smaller  factors  are  calculated  as  follows  (see  nomenclature  for  descriptions). 


CPF 


w 
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^  =  c^i  (WIOOO)  + 

w<  2,000  Ibf  “»  Cyypj  =  0.553  and  Cj^  =  0 

2,000  ^  w<  13,500 -♦  c^f^J  =  0.033  and  =  1.050 

where  S 

13,500  <20,000 /*/  -»  c^j  =  0.142  and  =  -0.420 

20,000  =  0.278  and  =  -3.115 

GF  =  1  +  0.05  Cqp  where  c^p.  =  \  if  tire  chains  are  used  or  0  if  not 

EF  =  1  +  0.05  Cpp  where  c^p  =  I  if  PWR  <  10  hpiton  or  0  if  not 

TF  =  1  +  0.05  Cjp  where  Cjp  =  \  if  manual  transmission  or  0  if  automatic 


The  data  used  during  the  development  of  MI  involved  tire  deflections  around  15 
percent  of  the  section  height  only;  therefore,  deflection  was  not  considered  in  MI.  When  data 
were  later  obtained  for  tires  with  much  higher  deflections,  the  performance  of  MI  as  a 
predictor  for  VCIi  greatly  diminished.  WES  developed  the  deflection  correction  factor 
(DCF)  to  take  the  effect  of  tire  deflection  into  account.  DCF  is  calculated  as  shown  below. 


DCF 


I  0.15 1“’ 

I  8*/ 


DCF  actually  acts  to  normalize  VCIj  to  a  16  percent  deflection  magnitude  and  is  applied  as 
a  multiplier  to  equations  that  were  originally  developed  to  predict  VCIj  as  a  function  of  MI 
alone.  The  standard  NRMM  equations  for  predicting  VCIj  as  a  function  of  MI  and  DCF  are 
shown  below  [Ahlvin  and  Haley  1992]. 


VCf  =  f(MI,DCF) 
M<  115 

where  < 

M>n5 


VCL  =  (  1 1.48  +  0.2  M  -  — I  DCF 

^  (  MI +  2.14} 

VCf  =  [4.1  MI  DCF 
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Figure  11.  Standard  NRMM  relationship  for  VCIj 

MI  and  DCF  were  calculated  for  each  of  the  observations  in  the  VCIi  database,  and 
the  data  were  plotted  with  the  standard  NRMM  equations  on  a  x-y  plot  as  shown  in  Figure 
1 1 .  Note  that  it  was  necessary  to  divide  VCIj  by  DCF  in  order  to  generate  the  simple  x-y 
plot.  The  figure  shows  that  the  data  are  described  reasonably  well  by  the  standard  NRMM 
relationship,  but  the  relationship  appears  to  be  very  unstable  beyond  an  MI  of  about  100. 

The  standard  NRMM  relationship  was  developed  using  only  part  of  the  data  contained 
in  the  VCIj  database.  New  relationships  developed  using  all  of  the  data  would  have  an 
immediate  advantage  over  the  NRMM  relationship  since  their  empirical  constants  would  be 
optimized  to  the  data.  This  situation  meant  that  fair  comparisons  between  the  standard 
NRMM  relationship  and  any  of  the  other  relationships  evaluated  during  this  research  could 
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only  be  made  if  the  NRMM  relationship  was  optimized  to  the  complete  data  set.  Therefore, 
a  modified  version  of  the  NRMM  relationship  was  developed  that  optimized  the  empirical 
constants  in  MI,  DCF,  and  the  prediction  equations  (i.e.,  curve-fits). 


Figure  12.  Apparent  origin  of  the  Weight  Factor 

Optimizing  the  NRMM  relationship  was  a  difficult  task  because  of  the  large  number 
of  empirical  constants  and  because  the  Weight  Factor  (WF),  which  is  one  of  the  smaller 
factors  within  MI,  is  a  series  of  four  curve-fits  whose  origin  was  not  clear.  A  thorough 
dissection  of  MI  using  the  measured  data  revealed  a  probable  source  for  the  origin  of  WF. 
Figure  12  shows  a  plot  of  CPF  divided  by  10  versus  the  average  axle  weight  for  all  of  the 
veil  observations.  The  four  WF  curve-fits  are  also  included  in  the  plot,  and  the  VCIj  data. 
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wath  the  exception  of  the  circled  outliers,  follow  trends  very  similar  to  the  WF  equations. 
This  figure  indicates  that  the  linear  equations  defining  WF  were  developed  from  trends  in  CPF 
versus  average  axle  weight.  Therefore,  WF  was  optimized  by  using  linear  regression  to 
develop  new  equations  based  on  the  observed  trends  in  CPF/10  versus  w  (note  that  the  circled 
outliers  were  ignored).  These  improved  WF  curves  are  shown  in  Figure  13. 


O  1-Pass  VCI  Data  — —  Optimized  Weight  Factor  Curves 

Figure  13.  Optimized  Weight  Factor 

Most  of  the  remaining  MI  factors  were  also  optimized.  None  of  the  observations 
were  influenced  by  GF  or  EF  since  none  had  tire  chans  and  all  had  horsepower-per-ton  ratios 
of  10  or  greater,  and  CPF  is  a  theoretical  parameter  where  Q.Sd  represents  tire  radius 
[Kennedy  and  Rush  1968].  Therefore,  these  three  factors  were  not  optimized.  DCF  and  the 
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equations  for  predicting  VCIi  were  optimized  as  well  so  that  the  final,  improved  NRMM 
relationship  was  fully  optimized  to  the  VCIj  database.  Numerical  optimization  techniques 
were  used  for  the  entire  optimization  process  with  the  exception  of  the  WF  portion.  The 


optimization  changes  made  in  the  NRMM  relationship  are  as  follows. 


MI'  = 

(CPF)(WF')  ^  _  cp\ 

(£F)  (7F0 

DCF'  = 

^  0.15^ 

(TEF')(GF)  J 

.  > 

TEF'  = 

8.3  +  1.13  i  WLF'  ~ 

w 

tl 

100 

2000  ’ 

6.9 

WF'  =  (w/1000)  + 

w<  2,338 /A/ 

^WFl 

=  0.467 

and 

^WF2  ~ 

0 

where  ' 

2,33%  ^w <12,564 Ibf 

^WFl 

=  0.044 

and 

^WF2 

0.988 

12,564  w 

^WFl 

=  0.093 

and 

^WF2 

0.372 

TF'  =  1  +  0.10Cjy.  where  Cjp  =  ^  if  manual  transmission  or  0  if  automatic 
VCf  =  f(M\DCF') 


where 


M'  ^  115 
MI'  >  115 


VCf 

VCf 


10.2  +  ^lllMl'  -  — — — 
Ml'  +  2.62  ) 

3.216  MI'  DCF' 


DCF' 


MI'  and  DCF'  were  calculated  for  each  of  the  observations  in  the  VCIj  database, 


and  the  data  were  plotted  with  the  optimized  NRMM  equations  on  a  x-y  plot  as  shown  in 
Figure  14.  Note  that  it  was  necessary  to  divide  VCIj  by  DCF'  in  order  to  generate  the 
simple  x-y  plot.  The  figure  shows  that  the  data  are  described  reasonably  well  by  the  opti¬ 
mized  NRMM  relationship.  There  is  little  improvement  in  the  \dsual  appearance  of  the 
optimized  relation  over  the  standard.  The  most  noticeable  change  is  that  the  break  in  slope 
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between  the  two  curve-fits  at  a  mobility  index  of  1 15  is  not  nearly  as  drastic.  It  ■mil  be  shown 
in  the  last  section  of  this  chapter  that  significant  improvement  in  quality  was  obtained. 


Figure  14.  Optimized  NRMM  relationship  for  VCIi 

The  other  three  existing  models  explored  for  VCIj  are  all  versions  of  the  Mean 
Maximum  Pressure  (MMP)  criterion  for  evaluating  soft  soil  performance  potential.  Where 
MI  is  the  soft  soil  performance  indicator  developed  by  and  for  the  United  States  military, 
MMP  is  the  soft  soil  performance  indicator  developed  by  and  for  the  United  Kingdom 
military.  The  original  MMP  for  wheeled  vehicles  was  developed  in  the  mid  1970s  by  D. 
Rowland  as  a  parameter  for  estimating  mean  ground  pressures  under  tires  that  would  be 
directly  comparable  to  the  MMP  he  developed  a  few  years  earlier  for  tracked  vehicles 
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[Larminie  1992],  Rowland  developed  MMP  for  wheels  based  in  large  part  on  the  results  of 
the  Numeric  Approach  research  performed  by  WES.  MMP  was  actually  a  recast  and  slightly 
altered  form  of  Freitag’s  Numeric  for  tires  in  clay  soils. 

MMP  is  used  more  as  an  indicator  of  vehicle  performance  potential.  The  trafficability 
research  conducted  by  WES,  which  led  to  the  development  of  the  standard  NRMM 
relationship,  focused  on  developing  a  system  of  measuring  soft-soil  performance,  and  this 
culminated  in  the  development  of  the  VCI  performance  metric  (i.e.,  minimum  soil  strength 
necessary  to  support  the  vehicle).  The  MI  and  DCF  prediction  parameters  were  developed 
for  the  sole  purpose  of  predicting  VCI.  VCI  was  used  as  the  basis  for  comparing  (and 
spedfying)  the  soft-soil  performance  since  it  allows  direct  comparison  of  different  vehicles, 
regardless  of  the  type  of  running-gear.  MMP  on  the  other  hand,  is  essentially  used  as  the 
performance  metric  where  higher  MMP  values  imply  worse  performance.  Correlations  appar¬ 
ently  have  been  made  between  MMP  and  soft-soil  performance  data,  but  no  sound 
relationship  has  ever  been  established  to  use  MMP  as  a  predictor  for  VCIj.  Therefore, 
relationships  had  to  be  developed  to  evaluate  the  three  versions  of  MMP. 

The  first  MMP  model  evaluated  was  based  on  Rowland’s  original  version  of  MMP 
which  was  calculated  as  shown  below. 


MMPj^ 


0.91  KGVW 
n  m  ^b/h 


where  for  all-drive  vehicles: 


m  =  2  ^K  =  3.66 
m  =  3  K  =  3.90 
/w  =  4  K  =  4.10 
ot  =  5  =>  K  =  4.32 
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This  formula  is  actually  slightly  different  from  the  MMP  formula  found  in  the  referenced 
publication.  The  true  formula  had  the  constant  2  in  place  of  the  average  number  of  tires  per 
axle  (n),  but  two  observations  in  the  VCIj  database  that  had  dual  tires  on  the  rear  (n=3.333) 
demonstrated  that  n  was  more  appropriate.  The  VCIi  data  were  plotted  versus  MMPr,  and 
a  linear  trend  was  apparent.  Therefore,  linear  regression  was  used  to  develop  a  mathematical 
equation  for  predicting  VCIi  as  a  function  of  MMPr.  The  resulting  equation  is  shown  below. 
FC/j  =  3.65  +  0  All  MMP 

The  equation  is  plotted  with  the  measured  data  in  Figure  15.  The  figure  shows  that  the  trend 
is  very  linear  and  that  MMPr  describes  the  VCIj  performance  data  well. 


Figure  15.  Relationship  between  Rowland’s  MMP  and  the  VCIi  data 
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The  second  MMP  model  was  based  on  a  slightly  modified  version  of  MMP  that  was 
proposed  by  Larminie  [1992],  Larminie  replaced  the  tire  section  height  with  diameter  in  the 
deflection  ratio,  and  this  change  was  proposed  based  on  an  evaluation  of  the  effect  that 
deflection  has  on  tire  contact  length.  He  also  adjusted  the  factor  K  to  keep  the  modified 
MMP  calculations  in  the  same  range  of  magnitudes  as  the  original  MMP  formula  since  MMP 
is  used  directly  for  comparing  vehicles.  The  MMP  values  for  this  model  were  calculated  as 
shown  below,  where  once  again  n  was  used  instead  of  the  constant  2  as  described  in  the 
previous  paragraph. 

_a97£:CT2L. 

where  for  all-drive  vehicles: 

m  =  2~*  K'=\.Z3 
/w  =  3  -  1.95 

m  =  A  K’  =  2.05 
m  =  5=*K'  =  2.\6 

The  veil  data  were  plotted  versus  MMPl,  and  once  again  a  linear  trend  was  apparent. 
Therefore,  linear  regression  was  used  to  develop  a  mathematical  equation  for  predicting  VCIi 
as  a  function  of  MMPl.  Th®  resulting  equation  is  shown  below. 

FC7i  =  3.68  +  OAZQMMP^ 

The  equation  is  plotted  with  the  measured  data  in  Figure  16.  The  figure  shows  that  the  trend 
is  very  linear  and  that  MMPl  describes  the  VCIj  performance  data  well.  The  figure  also 
demonstrates  that  MMPl  describes  the  data  in  a  fashion  nearly  identical  to  MMPr. 
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Figure  16.  Relationship  between  Larminie’s  MMP  and  the  VCIj  data 


The  third  version  of  MMP  was  based  on  the  latest  proposed  version  reported  by 
Maclaurin  [1997].  Maclaurin  proposed  significant  changes  in  the  form  and  use  of  MMP.  He 
agreed  with  Larminie’s  recommendation  that  the  diameter  should  be  used  instead  of  section 
hdght  in  the  deflection  ratio,  but  he  chose  to  mathematically  combine  the  d  from  the  6/d  ratio 
with  the  other  d  in  MMP  such  that  the  deflection  ratio  is  no  longer  shown.  He  did  not 
recommend  any  changes  be  made  to  the  K  factor,  instead  recommending  it  be  removed  from 
consideration.  The  proposed  new  MMP  formula  was  derived  from  the  results  of  recent 
drawbar  experiments  conducted  by  the  Defense  Evaluation  and  Research  Agency  (DERA) 
with  a  field  deployable,  computer  controlled,  single- wheel  tire  tester.  The  MMP  values  for 
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this  model  were  calculated  as  shown  below,  where  like  the  previous  two  versions,  n  was  used 
instead  of  the  constant  2. 

GVW 

MMP^  = - 5^111 - 

The  veil  data  were  plotted  versus  MNdP^,  and  once  again  a  linear  trend  was  apparent. 
Therefore,  linear  regression  was  used  to  develop  a  mathematical  equation  for  predicting  VCIi 
as  a  function  of  MNIP^.  The  resulting  equation  is  shown  below. 

FC/i  =  2.53  +  l35MMPj^ 

The  equation  is  plotted  with  the  measured  data  in  Figure  17.  The  figure  shows  that  the  trend 
is  very  linear  and  that  MMP^  describes  the  VCIi  performance  data  reasonably  well. 


Figure  17.  Relationship  between  Maclaurin’s  MMP  and  the  VCIj  data 
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As  mentioned  earlier,  the  ori^nal  MMP  was  actually  a  recast  and  slightly  altered  form 
of  Freitag’s Numeric  for  tires  in  clay  soils.  Therefore,  Freitag’s  Numeric  was  recast 

in  a  VCI  form  and  evaluated  as  a  predictor  for  VCIi.  This  model  was  evaluated  primarily  to 
see  how  it  compared  with  the  three  MMP  models.  Freitag’s  Numeric  is  calculated  as  shown. 

RCInbdfbY'^ 

aJ 

Recasting  the  numeric  in  a  VCI  form  simply  involved  removing  the  soil  strength  (i.e.,  RCI), 
since  predicting  a  soil  strength  was  the  new  goal,  and  inverting  the  formula  so  that  it  would 
be  in  units  of  pressure.  This  was  the  same  approach  Rowland  took  in  developing  the  original 
MMP.  The  recast  form  is  shown  below. 

^ 

nbd^ 

The  veil  data  were  plotted  versus  and,  like  the  MMP  models,  the  relationship  was 
linear.  Therefore,  linear  regression  was  used  to  develop  a  mathematical  equation  for 

predicting  VCIj,  and  the  resulting  equation  is  shown  below. 

FC/j  =  2.99  +  1.5^ 

The  equation  is  plotted  with  the  measured  data  in  Figure  1 8.  The  figure  shows  that  the  trend 
is  very  linear  and,  just  as  expected,  Freitag’s  Numeric  describes  the  VCIi  data  well. 
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Figure  18.  Relationship  between  Freitag’s  Numeric  and  the  VCIj  data 

New  Development  Models 

As  stated,  seven  major  types  of  new  development  models  were  explored.  These 
included  six  different  types  of  contact  pressure  models  and  a  numeric-type  model.  Only  two 
of  the  contact  pressure  models  and  the  numeric-type  model  will  be  discussed,  and  only  the 
best  variants  of  those  three  will  be  presented.  Before  the  individual  models  are  addressed, 
clarification  is  provided  on  what  is  meant  by  ‘contact  pressure  model’  and  ‘numeric-type 
model’  and  explanations  are  ^ven  as  to  why  these  types  of  models  were  the  focus  in  the  VCIj 
modeling. 

Contact  pressure  models  are  founded  on  the  theory  that  bearing  capacity  is  the 
controlling  factor  in  soil  response  at  VCIj  and  that  the  contact  pressure  imposed  by  the 
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vehicle  on  the  soil  drives  the  response.  These  models  use  logical,  though  possibly  simplistic, 
derivations  of  average  contact  area  (and  hence  pressure)  at  the  vehicle/terrain  interface. 
Different  ideas  and  levels  of  complexity  can  be  used  in  estimating  average  contact  area,  and 
hence  several  types  of  contact  pressure  models  were  explored.  In  actuality,  the  parameters 
developed  for  the  contact  pressure  models  during  this  research  do  not  attempt  to  calculate 
the  true  contact  area  for  tires  in  soil.  Instead,  they  attempt  to  cast  the  principal  tire 
dimensions  into  mathematical  forms  that  give  logical  consideration  to  the  influence  of  the  tire 
dimensions  on  the  true  contact  area.  This  approach  should  result  in  parameters  that  are  highly 
correlated  with  the  true  contact  area  and,  if  the  theory  is  correct,  VCIj  performance. 

Numeric-type  models  are  purely  empirical  with  no  pre-consideration  of  the  controlling 
mechanism.  Rather  than  using  a  logical  derivation  from  a  consideration  of  the  influence  that 
tire  dimensions  have  on  contact  area,  th^^  allow  the  data  to  completely  define  the  form  of  the 
prindpal  parameter.  They  are  loosely  founded  on  the  ideas  of  dimensional  analysis  and  simili¬ 
tude.  The  approach  with  numeric-type  models  is  to  develop  a  single,  simple  parameter  that 
casts  the  principal  tire  dimensions  into  a  mathematical  form  which  forces  the  units  to  work 
out  correctly  while  produdng  the  strongest  possible  correlation  with  the  performance  metric. 

The  focus  for  the  new  development  models  was  placed  on  contact  pressure  and 
numeric-type  models  partly  because  of  the  potential  that  was  demonstrated  by  the  existing 
MMP  formulas.  The  existing  MMP  formulas  were  developed  using  numeric-type  models,  and 

thdr  intent  was  to  predict  contact  pressure.  Another  reason  these  types  of  models  were  the 

\ 

focus  arose  fi’om  a  consideration  of  what  VCIj  physically  represents.  VCIj  is  a  zone  of  soil 
consistency  where  typical  vehicles  incur  high  sinkage  and  verge  on  drag^g  the 
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undercarriage,  which  makes  it  natural  to  infer  that  bearing  capacity  is  the  controlling  factor 
in  the  soil  response,  greatly  outweighing  traction.  It  is  probable  that  every  researcher  who 
has  ever  delved  into  trafficability  research  has  intuitively  thought  that  contact  pressure  was 
the  dominant  factor  from  the  vehicle  (i.e.,  loading)  side  of  the  VCIj  problem. 

The  first  new  development  model  is  a  contact  pressure  model.  It  uses  an  average 
contact  pressure  parameter  referred  to  as  CPI  for  the  principal  modeling  parameter.  CPI  is 
defined  in  general  terms  based  on  the  width  (  )  and  length  )  of  the  effective  contact 

area  under  the  tires  as  shown  below. 


The  basic  logic  behind  CPI  is  illustrated  graphically  in  Figure  19.  It  was  hypothesized  that 
the  effective  contact  width  remains  essentially  constant  and  can  be  represented  by  the  nominal 
section  width  of  the  tire.  It  was  also  hypothesized,  as  illustrated  in  Figure  19a,  that  each  of 
the  axles  has  an  effective  contact  length  ( Zj )  under  the  tires  that  originates  from  hard-surface 
deflection  and  that  the  first  axle,  which  is  incurring  the  majority  of  the  motion  resistance,  has 
additional  contact  length  (Z^)  along  the  front  edge  of  the  tires  due  to  sinkage.  These  hypo¬ 
theses  resulted  in  CPI  being  defined  as  shown  below. 

=  b  cmd  Z^  =  Zg  +  VmL^ 

_ w _ 

nb(L,^VmL,) 


CPI  = 


Figure  19.  Basic  logic  behind  the  average  contact  pressure  parameter  CPI 
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Figure  19b  shows  the  effective  geometry  used  to  mathematically  describe  the  two  contact 
lengths  Z-g  and  L^.  ^plying  the  long  chord  analogy  resulted  in  the  two  contact  lengths  being 
defined  as  shown  below. 


LC  = 

Zg  =  ZCg  =►  Zg  =  2)/db  - 


4  = 


ZC,-Zg 


L  =  VzLC  -yjdb-b^ 


where  i 


z<  —  -  b  =>  ZC 
2  * 

z  -  b  LC. 
2 


2sldib  +z)  -  (b  +zf 
d 


Once  the  basic  mathematical  form  of  the  CPI  parameter  had  been  established,  a  x-y 
plot  of  veil  versus  CPI  was  generated  using  the  observations  from  the  VCIj  database,  where 
clearance  height  (h^)  was  used  as  the  estimate  for  z.  The  plot  showed  that  a  linear  relation¬ 
ship  existed  and  that  CPI  had  good  potential  for  describing  VCIi.  Then  an  intensive  effort 
was  conducted  in  which  empirical  improvements  for  CPI  were  explored  through  the 
development  of  model  variants.  Many  different  ideas  were  explored  in  the  model  variants 
which  included  using  empirical  multipliers  within  the  CPI  parameter  to  modify  the  influence 
of  the  indivadual  parts  of  the  estimate,  adding  secondary  parameters  to  the  linear  rela¬ 
tionship  between  VCIi  and  CPI,  using  other  estimates  (besides  h^)  for  z,  and  others.  A  major 
source  of  ideas  for  the  model  variants  came  from  the  existing  VCIi  models  described  in  the 
last  section.  This  effort  produced  significant  improvements  to  the  relationship. 
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The  best  model  variant  used  as  the  estimate  for  z,  had  two  empirical  modifications 
within  the  CPI  parameter,  and  included  a  secondary  parameter  in  the  relationship.  The  first 
modification  reduced  the  influence  of  by  multiplying  a  constant  (Cj )  times  the  ViLC^  part 
of  the  calculation.  The  second  modification  involved  an  adjustment  to  the  effective 
contact  area  by  multiplying  aspect  ratio  raised  to  a  power  (/7j)  in  the  denominator  of  the  CPI 
formula.  The  secondary  parameter  included  in  the  relationship  was  a  simple  tire  construction 
factor  (TCP)  that  accounted  for  the  possibility  of  a  difference  in  performance  between  radial 
and  bias  ply  tires.  The  form  of  the  best  CPI  model  variant  is  shown  below. 


VCI,  =  Po 


TCP  =  where  Cj^  =  0  //  radial  or  1  if  bias  ply 

The  combined  linear  regression  and  numerical  optimization  process  was  used  to  determine 
the  empirical  constants  in  the  relationship.  The  process  resulted  in  Cj  being  0.74  and 
being  -0.47.  The  resulting  final  equations  are  shown  below  where  Cj  has  been  combined  vrith 
the  Vi  in  ViLC^  and  the  aspect  ratio  has  been  inverted  for  a  positive  exponent. 


VCf  =  2.88  +  1.78CP7;9  +  1.58  TCF 
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CPl  (5  and  TCP  were  calculated  for  each  of  the  observations  in  the  VCIj  database, 
and  the  data  were  plotted  with  the  equation  describing  the  relationship  on  a  x-y  plot  as  shown 
in  Figure  20.  Note  that  it  was  necessary  to  subtract  1.58  TCP  from  VCIj  in  order  to  generate 
the  simple  x-y  plot.  The  figure  shows  that  the  measured  VCIi  data  are  described  very  well 
by  the  best  CPI  relationship. 


Figure  20.  Best  CPI  relationship  for  VCIi 

The  second  new  development  model  is  also  a  contact  pressure  model.  It  uses  an 
average  contact  pressure  parameter  referred  to  as  CP2  for  the  principal  modeling  parameter. 
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CP2,  like  CPI,  is  defined  in  general  terms  based  on  the  width  and  length  (1^4)  of 
effective  contact  area  under  the  tires  as  shown  below. 


CP2  = 


w 


”  ^CA  ^CA 


The  basic  logic  behind  CP2  is  iUustrated  graphically  in  Figure  21 .  It  was  hypothesized,  once 
again,  that  the  effective  contact  width  remains  essentially  constant  and  can  be  represented  by 
the  nominal  section  width  of  the  tire.  It  was  also  hypothesized,  as  illustrated  in  Figure  21a, 
that  all  of  the  axles  have  an  effective  contact  length  under  the  tires  that  originates  fi-om 

hard-sur&ce  deflection  plus  an  empirical  correction  variable  (c).  These  hypotheses  resulted 
in  CP2  being  defined  as  shown  below. 

=  b  and 


CP2  = 


w 


nbL 


Figure  21b  shows  the  effective  geometry  used  to  mathematically  describe  the  contact  length 
L^+c-  Applying  the  long  chord  analogy  resulted  in  being  defined  as  shown  below. 


LC  =  2/0A/  - 


h,.  =  LC,,,  -  =  2^d(b*c)  -  (&*cf 


Once  the  basic  mathematical  form  of  the  CP2  parameter  had  been  established,  a 
couple  of  model  variants  were  developed.  The  best  variant  used  an  empirical  constant  (Cj ) 
times  tire  diameter  as  the  estimate  for  c  and  used  the  best  empirical  adjustments  that  were 
found  during  the  effort  to  empirically  improve  CPI.  The  empirical  adjustments  were:  (1) 
aspect  ratio  raised  to  a  power  (/?j)  which  was  used  as  an  adjustment  to  the  effective  contact 
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Figure  21.  Basic  logic  behind  the  average  contact  pressure  parameter  CP2 
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area  in  the  denominator  of  CP2  and  (2)  TCP  which  was  used  as  a  secondary  parameter  in  the 
relationship  between  VCIi  and  CP2  to  account  for  performance  differences  between  radial 
and  bias  ply  tires.  The  many  other  empirical  adjustments  that  were  evaluated  for  CPI  were 
not  re-evaluated  for  CP2  because  thdr  influence  on  CP2  would  have  been  very  similar  to  their 
influence  on  CPI.  The  form  of  the  final  CP2  model  variant  is  shown  below. 


VCI,  =  Po  +  PiCP2[  -H  p^rcF 


TCF  =  c^j,  where  =  0  //  radial  or  1  if  bias  ply 

The  combined  linear  regression  and  numerical  optimization  process  was  used  to  determine 
the  empirical  constants  in  the  relationship.  The  process  resulted  in  Cj  being  0.026  and 
being  -0.52.  The  resulting  final  equations  are  shown  below  where  the  aspect  ratio  has  been 
inverted  for  a  positive  exponent. 


FC/j  =  1.95  +  2.03  CP2(  +  1.95  TUF 

CP2[  and  TCF  were  calculated  for  each  of  the  observations  in  the  VCIj  database, 
and  the  data  were  plotted  with  the  equation  describing  the  relationship  on  a  x-y  plot  as  shown 
in  Figure  22.  Note  that  it  was  necessary  to  subtract  1 .95  TCF  fi-om  VCIi  in  order  to  generate 
the  simple  x-y  plot.  The  figure  shows  that  the  measured  VCIj  data  are  described  very  well 
by  the  best  CP2  relationship. 
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Figure  22.  Best  CP2  relationship  for  VCIj 

The  last  new  development  model  is  a  numeric-type  model.  As  previously  stated,  a 
numeric-type  model  allows  the  data  to  completely  define  the  form  of  a  simple,  principal 
parameter  (referred  to  as  a  numeric)  under  the  constraint  that  the  units  must  be  correct.  The 
veil  performance  metric  has  pressure  units  (i.e.,  force/length^).  Therefore,  the  numeric 
parameter  was  forced  to  also  have  pressure  units.  Several  approaches  were  used  for  develop¬ 
ing  the  numeric-type  models,  but  only  the  approach  used  in  developing  the  best  numeric 
variant  will  be  described. 

In  developing  the  best  VCIj  numeric  variant,  w/n  was  the  obwous  choice  for  including 
the  force  units.  The  approach  taken  to  include  the  length^  units  was  to  use  a  product  of  b, 
d,  and  three  dimensionless  tire  ratios  where  each  ratio  was  raised  to  a  power  p^,  and 
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Pj).  The  ratios  chosen  used  d  as  the  characteristic  length  in  dimensionless  ratios  with  each 
of  b,  h,  and  6  so  that  all  four  of  the  tire  shape  dimensions  could  be  included  in  the  final 
numeric.  Three  ratios  O-e.,  powers)  were  used  because  there  were  four  length  variables  less 
one  constraint.  This  modeling  approach  allowed  the  inclusion  of  b,  d,  h,  and  6  in  the  final 
numeric  with  each  rased  to  its  own  numerically  optimized  power  under  the  constraint  that 
the  summation  of  the  four  powers  must  equal  two. 

Once  the  basic  mathematical  form  of  the  numeric  was  decided,  a  statistical  model  was 
used  to  develop  the  relationship  between  VCIj  and  the  numeric.  A  linear  equation  was  used 
due  to  the  linear  behavior  demonstrated  in  the  MMP  models  and  the  two  average  contact 
pressure  models.  Also,  due  to  its  demonstrated  benefit  in  the  average  contact  pressure 
models,  TCF  was  used  as  a  secondary  parameter  in  the  relationship  to  account  for  perform¬ 
ance  differences  between  radial  and  bias  ply  tires.  The  form  of  the  final  numeric-type  model 
is  shown  below. 

VCI,  =  Po  ^2TCF 

N  _ _ ^ ^ _ 

nbd(bld)^'{dlhY^{bldY^ 

TCF  =  Cj^  where  =  0  if  radial  or  1  if  bias  ply 

The  combined  linear  regression  and  numerical  optintization  process  was  used  to  determine 
the  empirical  constants  in  the  relationship.  The  process  resulted  in  being  0.42,  being 
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0.66,  and  being  0.32.  The  resulting  final  equations  are  shown  below  where  the  numeric 
has  been  simplified  by  coUapsing  the  ratios  and  collecting  like  terms. 

^ 

0.92  50.32 

rC/j  =  1.56  +  2.02  +  2.03  TCF 

and  TCF  were  calculated  for  each  observation  in  the  VCIj  database,  and  the  data 
were  plotted  with  the  equation  describing  the  relationship  on  a  x-y  plot  as  shown  in  Figure 
23.  Note  that  it  was  necessary  to  subtract  2.03  TCF  from  VCIi  in  order  to  generate  the 
simple  x-y  plot.  The  figure  shows  that  the  measured  VCIi  data  are  described  very  well  by  the 
best  numeric  relationship. 


Figure  23.  Best  numeric  relationship  for  VCIj 
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Model  Comparisons 

Several  existing  and  new  VGIi  relationships  were  presented  in  the  preceding  sections 
of  this  chapter,  and  a  simple  x-y  plot  of  each  relationship  was  shown.  The  x-y  plots  allowed 
a  \dsual  interpretation  of  the  quality  with  which  each  relationship  described  the  measured  data, 
but  hi^er  accuracy  methods  of  assessing  quality  were  necessary  to  make  sound  comparisons 
between  the  relationships.  For  this  purpose,  the  standard  error  (5^)  and  adjusted  coefBcient 
of  determination  (Adj.  R^)  measures  of  quality  were  used.  These  quality  measures  were 
determined  during  the  development  of  each  of  the  many  model  variants,  and  they  were  used 
as  optimization  targets  in  models  that  required  numerical  optimization  as  previously  explained 
in  Chapter  HI.  Table  4  lists  these  measures  along  with  their  associated  degrees  of  freedom 
for  each  of  the  relationships. 

Several  comparisons  can  be  made  using  Table  4.  First,  the  table  demonstrates  that  the 
optimized  NRMM  relationship  described  the  data  significantly  better  (i.e.,  lower  and 
higher  R^)  than  the  standard  NRMM  relationship^.  The  table  also  demonstrates  that 
Rowland’s  original  version  of  MMP  described  the  data  better  than  both  Larminie  s  and 
Maclaurin’s  versions,  but  Frdtag’s  numeric  (fi-om  which  MMP  originated)  described  the  data 
better  than  any  of  the  three.  In  terms  of  the  new  relationships,  the  best  CP2  relationship 
described  the  data  slightly  better  than  the  best  CPI  relationship,  and  the  best  numeric 
relationship  ranked  a  close  third.  Most  significantly,  the  table  demonstrates  that  all  three  of 

*  -  The  degrees  of  fieedom  for  the  NRMM  relationships  were  not  the  same  for  all  of  the  observations 
twaiicft  the  wei^t  fector  equations  had  either  one  or  two  empirical  constants  and  because  the  two  equations 
describing  the  relationship  between  VCI,/DCF  and  MI  had  either  two  or  four  empiric^  constants.  This 
situation  resulted  in  there  being  four  possible  values  for  the  degrees  of  freedom  associated  with  each  individual 
observation. 
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the  new  relationships  describe  the  data  better  than  any  of  the  existing  relationships,  with  only 
the  newly  optimized  NRMM  relationship  coming  anywhere  close  to  the  quality  of  the  new 
relationships.  These  results  indicate  that  it  is  definitely  feasible  to  develop  new  VCIi  relation¬ 
ships  that  meet  the  dedsion  criteria  stated  in  the  objectives  for  this  research.  The  results  also 
indicate  that  contact  pressure  under  the  vehicle  is  the  driving  factor  in  VCIj  performance. 

Table  4.  Quality  comparison  for  the  VCIj  modeling  results 


Relationship 

Degrees  of 
Freedom 

Standard 
Error  (psi) 

Adjusted 

R-squared 

Standard  NRMM  (MI&DCF) 

70,  69, 

68,  or  67 

2.68 

90.3 

Optimized  NRMM 

70,  69, 

68,  or  67 

. i 

2.44 

91.9 

Rowland’s  MMP 

73 

2.91 

88.5 

Larminie’s  MMP 

73 

2.96 

88.1 

Maclaurin’s  MMP 

74 

3.01 

87.7 

Freitag’s  Numeric  in  a  VCI,  Form 

76 

2.82 

89.2 

Best  CPI  Variant 

74 

2.39 

92.2 

Best  CP2  Variant 

74 

2.37 

92.4 

Best  Numeric  Variant  for  VCI, 

73 

2.40 

92.2 

One  other  major  comparison  was  made  from  the  results  of  the  VCIi  modeling.  This 
comparison  involved  the  best  CP2  and  numeric  relationships.  It  was  noted  that  the  coeflB- 
cients  on  CP2[  and  both  equal  to  2.03  and  that  the  coefficient  on  TCP  in  both 
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relationships  was  essentially  equal  to  two.  This  observation  prompted  a  closer  look  at  the 
similarities  between  CP2  [  and  .  Table  5  shows  a  comparison  between  the  exponents  (or 
powers)  in  these  two  parameters,  and  it  indicates  that  the  two  parameters  are  very  similar  in 
form.  Figure  24  shows  a  x-y  plot  of  versus  CP2[,  and  it  indicates  that  the  two 
parameters  are  essentially  identical  in  terms  of  their  calculated  magmtudes  for  each  of  the 
observations  in  the  VCIj  database. 

Table  5.  Comparison  of  exponents  between  best  CP2  and  Numeric  relationships 


Exponents  in  Principal  Parameter  * 


Relationship 

d 

b 

h 

5 

Best  CP2  Variant 

1.00*’ 

1.52 

0.52 

0.50  *’ 

Best  Numeric  Variant  for  VCI, 

0.92 

1.42 

0.66 

0.32 

®  -  d,  b,  and  5  in  the  denominator;  h  in  the  numerator 
’’  -  approximate  values 


Two  major  conclusions  can  be  drawn  from  this  result.  First,  numeric-t5q)e  relation¬ 
ships  can  be  very  useful  in  areas  where  little  is  understood  about  the  controlling  mechanism. 
Second,  the  measured  VCIj  data  apparently  cast  into  a  form  that  is  nearly  identical  to 
CP2  [ .  Since  had  the  flexibility  to  mold  itself  into  the  best  possible  form  for  describing 
the  veil  data,  the  fact  that  it  molded  itself  into  a  form  nearly  identical  to  CP2  [  provides  even 
more  confidence  in  the  applicability  of  CP2  j  for  describing  VCIi  performance. 
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Figure  24.  Similarity  between  best  CP2  and  Numeric  parameters  for  VCIj 


CHAPTER  V 


DRAWBAR  MODELING 

Overview 

The  drawbar  modeling  consisted  of  two  major  areas  of  focus.  The  first  focus  was  on 
evaluating  existing  relationships  or  parameters  that  use  cone  penetrometer-based  soil 
measurements  and  relatively  simple  vehicle  characteristics  to  determine  drawbar  performance 
potential  in  clay  soils.  Five  existing  models  were  investigated  which  included  the  current 
NRMM  relationship  and  four  existing  drawbar  numerics.  The  second  focus  was  on  the 
development  of  new  relationships.  Four  major  types  of  new  models  were  developed  which 
included  a  bearing  ratio  model,  an  excess  bearing  ratio  model,  an  excess  RCI  ratio  model,  and 
a  numeric-type  model. 

For  all  of  the  major  types  of  new  models  that  were  developed,  more  than  one  variant 
was  investigated.  As  in  the  VCIj  modeling,  the  term  variant  is  used  here  to  describe  a  slight 
variation  on  the  same  basic  modeling  theme.  Slight  variations  in  the  drawbar  modeling 
included  such  things  as:  (1)  using  different  estimates  for  contact  pressure  in  the  bearing  ratio 
model  and  (2)  using  different  combinations  of  four  or  less  tire  shape  dimensions  in  the 
numeric-tjTJe  model.  Three  variants  of  the  bearing  ratio  model,  four  variants  of  the  excess 
bearing  ratio  model,  two  variants  of  the  excess  RCI  ratio  model,  and  sixteen  variants  of  the 
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numeric-type  model  were  developed.  Additionally,  two  model  variants  based  on  the  standard 
NRMM  relationship  were  investigated,  the  primary  purpose  being  to  optimize  the  current 
approach  to  the  field  drawbar  database. 

For  the  sake  of  clarity  and  concision,  all  of  the  individual  model  variants  are  not 
addressed.  Only  the  five  existing  models  and  the  best  variants  of  the  new  development 
models  are  compared  in  the  mmn  text  of  this  thesis.  However,  a  listing  of  all  the  model 
variants  explored  and  a  table  summarizing  their  quality  are  provided  in  Appendix  C. 

Existing  Models 

The  standard  NRMM  relationship  is  based  on  a  parameter  called  the  excess  soil 
strength  (RCIx).  RCIx  was  developed  based  on  observations  made  by  WES  researchers 
during  the  1960s  [Priddy  1995].  During  the  early  years  of  drawbar  measurements,  curves  of 
drawbar  versus  soil  strength  (in  RCI)  were  plotted  for  every  vehicle  configuration  tested. 
Once  several  vehicle  configurations  had  been  tested  over  a  fair  range  of  soil  strengths  (usually 
veil  to  about  150  or  200  RCI),  it  was  observed  that  drawbar  versus  RCI  curves  for  vehicles 
with  similar  characteristics  (e.g.,  wheeled  or  tracked)  were  very  similar  in  shape  and 
magnitude.  This  observation  led  to  the  realization  that  one  prediction  curve  could  be  used 
to  describe  the  performance  of  many  vehicles  with  similar  characteristics  if  the  drawbar  versus 
RCI  curves  for  those  vehicles  were  all  plotted  such  that  their  zero  drawbar  points  coincided. 
The  zero  drawbar  points  occur  where  the  soil  strength  is  equal  to  VCIj,  and  relationships  for 
predicting  VCIj  were  already  available.  Therefore,  the  drawbar  performance  potential  of 
many  vehicles  with  widely  varying  VCIj  values  but  similar  basic  characteristics  could  be 
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predicted  using  a  single  prediction  curve  by  subtracting  the  predicted  VCIi  value  for  each 
vehicle  from  the  actual  RCI  value  as  shown  below. 

RCIx  =  RCI  -  VCI^ 

As  explained  in  Chapter  I,  the  standard  NRMM  relationship  for  drawbar  performance 
of  wheeled  vehicles  on  CH  soils  separates  vehicles  into  two  prediction  categories:  vehicles 
with  CPF  ^  4  psi  and  vehicles  with  CPF  <  4  psi.  The  NRMM  relationship  was  based  on 
drawbar  performance  at  20  percent  slip  until  the  early  1990s,  but  it  is  now  based  on  drawbar 
performance  at  100  percent  slip.  This  research  concentrates  on  20  percent  slip  performance 
for  reasons  explained  in  Chapter  E.  Therefore,  the  relationship  shown  below,  which  is  based 
on  standard  NRMM  logic  explained  by  Priddy  [1995],  was  used  to  translate  the  standard 
NRMM  relationship  from  100  percent  slip  to  20  percent  slip  performance  magnitude. 

The  three  drawbar  coefficients  on  the  right-hand  side  of  the  above  equation  can  be  determined 
from  standard  NRMM  relationships,  and  the  equations  used  to  describe  these  relationships 
are  shown  below  [Priddy  1995]  (note  that  the  new  prediction  curves  proposed  in  the 
referenced  report  are  now  standard  for  NRMM). 
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The  standard  NRMM  relationship  translated  to  20  percent  slip  performance  magnitude,  which 
resulted  from  consideration  of  the  equations  just  mentioned,  is  shown  below. 
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RCIx  was  calculated  for  each  of  the  observations  in  the  field  drawbar  database.  The 
best  CP2  relationship,  which  was  demonstrated  in  Chapter  TV  as  the  best  overall  relationship 
for  veil,  was  used  to  predict  VCIi  in  the  RCIx  calculations.  The  data  were  plotted  with  the 
standard  NRMM  equations  (translated  to  20  percent  slip)  on  a  x-y  plot  as  shown  in  Figure 
25.  The  figure  shows  that  the  data  are  described  fairly  well  by  the  standard  NRMM 
relationship,  but  there  is  a  large  amount  of  scatter  and  the  data  are  not  distributed  very  well 
throughout  the  RCIx  range,  a  lot  of  the  data  being  clustered  at  about  270  RCIx- 

The  standard  NRMM  relationship  was  developed  using  only  part  of  the  data  in  the 
field  drawbar  database.  New  relationships  developed  using  all  of  the  data  would  have  an 
immediate  advantage  over  the  NRMM  relationship  since  their  empirical  constants  would  be 
optimized  to  the  data.  Also,  the  relationship  had  to  be  translated  from  100  percent  slip  to  20 
percent  slip  performance  magnitude  as  previously  described.  These  facts  meant  that  fair 
comparisons  between  the  standard  NRMM  relationship  and  any  of  the  other  relationships 
evaluated  during  this  research  could  only  be  made  if  the  NRMM  relationship  was  optimized 
to  the  complete  data  set.  Therefore  a  modified  version  of  the  NRMM  relationship  was 


developed  that  optimized  the  empirical  constants  in  the  prediction  equations  for  drawbar  as 
a  function  of  RCIx-  The  combined  linear  regression  and  numerical  optimization  process  was 
used  to  determine  the  empirical  constants  in  the  equations  (i.e.,  curve-fits),  and  the  resulting 
optimized  NRMM  relationship  is  shown  below. 
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Figure  25.  Standard  NEIMM  relationship  for  drawbar 
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The  field  drawbar  data  were  plotted  with  the  optimized  NRMM  equations  on  a  x-y 
plot  as  shown  in  Figure  26.  The  figure  demonstrates  that  the  new  optimized  curve-fits 
describe  the  field  drawbar  data  better  than  the  standard  curve-fits,  especially  for  the  data  with 
CPF  less  than  4  psi.  It  is  emphasized  that  the  optimized  curve-fit  for  CPF  less  than  4  psi  is 
not  recommended  for  use  beyond  this  study  because,  as  the  figure  shows,  there  are  only  three 
observations  in  the  high  drawbar  range  and  no  observations  beyond  an  RCIx  of  160. 


Figure  26.  Optimized  NRMM  relationship  for  drawbar 

The  other  four  existing  models  for  drawbar  are  all  based  on  existing  numerics.  The 
first  three  numerics,  which  were  developed  by  WES,  were  founded  on  the  dimensional 
analyris  research  conducted  by  Freitag  [1965]  that  was  described  in  Chapter  I.  Freitag  used 
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the  laboratory  drawbar  data  described  in  Chapter  n  and  his  assodated  laboratory  sinkage  data 

to  develop  a  numeric  that  related  well  to  the  motion  resistance,  sinkage,  torque  output  and 

drawbar  performances  measured  in  his  laboratory  research.  Although  Freitag’s  numeric  was 

developed  from  evaluation  of  laboratory  performance  measurements,  it  is  not  limited  to 

laboratory  use  and  was  evaluated  in  this  research  as  a  parameter  for  describing  field  drawbar 

data.  Freitag’s  numeric  is  calculated  as  shown  below. 
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The  field  drawbar  data  were  plotted  versus  Freitag’s  numeric  on  a  x-y  plot,  and  the 
resulting  trend  was  in  a  form  that  could  be  described  well  by  the  rectangular  hyperbola 
equation.  Therefore,  the  combined  linear  regression  and  numerical  optimization  process  was 
used  to  develop  a  rectangular  hyperbola  equation  for  describing  field  drawbar  performance 
as  a  function  of  The  resulting  equation  is  shown  below  and  is  plotted  with  the  field 
drawbar  data  in  Figure  27. 
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The  figure  demonstrates  that  the  trend  is  hyperbolic  and  that  Freitag’s  numeric  describes  the 
field  drawbar  data  reasonably  well  with  the  exception  of  the  data  having  CPF  less  than  4  psi 
in  the  low  drawbar  range.  The  figure  also  demonstrates  that  there  is  a  much  better  distribu¬ 
tion  of  data  throughout  the  full  range  for  there  was  for  RCIx. 
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Figure  27.  Relationship  between  Freitag’s  Numeric  and  the  field  drawbar  data 


The  other  two  WES  numerics  involved  modifications  to  Freitag’s  original  numeric. 
Tumage  [1972]  showed  that  a  modified  form  of  Freitag’s  numeric  described  the  laboratory 
performance  of  several  tires  with  greatly  varying  shapes  better  than  Freitag’s  numeric. 
Tumage’s  data  included  the  four  tires  evaluated  in  Freitag’s  original  research  plus  a  few  other 
tires  with  different  shapes,  and  this  larger  mix  of  tire  shapes  led  him  to  incorporate  an 
additional  consideration  for  the  ratio  of  tire  section  width-to-diameter.  Tumage’s  numeric 
is  calculated  as  shown  below. 
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Figure  28.  Relationship  between  Tumage’s  Numeric  and  the  field  drawbar  data 


The  field  drawbar  data  were  plotted  versus  Tumage’s  numeric  on  a  x-y  plot,  and  the 
resulting  trend  was  in  a  form  that  could  be  described  well  by  the  rectangular  hyperbola 
equation.  Therefore,  the  combined  linear  regression  and  numerical  optimization  process  was 
used  to  develop  a  rectangular  hyperbola  equation  for  describing  field  drawbar  performance 
as  a  function  of  .  The  resulting  equation  is  shown  below. 
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A  x-y  plot  with  the  field  drawbar  data  and  the  equation  for  drawbar  versus  is  shown  in 
Figure  28.  The  figure  demonstrates  that  the  trend  is  hyperbolic  and  that  Tumage’s  numeric 
describes  the  field  drawbar  data  reasonably  well.  Tumage’s  numeric  is  shown  to  also  have 
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problems  describing  the  data  with  CPF  less  than  4  psi  in  the  low  drawbar  range,  but  it  does 
appear  to  be  slightly  better  than  Freitag’s  numeric  in  this  range.  The  data  are  shown  to  be 
well  distributed  throughout  the  full  range  of 

The  last  WES  numeric  was  developed  primarily  by  Smith  in  collaboration  with 
Tumage  [1995],  Smith’s  numeric  was  similar  to  Tumage’s  numeric  in  content  but  much 
different  in  form.  Smith  cast  the  ratios  of  tire  deflection-to-section  height  and  section  width- 
to-diameter  in  forms  different  from  both  Freitag’s  and  Tumage’s  numerics.  According  to 
Tumage  [1995],  Smith’s  numeric  should  be  the  best  WES  numeric  for  describing  field 
performance.  Smith’s  numeric  is  calculated  as  shown  below. 

_  RCln  bd(  ^  b 

The  field  drawbar  data  were  plotted  versus  Smith’s  numeric  on  a  x-y  plot,  and  the 
resulting  trend  was  in  a  form  that  could  be  described  well  by  the  rectangular  hyperbola 
equation.  Therefore,  the  combined  linear  regression  and  numerical  optimization  process  was 
used  to  develop  a  rectangular  hyperbola  equation  for  describing  field  drawbar  performance 
as  a  function  of  The  resulting  equation  is  shown  below  and  is  plotted  with  the  field 

drawbar  data  in  Figure  29. 


The  figure  demonstrates  that  the  trend  is  hjqierbolic  and  that  Smith’s  numeric  describes  the 
field  drawbar  data  reasonably  well.  Smith’s  numeric  is  shown  to  also  have  slight  problems 
describing  the  data  with  CPF  less  than  4  psi  in  the  low  drawbar  range,  but  it  appears  to  be 
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better  than  both  Freitag’s  and  Tumage’s  numerics  in  this  range.  ‘s  also  shown  to 

produce  a  good  distribution  throughout  the  full  range  of  data. 


Figure  29.  Relationship  between  Smith’s  Numeric  and  the  field  drawbar  data 


The  last  existing  model  is  based  on  a  numeric  for  drawbar  that  was  recently  proposed 
for  use  in  the  United  Kingdom  by  Maclaurin  [1997].  Maclaurin’s  numeric  is  actually  loosely 
founded  on  Freitag’s  original  numeric,  but  its  development  was  accomplished  through  a 
significantly  different  approach.  Unlike  Tumage’s  and  Smith’s  numerics  which  involved 
modifications  to  Freitag’s  numeric,  Maclaurin’s  numeric  is  completely  original.  Maclaurin 
derived  his  numeric  fi-om  the  results  of  recent  drawbar  experiments  conducted  by  DERA  with 
a  field  deployable,  computer  controlled,  single- wheel  tire  tester.  The  formula  used  for  this 
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drawbar  model  is  actually  slightly  different  from  the  formula  found  in  the  referenced  publica¬ 
tion.  The  true  formula  had  the  constant  2  for  the  number  of  tires  per  axle,  but  the  variable 
n  was  used  here  for  reasons  discussed  in  Chapter  IV  pertaining  to  Maclaurin’s  MMP. 
Maclaurin’s  numeric  was  calculated  as  shown  below. 


The  field  drawbar  data  were  plotted  versus  Maclaurin’s  numeric  on  a  x-y  plot,  and  the 
resulting  trend  was  in  a  form  that  could  be  described  well  by  the  rectangular  hyperbola 
equation.  Therefore,  the  combined  linear  regression  and  numerical  optimization  process  was 
used  to  develop  a  rectangular  hyperbola  equation  for  describing  field  drawbar  performance 
as  a  function  of  Nj^.  The  resulting  equation  is  shown  below  and  is  plotted  with  the  field 
drawbar  data  in  Figure  30. 
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The  figure  demonstrates  that  the  trend  is  hyperbolic  and  that  Maclaurin’s  numeric  describes 
the  field  drawbar  data  reasonably  well.  Just  like  the  three  WES  numerics,  Maclaurin’s 
numeric  is  shown  to  have  problems  describing  the  data  with  CPF  less  than  4  psi  in  the  low 
drawbar  range.  just  like  the  other  existing  numerics,  is  shown  to  produce  a  good  distri¬ 
bution  throughout  the  full  range  of  data. 
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Figure  30.  Relationship  between  Maclaurin’s  Numeric  and  the  field  drawbar  data 

New  Development  Models 

As  stated,  four  major  types  of  new  development  models  were  explored  in  the  drawbar 
modeling.  These  included  a  bearing  ratio  model,  an  excess  bearing  ratio  model,  an  excess 
RCI  ratio  model,  and  a  numeric-type  model.  Only  the  numeric-type  model  will  be  discussed, 
and  only  the  major  variants  of  the  numeric-type  model  will  be  presented.  Before  the  individ¬ 
ual  model  variants  are  addressed,  clarification  is  provided  on  what  is  meant  by  ‘bearing  ratio 
model’,  ‘excess  bearing  ratio  model’,  and  ‘excess  RCI  ratio  model’  (see  first  few  sections  on 
new  development  models  in  Chapter  IV  for  clarification  on  numeric-type  model ). 
Explanations  are  also  provided  as  to  why  these  types  of  models  were  explored  in  the  drawbar 
modeling  effort  and  why  the  numeric-type  model  is  the  only  one  presented. 
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The  beaiing  ratio  model  is  based  on  a  dimensionless  parameter  that  is  the  ratio  of  soil 
strength  over  contact  pressure.  Therefore,  it  represents  the  ratio  of  resistance  over  loading 
for  the  bearing  interaction.  The  bearing  ratio  parameter  closely  resembles  the  numeric 
parameter  for  drawbar,  but  the  bearing  ratio  parameter  uses  a  logically  derived  expression  for 
contact  pressure  similar  to  the  logical  derivations  of  contact  pressure  used  for  the  VCIj 
modeling  in  Chapter  IV.  In  fact,  the  variants  explored  for  the  bearing  ratio  model  used  the 
same  expressions  for  contact  pressure  that  resulted  from  the  VCIj  modeling  (i.e.,  CP2,  etc.). 
The  excess  bearing  ratio  model  is  essentially  the  same  with  the  exception  that  it  uses  excess 
soil  strength  (RCIx)  rather  than  the  actual  soil  strength  (RCI).  The  excess  RCI  ratio  model 
is  somewhat  similar,  but  it  uses  the  ratio  of  excess  soil  strength  over  VCIj.  All  three  of  these 
types  of  models  are  based  in  part  on  the  ideas  utilized  in  the  standard  NRMM  drawbar 
relationship  and  Freitag’s  original  dimensional  analysis  research. 

The  bearing  ratio  model  was  explored  for  drawbar  because  of  the  decent  level  of 
success  demonstrated  in  the  existing  numerics  and  the  fact  that  Freitag’s  numeric  worked 
reasonably  well  describing  VCIj  which  was  proven  to  be  almost  entirely  controlled  by  contact 
pressure.  However,  the  bearing  ratio  model  did  not  work  very  well  (as  can  be  seen  in 
Appendbc  C),  and  this  led  to  the  idea  of  trying  excess  bearing  ratio.  It  was  hypothesized  that 
the  collapang  effect  of  RCIx  on  drawbar  versus  RCI  might  produce  a  beneficial  effect  on  the 
bearing  ratio  model  in  a  similar  feshion.  The  modeling  demonstrated  that  this  hypothesis  was 
false,  and  in  fact,  the  excess  bearing  ratio  model  was  worse  than  the  bearing  ratio  model. 

The  excess  RCI  ratio  model  was  explored  due  to  the  success  demonstrated  in  the 
RCIx  models  (standard  and  optimized).  Dividing  RCIx  by  VCIi  has  two  logical  benefits. 
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First,  it  makes  the  modeling  parameter  dimensionless  forcing  the  units  in  the  relationship  to 
be  correct  (recall  that  the  drawbar  coeflBdent  performance  metric  is  dimensionless).  Second, 
it  allows  the  principal  vehicle  characteristics  to  have  influence  on  drawbar  performance 
throughout  the  full  range  of  soil  strengths  rather  than  just  at  zero  drawbar.  The  excess  RCI 
ratio  model  demonstrated  better  potential  than  the  bearing  ratio  and  excess  bearing  ratio 
models,  but  it  did  not  work  very  well  either,  falling  short  of  the  quality  demonstrated  in  the 
optimized  RCIx  model. 

The  numeric-type  model  was  explored  due  to  the  level  of  success  demonstrated  in 
the  existing  numerics  and  due  to  the  success  of  the  numeric-type  model  in  the  VCIi  modeling. 
As  will  be  shown,  the  numeric-type  model  was  veiy  successful  in  the  drawbar  modeling.  The 
numeric-type  model  not  only  worked  well  describing  the  field  drawbar  data,  it  also  provided 
some  potential  insight  into  the  controlling  mechamsm  in  drawbar  performance.  Therefore, 
the  m^or  variants  of  the  numeric-type  model  were  the  focus  of  the  discussion  for  the  drawbar 
modeling  effort. 

As  described  in  Chapter  IV,  a  numeric-type  model  allows  the  data  to  completely 
define  the  form  of  a  simple,  principal  parameter  (referred  to  as  a  numeric)  under  the  constraint 
that  the  parameter’s  units  must  be  correct.  The  drawbar  coefiBcient  performance  metric  is 
dimensionless,  and  hence  the  numeric  parameter  must  also  be  dimensionless.  The  soil  part 
of  the  numeric  was  soil  strength  measured  in  RCI  which  introduces  pressure  units  (i.e., 
force/length^.  For  the  vehicle  part  of  the  numeric,  w/n  was  used  for  including  force  units  and 
the  prindpal  tire  shape  dimensions  were  used  to  include  the  length^  units.  This  basic  approach 
was  used  in  all  of  the  individual  numeric  variants.  The  only  differences  in  the  variants  were 
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the  number  of  prindpal  tire  shape  dimensions  included  and  the  approach  used  to  develop  their 
optimum  arrangement. 

For  each  of  the  numeric-type  model  variants,  the  same  process  was  used  to  develop 
the  relationship  between  the  drawbar  performance  metric  and  the  numeric.  Once  the  basic 
mathematical  form  of  the  drawbar  numeric  was  decided,  a  statistical  model  was  used  to 
develop  the  relationship.  The  rectangular  hyperbola  equation  was  used  in  the  relationships 
because  of  its  strength  in  describing  the  relationship  for  all  of  the  existing  models.  Due  to  the 
performance  differences  that  were  observed  between  bias  ply  and  radial  tire  configurations 
in  the  VCIj  modeling,  tire  construction  was  considered  in  the  drawbar  modeling,  but  there 
was  no  ^parent  difference  in  performance  between  bias  ply  and  radial  tire  configurations  in 
the  field  drawbar  database. 

The  first  numeric-type  model  variant  used  an  approach  for  the  vehicle  part  that  was 
identical  to  the  approach  used  for  the  best  numeric  variant  for  VCI,.  The  approach  was  to 
use  a  product  of  b,  d,  and  three  dimensionless  tire  ratios  where  each  ratio  was  raised  to  a 
power  (Pi ,  Pj’  ratios  chosen  used  d  as  the  characteristic  length  in  dimension¬ 

less  ratios  with  each  of  b,  h,  and  6  so  that  all  four  of  the  principal  tire  shape  dimensions  could 
be  included  in  the  final  numeric.  Three  ratios  (i.e.,  powers)  were  used  because  there  were 
four  length  variables  less  one  constraint.  This  modeling  approach  allowed  inclusion  of  b,  d, 
h,  and  6  in  the  final  numeric  with  each  raised  to  its  own  numerically  optimized  power  under 
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the  constraint  that  the  summation  of  the  four  powers  must  equal  two.  The  form  of  the  first 
numeric-type  model  variant  is  shown  below. 


=  B  +  Pi 

=  ^(nbd(,b/d)'’'(d/h)’'‘(6/d)'^^ 


The  combined  linear  regression  and  numerical  optimization  process  was  used  to 
determine  the  empirical  constants  in  the  relationship.  The  process  resulted  in  being  -0.90, 
/?2  bang  0.82,  and  being  0.45.  The  resulting  final  equations  are  shown  below  where  the 
numeric  has  been  simplified  by  collapsing  the  ratios  and  collecting  like  terms. 
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Nj^j  was  calculated  for  each  observation  in  the  field  drawbar  database,  and  the  data  were 
plotted  with  the  equation  describing  the  relationship  on  a  x-y  plot  as  shown  in  Figure  3 1 .  The 
figure  shows  that  the  field  drawbar  data  are  described  very  well  by  the  first  numeric-type 
model  variant,  and  unlike  all  the  existing  numerics,  N^j  had  no  problem  describing  the  data 
with  CPF  <  4  psi  in  the  lower  drawbar  range. 
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Figure  31.  Relationship  between  new  numeric  variant  1  and  the  field  drawbar  data 

Due  to  the  small  exponent  on  section  width  (i.e.,  0.10)  in  N^j,  it  was  hypothesized 
that  section  height  may  have  caused  a  reduced  influence  for  section  width  as  a  result  of  the 
interrelation  between  these  two  variables  for  the  majority  of  field  drawbar  observations  (see 
Chapter  E).  For  this  reason,  the  second  numeric-type  model  variant  was  developed  without 
section  height  included.  The  fact  that  Maclaurin’s  numeric  did  not  consider  section  height 
also  supported  the  removal  of  this  variable.  The  approach  was  to  use  a  product  of  b,  d,  and 
two  dimensionless  tire  ratios  where  both  ratios  were  raised  to  a  power  (pj  and  p^).  The 
ratios  chosen  used  d  as  the  characteristic  length  in  dimensionless  ratios  with  both  b  and  5. 
Two  ratios  (i.e.,  powers)  were  used  because  there  were  three  length  variables  less  one 
constraint.  This  modeling  approach  allowed  inclusion  of  b,  d,  and  6  in  the  final  numeric  with 
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each  raised  to  its  own  numerically  optimized  power  under  the  constraint  that  the  summation 
of  the  three  powers  must  equal  two.  The  form  of  the  second  numeric-type  model  variant  is 
shown  below. 

=  6  +  Pi 

^20  ^D4  *  P2 

The  combined  linear  regression  and  numerical  optimization  process  was  used  to 
determine  the  empirical  constants  in  the  relationship,  and  the  process  resulted  in  /7j  being  -1.0 
and  /?2  being  0.50.  The  resulting  final  equations  are  shown  below  where  the  numeric  has 
been  simplified  by  collapsing  the  ratios  and  collecting  like  terms. 

ST'  .  0.69  -  - 

was  calculated  for  each  observation  in  the  field  drawbar  database,  and  the  data  were 
plotted  vwth  the  equation  describing  the  relationship  on  a  x-y  plot  as  shown  in  Figure  32.  The 
figure  shows  that  the  field  drawbar  data  are  described  very  well  by  this  model  variant,  and 
is  shown  to  describe  the  data  with  a  form  and  quality  similar  to  that  of  Nj^j . 

The  resulting  form  of  was  surprising.  Section  height  was  removed  from  consid¬ 
eration  in  the  numeric  under  the  hypothesis  that  section  vddth  would  then  demonstrate  a 
stronger  influence  (i.e.,  a  larger  exponent),  but  instead  the  influence  of  section  width 
completely  disappeared.  This  result  indicates  that  section  width  plays  a  very  small  role  in 
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drawbar  performance,  and  hence  contact  pressure  is  not  the  controlling  factor  in  drawbar 
performance.  Tlus  conclusion  is  supported  by  Tumage’s  [1972]  findings  that  a  bearing  ratio 
parameter  using  actual  measured  hard-surface  contact  area  did  not  adequately  describe 
drawbar  performance  in  clay  soils,  and  it  is  also  supported  by  the  poor  results  obtained  in  this 
research  for  the  bearing  ratio  models.  Obviously,  this  conclusion  is  not  valid  at  the  zero 
drawbar  point  (i.e.,  VCIj),  but  these  results  indicate  that  a  transition  in  the  controlling 
mechanism  occurs  very  shortly  beyond  zero  drawbar  and  that  something  other  than  contact 
pressure  controls  throughout  the  majority  of  the  drawbar  range. 


Figure  32.  Relationship  between  new  numeric  variant  4  and  the  field  drawbar  data 
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The  third  numeric-type  model  variant  was  developed  based  on  a  closer  look  at  the 
form  of  Nj^.  It  was  observed  that  could  be  rewritten  as  shown  below. 

-  ^[ndm] 

It  was  also  recognized  that  is  very  similar  to  the  estimate  for  hard-surface  contact  length 
that  results  from  ^plication  of  the  long  chord  analogy  (i.e.,  Lg  =  2yjdh  -  6^  ).  Therefore, 
it  was  hypothesized  that  yfdb  could  be  replaced  by  the  logically  derived  Lg  in  the  numeric 
parameter  \rithout  affecting  the  form  and  quality  of  the  relationship.  The  form  of  the  third 
numeric-type  model  variant  is  shown  below. 

Z)^c/  _  ^  p, 

¥  +B 

"  20  "2 

The  combined  linear  regression  and  numerical  optimization  process  was  used  to 
determine  the  empirical  constants  in  the  relationship,  and  the  resulting  final  equations  are 
shown  below. 


Nn8  = 


RCI 


"  ^D8 


^wt/Lgj  where  Zg  =  l^jdb  -  6^ 


DRCi  5  0 

—  =  0.692  - 


Nr.o  was  calculated  for  each  observation  in  the  field  drawbar  database,  and  the  data  were 

Do 


plotted  with  the  equation  describing  the  relationship  on  a  x-y  plot  as  shown  in  Figure  33.  The 
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figure  verifies  that  describes  the  data  well  and  with  a  form  and  quality  nearly  identical 
to  that  of 


Figure  33.  Relationship  between  new  numeric  variant  8  and  the  field  drawbar  data 

An  interesting  conclusion  resulted  from  the  third  numeric-type  model  variant.  It 
appears  that  the  product  of  tire  diameter  and  hard-surface  contact  length  represents  the 
vehicle  side  of  the  controlling  mechanism  in  drawbar  performance.  This  conclusion  makes 
more  sense  if  it  is  realized  that  tire  diameter  could  be  replaced  with  tire  radius  without  any 
effect  on  the  form  and  quality  of  the  relationship.  The  product  of  tire  radius  and  hard-surface 
contact  length  indicates  that  the  controlling  mechanism  in  drawbar  performance  may  be  the 
length  of  the  torque  moment  arm  (i.e.,  active  radius)  times  the  length  of  soil  in  shear. 
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In  order  to  gain  greater  confidence  in  the  validity  of  as  a  predictor  for  drawbar 
performance,  additional  drawbar  data  were  explored.  Freitag’s  numeric,  which  includes  a 
strong  influence  fi'om  section  width  and  was  developed  entirely  fi’om  laboratory  data,  was 
shown  to  describe  the  field  drawbar  data  fairly  well.  Therefore,  Nj^g,  which  completely 
ignores  section  width  and  was  developed  entirely  fi'om  field  data,  was  evaluated  using  labora¬ 
tory  drawbar  data.  The  laboratory  drawbar  data  (described  in  Chapter  H)  fi'om  Freitag’s 
original  dimensional  analysis  research  were  used  for  the  evaluation. 

Before  conducting  the  evaluation  of  N^g  with  the  laboratory  data,  the  difference 
between  field  and  laboratory  drawbar  measurements  was  observed.  Figure  34  shows  a  x-y 
plot  of  drawbar  coeflScient  versus  Freitag’s  numeric.  The  plot  contains  the  laboratory 
drawbar  measurements  and  the  curve-fit  developed  earlier  for  describing  the  relationship 
between  the  field  drawbar  measurements  and  Freitag’s  numeric.  The  figure  shows  that  the 
laboratory  drawbar  measurements  fallow  a  trend  that  is  very  different  fi'om  the  trend  in  the 
field  drawbar  measurements  in  terms  of  form  and  magnitude.  Based  on  a  review  of  Figure 
27,  it  is  probable  that  the  discrepancy  below  a  drawbar  coefBcient  of  0.40  is  probably  not  as 
large  as  indicated  in  Figure  34,  but  the  enormous  discrepancy  above  a  drawbar  coefBcient  of 
0.40  is  undeniable.  Despite  the  difference  between  drawbar  performances  observed  in  the  two 
environments,  the  quality  with  which  Freitag’s  numeric  described  the  field  data  indicates  that 
the  same  modeling  parameter  may  describe  the  data  well  in  both  environments,  although  with 
very  different  curve-fits. 
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Figure  34.  Comparison  between  laboratory  and  field  drawbar  measurements 

For  evaluating  Nj^g  in  the  laboratory  environment,  comparison  x-y  plots  were  devel¬ 
oped  as  shown  in  Figure  35.  The  upper  plot  (Figure  35a)  shows  the  laboratory  drawbar  data 
plotted  versus  Freitag’s  numeric,  and  the  lower  plot  (Figure  35b)  shows  the  same  data  plotted 
versus  Nj^g.  Although  the  scales  are  much  different  between  the  upper  and  lower  plots,  the 
absolute  horizontal  and  vertical  spread  of  the  data  are  essentially  the  same.  Therefore,  an 
adequate  visual  assessment  for  the  difference  in  quality  between  the  two  numerics  can  be 
made  using  the  figure. 

Figure  35  demonstrates  that  Freitag’s  numeric  describes  the  laboratory  drawbar  data 
better  than  Nj^g.  The  quality  of  N^g  is  not  too  bad  though,  and  it  collapses  the  data  for  three 
of  the  tires  into  a  central  trend  very  similar  in  quality  and  form  to  that  of  Freitag’s  numeric. 


Drawbar  Coefficient  (Lab)  P  Drawbar  Coefficient  (Lab) 


Figure  35.  New  numeric  variant  8  versus  Freitag’s  Numeric  in  laboratory  drawbar  data 
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However,  the  drawbar  performance  of  the  4.00-20  tire  follows  a  significantly  different  trend 
when  Nj^  is  used.  The  4.00-20  tire  had  a  veiy  small  carcass  area  relative  to  its  diameter,  and, 
as  illustrated  in  Figure  6,  it  had  a  very  small  section  height  relative  to  its  diameter  when 
compared  to  the  other  tires  fi'om  all  three  of  the  databases.  This  observation  led  to  the 
hypothesis  that  the  ratio  of  tire  section  height  to  diameter  could  be  used  as  a  potentially 
improving  modification  to  the  numeric  parameter.  It  was  noted  that  this  modification  should 
have  no  detrimental  effect  on  the  quality  of  the  numeric  in  describing  the  field  data  since  all 
of  the  tires  fi'om  the  field  drawbar  database  had  h/d  ratios  that  were  very  similar  in  magnitude 
(note  that  this  would  not  have  been  the  case  for  b/d). 

The  last  numeric-type  model  variant  was  developed  for  evaluating  the  potential  benefit 
of  the  h/d  modification.  This  model  variant  used  a  principal  parameter  that  was  the  product 
of  N^g  and  h/d  raised  to  a  power  (pj ).  The  form  of  this  model  variant  is  shown  below. 


Atrial  and  error  approach  was  used  to  determine  the  value  of  that  best  collapsed  all  four 
of  the  laboratory  tires  into  a  central  trend,  and  a  value  of  0.6  was  visually  assessed  to  work 
best.  Figure  36  shows  comparison  x-y  plots  in  a  layout  identical  to  that  of  Figure  35.  Once 
again,  the  absolute  horizontal  and  vertical  spread  of  the  data  in  the  upper  plot  (Figure  36a) 
and  the  lower  plot  (Figure  36b)  are  essentially  the  same.  Therefore,  Figure  36  allows  a  visual 
assessment  for  the  difference  in  quality  between  Freitag’s  numeric  and  to  be  made.  The 
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Freitag's  Numeric 


Performance  of  Freitag’s  Numeric 


Figure  36.  New  numeric  variant  16  versus  Freitag’s  Numeric  in  laboratory  drawbar  data 
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figure  demonstrates  that  there  is  no  significant  difference  in  quality  and  form  between  the  two 
numerics. 


Once  the  appropriate  exponent  for  the  h/d  modification  had  been  determined  using  the 
laboratory  data,  the  combined  linear  regression  and  numerical  optimization  process  was  used 
to  determine  the  other  empirical  constants  in  the  relationship  using  the  field  data.  The 
resulting  final  equations  are  shown  below. 
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N^jg  was  calculated  for  each  observation  in  the  field  drawbar  database,  and  the  data  were 
plotted  with  the  equation  describing  the  relationship  between  drawbar  and  Nj^jg  on  a  x-y  plot 
as  shown  in  Figure  37.  The  figure  verifies  that  Nj^jg  describes  the  data  well  and  with  a  form 
and  quality  nearly  identical  to  that  of  N^jg . 
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Figure  37.  Relationship  between  new  numeric  variant  16  and  the  field  drawbar  data 


Model  Comparisons 

Several  existing  and  new  drawbar  relationships  were  presented  in  the  preceding 
sections  of  this  chapter,  and  a  simple  x-y  plot  of  each  relationship  was  provided.  The  x-y 
plots  allowed  a  visual  interpretation  of  the  quality  with  which  each  relationship  described  the 
measured  data,  but  higher  accuracy  methods  of  assessing  quality  were  necessary  to  make 
sound  comparisons  between  the  various  relationships.  For  this  purpose,  the  standard  error 
(SJ  and  adjusted  coeffident  of  determination  (Adj.  measures  of  quality  were  used.  These 
quality  measures  were  determined  during  the  development  of  each  of  the  many  model 
variants,  and  they  were  used  as  optimization  targets  in  models  that  required  numerical 
optimization  as  described  in  Chapter  HI.  Table  6  lists  these  measures  along  with  their 
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assodated  degrees  of  freedom  for  each  of  the  relationships.  The  quality  measures  shown  in 
the  table  are  based  on  a  consideration  of  the  field  drawbar  data  only  (i.e.,  no  laboratory  data). 

Table  6.  Quality  comparison  for  the  drawbar  modeling  results  (field  data  only) 


Relationship 

Degrees  of 
Freedom 

Standard 

Error 

1 

Adjusted 

R-squared 

Standard  NRMM  (RCI^) 

81 

0.097 

77.8 

Optimized  NRMM 

81 

0.079 

85.3 

Freitag’s  Numeric 

80 

0.095 

78.7 

Tumage’s  Numeric 

80 

0.085 

82.8 

Smith’s  Numeric 

79 

0.082 

84.1 

Maclaurin’s  Numeric 

78 

0.090 

80.6 

New  Numeric  Variant  1 

78 

0.063 

90.5 

New  Niuneric  Variant  4 

79 

0.063 

90.5 

New  Numeric  Variant  8 

80 

0.063 

90.6 

New  Numeric  Variant  16 

79 

0.064 

90.2 

Several  comparisons  can  be  made  using  Table  6.  First,  the  table  demonstrates  that  the 
optimized  NRMM  relationship  described  the  data  significantly  better  (i.e.,  lower  and 
higher  than  the  standard  NRMM  relationship.  For  the  three  existing  WES  numerics,  the 
table  shows  that  Tumage  made  significant  improvement  to  Freitag’s  numeric  and  that  Smith, 
in  collaboration  with  Tumage,  made  slightly  more  improvement.  Maclaurin’s  numeric,  which 
was  the  existing  numeric  developed  recently  by  DERA,  described  the  data  better  than 
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Freitag’s  numeric  but  not  as  well  as  Tumage’s  numeric.  In  terms  of  the  new  relationships, 
all  of  the  new  drawbar  numerics  essentially  described  the  data  with  equal  quality.  Drawbar 
numeric  variant  8,  which  included  the  logical  ^ ,  ranked  best,  and  drawbar  numeric  variant 
16,  which  included  consideration  of  the  laboratory  drawbar  data,  ranked  a  very  close  last. 
Most  significantly,  the  table  demonstrates  that  the  new  drawbar  numerics  describe  the  data 
much  better  than  any  of  the  existing  relationships.  The  optimized  NRMM  relationship 
described  the  data  better  than  any  of  the  existing  numerics,  but  its  quality  was  significantly 
worse  than  the  new  numerics.  These  results  indicate  that  it  is  definitely  feasible  to  develop 
new  drawbar  relationships  that  meet  the  decision  criteria  stated  in  the  objectives  for  this 
research.  The  results  also  indicate  that  the  driving  factor  in  drawbar  performance  may  be  the 
product  of  the  torque  moment  arm  at  the  tires  (i.e.,  active  radius)  and  the  length  of  soil  in 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

Summary 

The  current  traction  relationships  in  NRMM  describe  the  performance  of  real  vehicles 
in  naturally  occurring  soils  reasonably  well,  but  they  have  some  major  weaknesses.  The 
weaknesses  stem  primarily  from  an  inadequate  and  sometimes  illogical  consideration  of  the 
effects  that  vehicle  characteristics  have  on  traction  performance.  This  research  sought  to 
demonstrate  the  feasibility  of  developing  better  traction  relationships  for  NRMM  through 
exploitation  of  existing  performance  data.  The  approach  used  to  demonstrate  the  feasibility 
was  to  conduct  a  statistical  modeling  effort  for  developing  some  new  traction  relationships. 
The  applied  modeling  effort  encompassed  a  systematic  and  extensive  exploration  of  new  and 
old  traction  modeling  concepts  that  rely  on  relatively  simple  characterizations  of  the  vehicle 
and  the  soil.  The  effort  concentrated  on  VCIj  and  drawbar  performance  characteristics  of 
pneumatic-tired  vehicles  interacting  with  non-slippery,  highly  plastic  clay  soils  (CH  classi¬ 
fication  by  uses). 

Readily  accessible  sources  of  VCIi,  field  drawbar,  and  laboratory  drawbar  perform¬ 
ance  data  were  acquired.  The  range  and  distribution  for  each  of  the  critical  system  variables 
assodated  with  the  performance  data  were  observed,  and  possible  interrelations  between  the 
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system  variables  were  explored.  The  acquired  data  provided  a  substantial  number  of  observa¬ 
tions  from  sources  that  spanned  a  period  of  about  three  decades,  and  they  incorporated  a 
broad  range  of  values  with  decent  distributions  for  the  critical  system  variables.  The  data  also 
demonstrated  that  there  are  some  trends  in  tire  shape  that  predominate  in  off-road  tire  design. 

Throughout  the  modeling  effort,  some  basic  concepts  were  considered.  One  of  these 
concepts  was  a  simplified  force  analysis  of  tire/soil  interaction  that  provides  the  fundamental 
logic  behind  the  NRMM  traction  algorithms.  Another  concept  was  a  mathematical  analogy 
for  the  effect  of  tire  deflection  on  contact  area  at  the  tire/soil  interface.  Some  basic  analytical 
tools  were  used  throughout  the  research  as  well.  These  tools  were  necessary  to  develop  the 
new  relationships,  and  they  included  least  squares  linear  regression,  iterative  numerical 
optimization  routines,  and  a  combined  process  of  both. 

In  the  veil  modeling,  several  existing  and  new  models  were  explored.  The  existing 
models  included  the  standard  NRMM  relationship  (MI  with  DCF),  three  versions  of  MMP, 
and  Frdtag’s  numeric  cast  in  a  VCI  form.  The  new  models  included  various  contact  pressure 
models  and  a  numeric-type  model.  The  contact  pressure  models  were  based  on  the  theory 
that  tire  contact  pressure  dominates  VCIi  performance.  The  numeric-type  model  allowed  the 
data  to  completely  define  the  form  of  a  simple,  single  system  parameter  referred  to  as  a 
numeric.  Several  variants  were  explored  for  most  of  the  new  types  of  VCIi  models. 

The  results  of  the  VCIi  modeling  demonstrated  that  several  new  relationships 
described  the  data  better  than  any  of  the  existing  relationships.  The  two  best  variants  of  the 
contact  pressure  models  provided  the  best  relationships,  and  the  best  numeric  variant  ranked 
a  close  third.  The  standard  NRMM  relationship  was  better  than  all  three  versions  of  MMP 
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and  Frdtag’s  numeric  cast  in  a  VCI  form,  but  even  a  modified  version  of  the  NRMM  relation¬ 
ship  that  was  fully  optimized  to  the  VCIj  data  could  not  match  the  quality  of  the  best  new 
relationships. 

Several  existing  and  new  models  were  also  explored  in  the  drawbar  modeling.  The 
existing  models  included  the  standard  NRMM  relationship  (RCIx),  three  WES  numerics 
(Freitag,  Tumage,  and  Smith),  and  a  recent  DERA  numeric  (Maclaurin).  The  new  models 
included  a  bearing  ratio  model,  an  excess  bearing  ratio  model,  an  excess  RCI  ratio  model,  and 
a  numeric-type  model.  The  bearing  ratio  model  was  based  on  the  ratio  of  soil  strength  in  RCI 
over  contact  pressure,  and  the  excess  bearing  ratio  model  used  the  ratio  of  RCIx  over  contact 
pressure.  The  RCI  ratio  model  used  the  ratio  of  soil  strength  in  RCI  over  VCIj,  and  the 
numeric-type  model  allowed  the  data  to  completely  define  the  form  of  a  simple,  single  system 
parameter  referred  to  as  a  numeric.  Several  variants  were  explored  for  all  of  the  new  types 
of  drawbar  models. 

Of  the  existing  models,  the  standard  NRMM  relationship  described  the  field  drawbar 
data  worst  and  Smith’s  numeric  did  best.  However,  a  modified  version  of  the  NRMM  rela¬ 
tionship  that  was  fully  optimized  to  the  field  drawbar  data  demonstrated  that  RCIx  could 
describe  the  data  better  than  any  of  the  four  existing  numerics.  Of  the  new  models,  only  the 
numeric-type  model  proved  useful,  and  nearly  all  of  the  new  numeric  variants  described  the 
field  drawbar  data  better  than  the  optimized  NRMM  relationship  and  all  of  the  existing 
numerics.  In  addition  to  describing  the  field  drawbar  data  better  than  RCIx  and  the  existing 
numerics,  one  of  the  new  numerics  even  described  the  laboratory  drawbar  data  vdth  a  quality 
essentially  equal  to  that  of  Freitag’ s  numeric. 
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Conclusions 

Based  on  the  results  obtained  from  the  modeling  efibrt,  several  significant  conclusions 

were  made  as  follows. 

1 .  The  fea^ility  exists  to  develop  better  traction  relationships  by  utili2dng  existing  traction 
performance  data.  Moreover,  better  relationships  can  be  developed  that  are  logical  in 
arrangement,  have  theoretical  foundations,  are  dimensionally  correct,  provide  a  more 
sound  understanding  of  the  role  that  vehicle  characteristics  have  on  traction  performance 
throughout  the  full  range  of  soil  strengths,  and  describe  the  existing  data  better  than  the 
current  NRMM  relationships. 

2.  veil  performance  is  completely  dominated  by  soil  bearing  capacity  with  tire  contact 
pressure  acting  as  the  driving  force.  Any  parameter  that  is  highly  correlated  with  the  true 
average  contact  pressure  under  the  tires  on  a  wheeled  vehicle  can  be  used  as  a  sound 
predictor  of  VCIj  performance. 

3.  There  is  a  difference  in  VCIj  performance  resulting  strictly  from  a  difference  in  tire 
construction  such  that  radial  tires  outperform  bias  ply  tires  by  about  2  psi  on  average. 

4.  The  best  CP2  relationship  is  the  best  relationship  for  predicting  VCIi  performance  out  of 
all  the  relationships  that  were  evaluated. 

5.  There  is  a  change  in  the  controlling  mechanism  such  that  drawbar  performance  beyond 
a  soil  strength  of  VCIi  is  not  controlled  by  tire  contact  pressure. 

6.  It  appears  that  the  dominant  factor  in  drawbar  performance  may  be  the  product  of  the 
torque  moment  arm  at  the  tires  (i.e.,  active  radius)  and  the  length  of  soil  in  shear. 
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7.  Numeric-type  models  can  be  extremely  useful  for  developing  relationships  in  areas  where 
little  is  known  about  the  controlling  mechanism. 

8.  A  shortcoming  in  the  Numeric  Approach  laboratory  research  resulted  from  the  fact  that 
an  attempt  was  made  to  develop  a  single  numeric  that  would  represent  the  controlling 
vehicle  influence  for  not  only  drawbar  but  also  for  sinkage,  motion  resistance,  and  torque 
output.  This  research  demonstrated  that  the  controlling  vehicle  influence  is  not  necessar¬ 
ily  the  same  for  all  four  of  the  performance  considerations. 


CHAPTER  Vn 


DIRECTIONS  FOR  FURTHER  RESEARCH 

The  results  of  this  thesis  have  demonstrated  great  promise  for  substantial  improvements 
to  the  traction  prediction  c^abilities  in  NRMM.  However,  much  more  research  is  necessary 
before  modifications  to  NRMM  should  be  considered.  Further  research  should  be  conducted 
for  veil  performance  of  wheeled  vehicles  in  non-slippery  CH  soils.  Similar  research  should 
also  be  conducted  for  VCIi  performance  of  wheeled  vehicles  in  the  other  major  USCS  soil 
types  and  for  tracked  vehicles  in  all  the  major  USCS  soil  types.  There  are  existing  data  fi'om 
previous  WES  studies  that  could  be  used  for  these  efforts.  Further  research  should  be 
conducted  for  traction  performance  above  VCIi  as  well,  but  the  efforts  should  take  a  much 
different  direction  than  the  current  NRMM  implementation. 

Further  research  should  be  conducted  to  gain  a  better  understanding  of  the  manner  in 
which  the  prindpal  tire  variables  affect  hard-surface  contact  length,  loaded  section  width,  and 
hard-surface  contact  area.  It  is  apparent  that  the  better  a  parameter  describes  hard-surface 
contact  area,  the  better  it  will  describe  VCIj  performance.  Efforts  should  also  be  conducted 
to  leverage  past  research  into  the  nature  of  stress  distributions  under  tires  in  soil  and  into  the 
deflected  shape  of  tires  in  soil.  WES  has  published  research  in  this  area,  and  information  fi’om 
these  and  other  similar  studies  could  probably  provide  useful  insight. 
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Potential  also  exists  for  improving  the  soil  ade  of  the  VCIj  problem.  A  review  of  the  root 
multi-pass  data  for  each  of  the  observations  used  during  this  research  should  be  conducted. 
A  digital  database  that  contdns  the  root  multi-pass  data  for  most,  if  not  all,  of  the  vehicle 
configurations  is  currently  available  at  WES,  and  this  database  would  greatly  simplify  the 
review.  Better  and  simpler  critical  layer  concepts  would  be  the  primary  goal.  A  secondary 
goal  would  be  to  verify  the  VCIj  measurements  for  all  of  the  observations  prior  to  establish¬ 
ing  new  relationships  for  NRMM  (some  of  the  observations  appear  questionable). 

The  current  NRMM  implementation  for  traction  beyond  VCIj  uses  drawbar  relationships 
and  motion  resistance  relationships  to  predict  tractive  force  potential  for  a  particular  vehicle/ 
terrain  scenario.  This  implementation  probably  originated  from  the  fact  that  drawbar  and 
motion  resistance  are  measurable  quantities  that  were  being  researched  before  a  speed- 
predicting  mobility  model  was  ever  developed.  It  probably  also  had  origins  fi’om  a  considera¬ 
tion  that  motion  resistance  on  powered  traction  elements  could  be  significantly  different  fi'om 
that  on  non-powered  traction  elements.  A  better,  more  logical  implementation  for  NRMM 
would  be  to  use  tractive  force  relationships  instead  of  drawbar  relationships. 

Tractive  force  is  a  terrain  reaction  force  that  occurs  at  the  vehicle/terrain  interface 
whereas  drawbar  is  a  resultant  effect  on  the  motion  (or  towing  ability)  of  the  vehicle.  The 
principal  vehicle  and  terrain  variables  have  a  direct  influence  on  tractive  force  and  motion 
resistance  and  only  a  secondary  influence  on  drawbar.  Tractive  force  relationships  could  be 
developed  by  adding  drawbar  and  motion  resistance  data  prior  to  developing  relationships 
rather  than  adding  the  two  relationships  as  is  done  in  NRMM  now.  The  error  associated  with 
using  motion  resistance  data  from  non-powered  traction  elements  to  estimate  motion 
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resistance  for  powered  traction  elements  would  probably  be  small  relative  to  the  other  sources 
of  error  (or  variability)  in  traction  measurements.  This  approach  would  probably  produce 
much  more  logical  relationships  and  a  more  sound  consideration  of  the  influence  that  vehicle 
characteristics  have  on  traction  performance. 

An  example  of  problems  that  result  fi-om  modeling  drawbar  becomes  apparent  fi-om  a 
closer  look  at  the  data  fi-om  Figure  33.  Figure  38  shows  a  closeup  view  of  part  of  the  data 
fi-om  Figure  33.  The  four  points  delineated  in  the  figure  are  fi-om  a  single  vehicle  configura¬ 
tion  with  differences  in  tire  deflection  only.  The  figure  shows  that  the  best  drawbar  numeric 
indicates  that  drawbar  should  always  increase  with  increased  deflection,  but  the  data  indicate 
what  experience  has  shown  [Jones  1992],  which  is  that  drawbar  on  high  soil  strengths 
(essentially  no  sinkage)  typically  increases  when  deflection  is  increased  from  low  magnitudes 
but  then  begins  to  decrease  as  the  deflections  get  large.  This  can  be  attributed  to  internal 
motion  resistance  of  the  tires  increasing  at  a  faster  rate  than  tractive  force  when  tire  deflec¬ 
tions  get  too  large.  This  observation  is  supported  by  Tumage’s  [1995]  discussion  of  the 
importance  of  considering  internal  motion  resistance  when  modeling  drawbar.  It  is  probable 
that  tractive  force  does  always  increase  wth  increased  tire  deflection. 

The  last  point  of  discussion  for  further  research  is  based  on  the  critical  layer  concept. 
Critical  layer  concepts  are  currently  only  applied  for  VCIj  modeling.  Traction  b^ond  VCIj 
is  currently  modeled  based  on  the  0-6  inch  soil  layer.  This  layer  is  probably  appropriate  for 
traction  throughout  the  majority  of  the  soil  strength  range  beyond  VCIj,  but  in  a  small  zone 
of  strengths  just  beyond  VCIi,  a  deeper  layer  is  probably  more  appropriate.  Further  research 
for  traction  beyond  VCIi  should  give  consideration  to  this  issue. 
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Figure  38.  Sample  of  problems  due  to  modeling  drawbar  rather  than  tractive  force 
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APPENDIX  A 

TABLES  PRESENTING  DETAILED  DESCRIPTIONS 
OF  THE  OBSERVATIONS  USED  DURING 
THE  TRACTION  MODELING 
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Table  Al.  General  description  of  the  observations  used  for  the  VCI,  modeling 
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Data 

Reference 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

impublished 

unpublished 

00 

1 

TRGL-93-8 

MP  GL-91-8 

MP  GL-91-8 

MP  GL-91-8 

Measured 

VCI, 

25 

25 

25 

25 

24 

23 

32 

30 

24 

22 

22 

ON 

20 

22 

24 

Critical 

Layer 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

3-9 

3-9 

3-9 

6-12 

6-12 

6-12 

6-12 

6-12 

GVW 

Ibf 

35070 

22310 

35770 

23385 

35060 

23585 

27100 

27100 

6300 

6300 

6300 

6300 

35780 

35780 

21325 

21325 

21325 

Tire  Characteristics 

Chains 

(Y/N) 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

x;  ^ 

o 

36 

j 

33 

35 

35 

26 

35 

VO 

fn-H 

25 

30 

lo 

30 

32 

24 

42 

20 

Inflation 

psi 

22 

24 

26 

24 

32 

20 

52 

29 

11.3 

cn 

00* 

41 

14.5 

28 

B 

22 

o 

Description 

395/80R20  Goodyear  MV/T 

- — — _ 1 

395/80R20  Goodyear  MV/T 

395/80R2P  Michelin  Pilot  XML 

395/80R20  Michelin  Pilot  XML 

395/80R20  Goodyear  MV/T 

395/80R20  Goodyear  MV/T 

Michelinll.00R16XL 

Michelin  11. 00R16  XL 

Armstrong  9.50x16.50  Bias 

Armstrong  9.50x16.50  Bias 

1200x500-508  lOPRBias 

1200x500-508  lOPRBias 

oz-oozi 

12.00-20 

or-oozi 

Vehicle 

Identification 

TCMMTV 

TCM  LMTV 

S&S  MTV 

S&S  LMTV 

TTCMTV 

TTCLMTV 

LAV  84 

LAV  84 

M1008CUCV 

M1008  CUCV 

M1008  CUCV 

M1008  CUCV 

FSU  SA-8 

FSU  SA-8 

ZIL-131 

ZIL-131 

ZIL-131 

Item 

No. 

B 

VO 

B 

00 

0\ 

o 

piH 

cs 

cn 

16 

00 

Table  Al.  (continued) 
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Data 

Reference 

MPGL-91-8  • 

MPGL-91-8  * 

CR  No.  3-152  * 

CR  No.  3-152  * 

TRGL-91-4  * 

TRGL-91-4  • 

TRGL-91-4  * 

TRGL-90-3 

TR  GL-90-3 

TRGL-90-3 

TRGL-90-3 

TRGL-90-20 

TRGL-90-1 

TRGL-89-9 

TR  GL-89-9 

MP  GL-84-6 

MP  GL-84-6 

Measured 

VCI, 

26 

«N 

ON 

ON 

CO 

CO 

ON 

CO 

CO 

cs 

irv 

cs 

29 

29 

CO 

cs 

CS 

CO 

CS 

CO 

CS 

29 

28 

Critical 

Layer 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

GVW 

Ibf 

18225 

18225 

19410 

19410 

86255 

80655 

80795 

50250 

50250 

50250 

50250 

22690 

13920 

32210 

32190 

32485 

32485 

Tire  Characteristics 

Chains 

(Y/N) 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

6lh 

percent 

cs 

29 

fTi 

CO 

«o 

CS 

es 

CO 

29 

lO 

CO 

29 

CS 

CS 

CO 

CO 

CO 

to 

CO 

cs 

Inflation 

psi 

20 

35.5 

16/22 

42 

42 

20/20 

30/30 

20/20 

30/30 

30 

19/23 

60/30 

25/25 

Description 

11.00x20 

O 

o 

o 

o 

o 

o 

o 

o 

o 

425/95R20  XZL 

17.5R25 

425/95R20  XZL 

15.00-21  TO 

15.00-21  TO 

15.00-21  TO 

15.00-21  TO 

12.5R16.5  XL 

12.5R20XL 

14.00R20  XL 

14.00R20  XL 

11.00x20  NDCC  Goodyear 

1 1.00x20  NDCC  Goodyear 

PI 

A 

K 

If 

vcnicic 

Identification 

M35A2 

M35A2 

M35A1 

M35A1 

OSHKOSH  PLS 

PACCAR  PLS 

1 

GMPLS 

TATRA  813 

TATRA  813 

TATRA  813 

TATRA  813 

FAMV 

TRM2000 

TECH  DEMO 

M923A1 

M923 

M923 

Item 

No. 

o\ 

20 

<s 

22 

CO 

CN 

25 

27 

29 

30 

CO 

CS 

CO 

CO 

CO 

CO 

36 

37 

38 

39 

Table  Al.  (continued) 
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Data 

Reference 

MPGL-84-6 

MP  GL-84-6 

MPGL-84-6  • 

MP  GL-84-6 

MP  GL-84-6 

MP  GL-84-6 

00 

1 

MP  GL-84-6 

TRGL-86-10 

TRGL-85-6 

TRGL-85-6 

TRGL-85-4 

TRGL-85-4 

TRGL-85-4  * 

TRGL-85-4 

TRGL-85-4 

unpublished 

Measured 

VCI, 

29 

22 

28 

23 

m 

<s 

00 

tn 

00 

30 

26 

32 

27 

36 

29 

27 

27 

Critical 

Layer 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

i 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

6-12 

i 

6-12 

6-12 

6-12 

GVW 

Ibf 

32260 

32260 

32105 

32105 

33585 

33585 

33790 

33790 

16170 

38620 

38620 

62250 

57770 

60145 

60145 

60145 

56290 

Tire  Characteristics 

Chains 

(Y/N) 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

L _ _ _ 

No 

6/h 

percent 

cs 

47 

22 

40 

22 

37 

rs 

40 

18*’ 

22 

40 

i _ _  .. 

29 

27 

20 

28 

34 

34 

Inflation 

psi 

30 

O 

30 

O 

00 

25 

00 

j 

44 

fO 

20/30 

20/20 

35/40 
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Table  A2.  Principal  vehicle  characteristics  associated  with  the  VCIj  observations 
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Average  per-Axle  Dimensions 
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Average  per-Axle  Dimensions 
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Table  A3.  General  description  of  the  field  observations  used  for  the  drawbar  modeling 
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0.58 

0.65 

0.63 

090 

0.71 

0.72 

0.71 

0.70 

0.65 

Z90 

00 

o 

o 

0-6 

RCI 

09 

74 

300 

300 

300 

298 

300 

300 

300 

300 

300 

300 

300 

300 

325 

228 

21 

Slip 

% 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

GVW 

Ibf 

o 

o 

ON 

1900 

7570 

7570 

8420 

7570 

8420 

8420 

8420 

8420 

6305 

6305 

6305 

6305 

28840 

o 

00 

00 

19410 

Tire  Characteristics 

Chains 

(Y/N) 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

No 

ilh 

percent 

«l 

R 

14.1 

14.1 

15/18 

16.7 

22/25 

26/31 

31/38 

33/43 

22/16 

26/20 

32/24 

34/28 

38.5 

24.9 

35 

Inflation 

psi 

I 

unknown 

20/22 

zzm 

20/30 

20/22 

14/22 

10/15 

7/10 

6/8 

14/22 

10/15 

7/10 

6/8 

15 

30 

13 

Description 

7.50-10 

7.50-10 

Bias 

Bias 

36x12.50-16.5 

Radial 

37xl2.50R16.5 

37xl2.50R16.5 

37xl2.50R16.5 

37xl2.50R16.5 

37xl2.50R16.5 

37X12.50R16.5 

37xl2.50R16.5 

37xl2.50R16.5 

12.50R20  XL 

12.50R20  XL 

o 

o 

o 

Vehicle 

Identification 

M274  Mule 

M274  Mule 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M1025  HMMWV 

M1025  HMMWV 

M1025  HMMWV 

M1025  HMMWV 

LAV  25 

LAV  25 

M35A1 

Item 

No. 

37 

00 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 
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Data 

Reference 

CR  No.  3-152 

CRNo.  3-152 

CRNo.  3-152 

CRNo.  3-152 

CRNo.  3-152 

CRNo.  3-152 

CRNo.  3-152 

TRM-70-11 

TRGL-91-2 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

unpublished 

Measurements 

1 

d 

m 

d 

© 

in 

d 

NO 

d 

r- 

<s 

d 

d 

d 

00 

in 

d 

m 

NO 

d 

NO 

© 

d 

m 

d 

d 

d 

d 

m 

d 

cn 

m 

d 

m 

d 

^  y 

®  oi 

<s 

NO 

m 

00 

cn 

00 

cn 

© 

© 

NO 

cs 

t-* 

fiM 

5 

m 

in 

© 

cs 

cn 

m 

in 

00 

00 

ON 

ON 

m 

S'? 

o 

<s 

© 

rs 

© 

© 

fS 

© 

cs 

© 

© 

CS 

cs 

cs 

© 

cs 

© 

cs 

© 

cs 

© 

cs 

S 

© 

cs 

© 

cs 

© 

cs 

s 

^  Is 

5 

o 

ON 

© 

© 

ON 

© 

On 

© 

ON 

© 

ON 

© 

ON 

m 

cs 

cs 

00 

© 

CS 

VO 

© 

00 

00 

NO 

© 

00 

00 

NO 

© 

00 

00 

NO 

S 

00 

NO 

© 

00 

00 

NO 

© 

00 

00 

VO 

s 

00 

NO 

© 

00 

00 

VO 

Tire  Characteristics 

Chains 

(Y/N) 

IZ 

:z 

'Z 

:z 

:z 

:z 

'Z 

'Z 

'Z 

"Z 

IZ 

IZ 

IZ 

IZ 

IZ 

fi/h 

percent 

»o 

C<) 

m 

m 

m 

m 

m 

m 

m 

*n 

m 

cs 

oo 

cn 

cs 

18.4/11.0 

18.4/11.0 

18.4/11.0 

18.4/11.0 

18.4/11.0 

18.4/11.0 

18.4/11.0 

© 

00 

p*4 

Inflation 

psi 

cn 

tn 

in 

in 

<n 

in 

m* 

cn 

d 

cs 

m 

cs 

NO 

NO 

VO 

VO 

NO 

NO' 

NO 

VO 

Description 

o 

o 

p 

© 

© 

© 

1 

© 

© 

© 

© 

© 

© 

© 

© 

© 

;  © 

\  © 

1  © 

11.00x20 

© 

© 

MI  24R21  XL 

NO 

cn 

cs 

NO 

cn* 

cs 

NO 

«i^ 

d 

cs 

NO 

2 

cs 

VO 

•a 

cn 

cs 

VO 

cn 

cs 

NO 

•a 

cn 

cs 

'  NO 

f-H 

i  d 
cs 

1 

• 

¥ 

vcnicic 

Identification 

< 

tn 

tn 

s 

< 

s 

< 

>o 

< 

<n 

r»i 

2 

< 

■o 

i 

i 

m 

cn 

2 

DA-1500 

Timbeijack  404 

Timbeijack  404 

Timbeijack  404 

Timbeijack  404 

Timbeijack  404 

Timbeijack  404 

Timbegack  404 

Timbeijack  404 

Item 

No. 

in 

! 

m 

m 

NO 

in 

5n 

00 

m 

ON 

m 

© 

NO 

NO 

CS 

VO 

VO 

NO 

m 

VO 

© 

© 

r- 

VO 

00 

VO 

ON 

NO 

© 

r- 
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■s 

«? 

0? 

PM 

fM 

u 

p8 

m 

cn 

1 

1 

1 

1 

1 

1 

1 

1 

• 

1 

• 

1 
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s 

s 

£ 
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162 
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CO 

m 

CO 

t- 

r- 

so 

© 

cs 

CO 

CO 

8 

V 

S 

0 

© 
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© 

00 

r- 

CO 

00 

00 

© 

VO 

© 

© 

CO 

cs 

<s 

CS 

cs 

cs 

cs 

cs 
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CS 

CO 

cs 

cs 

K 

► 
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© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

CO 

CO 

CO 

00 

00 

00 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 
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00 

r- 
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© 
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© 

0 

© 
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© 
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00 

00 
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00 

00 
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os 
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cn 

m 
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0 

0 

0 
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0 
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0 
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0 
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0 
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gb 

Z 

z 

z 

z 

z 

z 

z 

z 

z 

Z 

z 
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0 

8 
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so 

s 
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© 

© 

0 
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© 

© 

© 

© 

© 

© 

© 
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0^ 

cs 

tn 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

CS 

cs 
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Oi 

00 

•c 

8 

o 

b 

8 

© 

00 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

0\ 
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Os 

6 

c 

8 

cs 

H 

H 
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H 

H 

H 

W' 

o 
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nn 

n 

b 

H 

8 

00 

© 

00 

© 
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< 

c 
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< 

< 
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<c 

< 
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< 

o 

•0 

*9 

VO 

© 

so 

© 

VO 

SO 

© 

© 

so 

VO 

VO 

VO 

VO 

b> 

O. 

© 

© 

© 

© 

1 

1 

1 

1 

■ 

1 

1 

1 

■ 

1 

1 

I 

•c 

<o 

m 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

I-H 

fw 

X 

X 

CO 

CO 

CO 

an 

© 

© 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

0^ 

fs 

© 

© 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

cs 

CS 

00 

00 

00 

O 

cs 

cs 

8 

rj- 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

00 

00 

00 

e 

00 

00 

I 

il'-S 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

< 

< 
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© 

© 

© 

© 

© 

© 

© 

© 

© 

00 

00 

00 

1 

p 

•M  .2 

JS  c 
>  8 

8 

•c? 

0 
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03 

i 

i 

s 

i 

s 

1 

I 

i 

i 

i 

"O 

s 

tb 

b 

n 

a 

§ 

% 

§ 

NN 

P 

1 
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V>M 

cs 

CO 

»o 

VO 

00 

o\ 

© 

CS 

CO 

lO 

© 
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NN  ^ 

r- 

r- 
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r- 

r- 

r- 
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00 

00 

•00 

00 

00 

00 

-  assumed  value 
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Calculated 

CPF 

psi 

17.8 

17.8 

17.8 

00 

15.9 

15.9 

15.9 

15.9 

18.8 

18.8 

18.8 

00 

00 

801 

801 

00 

o 

801 

00 

o 

Average  per-Axle  Dimensions 

B 

o 

o 

o 

p 

66Z 

2.99 

3.06 

3.06 

00 

o 

2.08 

3.93 

3.93 

2.70 

OLZ 

2.87 

2.87 

2.87 

2.87 

2.87 

B 

14.5 

14.5 

14.5 

14.5 

roi 

roi 

roi 

roi 

13.9 

13.9 

13.9 

13.9 

8.2 

8.2 

8.2 

CS 

00 

8.2 

B 

16.1 

16.1 

16.1 

16.1 

601 

601 

601 

601 

17.1 

17.1 

17.1 

17.1 

10.2 

10.2 

ZOl 

10.2 

10.2 

B 

52.8 

52.8 

52.8 

52.8 

38.7 

38.7 

38.7 

38.7 

51.4 

51.4 

51.4 

51.4 

00 

CS 

m 

32.8 

32.8 

32.8 

32.8 

e 

cs 

cs 

<s 

CS 

fS 

CS 

CS 

cs 

CS 

CS 

CS 

CS 

^  £ 

15094 

15094 

15094 

15094 

6724 

6724 

6724 

6724 

16500 

16500 

16500 

16500 

3620 

3620 

3620 

3620 

3620 

Tire 

Construction 

0=radial 

l=bias 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

- 

- 

- 

S 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

CS 

CS 

CS 

CS 

CS 

Vehicle 

Identification 

M977  HEMTT 

M977  HEMTT 

M977  HEMTT 

M977  HEMTT 

LAV  25 

LAV  25 

LAV  25 

LAV  25 

MK48  LVS 

MK48  LVS 

MK48  LVS 

MK48  LVS 

r- 

cn 

2 

M37 

M37 

M37 

M37 

Item 

No. 

r4 

B 

NO 

B 

00 

ON 

o 

CS 

NO 

r-- 

00 
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Calculated 

CPF 

psi 

801 

10.8 

801 

801 

00 

o 

00 

o 

00 

o 

801 

00 

o 

^4 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

Average  per-Ade  Dimensions 

m 

2.87 

2.87 

2.87 

2.87 

1.07 

1.23 

1.23 

1.23 

1.23 

0.93 

0.93 

0.93 

£60 

0.93 

£60 

0.93 

0.93 

m 

8.2 

8.2 

8.2 

8.2 

00* 

8.2 

00 

8.2 

8.2 

6.2 

6.2 

6.2 

6.2 

6.2 

6.2 

6.2 

6.2 

m 

ZOl 

10.2 

10.2 

10.2 

10.2 

10.2 

10.2 

10.2 

10.2 

8.5 

8.5 

8.5 

00 

8.5 

8.5 

8.5 

8.5 

m 

32.8 

32.8 

32.8 

32.8 

32.8 

32.8 

32.8 

32.8 

32.8 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

e 

es 

tN 

es 

cs 

<s 

<s 

<s 

cs 

CS 

CS 

CS 

CS 

w 

Ibf 

3620 

3620 

3620 

3620 

3620 

3620 

3620 

3620 

3620 

950 

950 

950 

950 

950 

950 

950 

950 

Tire 

Construction 

0=radiai 

l=bias 

- 

- 

- 

- 

- 

- 

- 

- 

1—4 

E 

<s 

cs 

<s 

fS 

<s 

<s 

CS 

CS 

CS 

CS 

CS 

Vehicle 

Identification 

M37 

M37 

M37 

M37 

M37 

M37 

M37 

M37 

M37 

M274  Mule 

M274  Mule 

M274  Mule 

M274  Mule 

M274  Mule 

M274  Mule 

M274  Mule 

M274  Mule 

Item 

No. 

o\ 

20 

22 

23 

24 

25 

26 

28 

29 

30 

31 

32 

33 

34 

35 

36 
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Calculated 

CPF 

psi 

5.0 

5.0 

00 

8.3 

9.3 

8.3 

9.2 

9.2 

9.2 

9.2 

6.9 

6.9 

6.9 

69 

13.8 

13.8 

13.3 

Average  per-Axle  Dimensions 

B 

C60 

0.93 

1.27 

1.27 

1.49 

1.57 

IZZ 

2.68 

3.24 

3.57 

1.79 

9VZ 

2.63 

ON 

IN 

3.62 

2.34 

3.57 

B 

6.2 

6.2 

9.0 

9.0 

9.0 

9.4 

9.4 

ON 

On 

On 

9.4 

9.4 

9.4 

ON 

ON 

ON 

10.2 

B 

8.5 

8.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.5 

12.9 

12.9 

ZTl 

B 

22.5 

22.5 

36.3 

36.3 

36.3 

36.7 

36.7 

36.7 

36.7 

36.7 

36.7 

36.7 

36.7 

36.7 

40.5 

40.5 

43.6 

fS 

<s 

IN 

IN 

IN 

IN 

IN 

IN 

IN 

IN 

IN 

CN 

IN 

IN 

^  £ 

950 

950 

3785 

3785 

4210 

3785 

4210 

4210 

4210 

4210 

3153 

3153 

3153 

3153 

7210 

7210 

6470 

Tire 

Construction 

0==radial 

l=»bias 

- 

© 

o 

o 

o 

o 

o 

o 

o 

o 

o 

- 

S 

<s 

IN 

IN 

IN 

IN 

IN 

IN 

IN 

IN 

IN 

IN 

B 

B 

m 

Vehicle 

Identification 

M274  Mule 

M274  Mule 

M998HMMWV 

M998HMMWV 

M998HMMWV 

M998HMMWV 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M998  HMMWV 

M1025  HMMWV 

M1025  HMMWV 

M1025  HMMWV 

M1025  HMMWV 

LAV  25 

LAV  25 

M35A1 

Item 

No. 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

52 

53 
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Calculated 

CPF 

psi 

13.3 

13.3 

13.3 

13.3 

13.3 

13.3 

13.3 

VZl 

12.2 

5.8 

00 

5.8 

5.8 

5.8 

00 

5.8 

Average  per-AxIe  Dimensions 

B 

3.57 

3.57 

3.57 

3.57 

3.57 

1.53 

1.53 

2.55 

3.74 

2.50 

2.50 

2.50 

2.50 

2.50 

2.50 

2.50 

B 

10.2 

10.2 

10.2 

10.2 

10.2 

10.2 

10.2 

10.2 

15.7 

16.7 

16.7 

16.7 

16.7 

16.7 

16.7 

16.7 

B 

11.2 

11.2 

11.2 

11.2 

11.2 

11.2 

11.2 

11.2 

23.4 

23.2 

23.2 

23.2 

23.2 

23.2 

23.2 

23.2 

B 

SO 

43.6 

43.6 

43.6 

43.6 

43.6 

43.6 

43.6 

54.8 

62.5 

62.5 

62.5 

62.5 

62.5 

62.5 

62.5 

a 

■ 

«s 

CS 

fS 

<s 

CS 

6470 

6470 

6470 

6470 

6470 

6470 

6470 

6075 

15618 

8440 

8440 

8440 

8440 

8440 

8440 

8440 

Tire 

Construction 

0=radial 

l=bias 

- 

- 

- 

^4 

- 

- 

- 

o 

- 

- 

- 

- 

flpH 

S 

cn 

CO 

CO 

CO 

CO 

CO 

B 

<s 

CS 

CS 

Vehicle 

Identification 

M35A1 

M35A1 

M35A1 

M35A1 

M35A1 

M35A1 

M35A1 

M35A2 

DA-1500 

Timbeijack  404 

Timbegack  404 

Timbeqack  404 

Timbeijack  404 

Timbegack  404 

o 

•c? 

1 

H 

Timbegack  404 

Item 

No. 

54 

55 

56 

57 

58 

59 

09 

VO 

62 

63 

64 

65 

99 

67 

68 

69 

1 

Timbeijack  404  2  1  8440  2  62.5  23.2  16.7  2.50 


Table  A4.  (continued) 


-135- 


Calculated 

CPF 

psi 

5.8 

15.0 

15.0 

cS 

^4 

VZ 

VZ 

2.1 

2.1 

2.1 

2.1 

VZ 

ri 

3.2 

3.2 

3.2 

Average  per-Axle  Dimensions 

8 

in. 

2.50 

2.40 

3.20 

2.40 

2.40 

2.40 

2.40 

2.40 

2.40 

2.40 

2.40 

2.40 

OPZ 

2.80 

2.80 

2.80 

h 

in. 

16.7 

001 

001 

12.0 

12.0 

12.0 

12.0 

12.0 

o 

cs 

12.0 

12.0 

12.0 

12.0 

14.0 

14.0 

14.0 

b 

in. 

23.2 

18.3 

18.3 

o 

o 

40.0 

40.0 

o 

o 

O 

o' 

o 

o* 

O 

O 

40.0 

o 

o’ 

O 

o 

31.0 

31.0 

31.0 

d 

in. 

62.5 

43.5 

43.5 

42.0 

42.0 

42.0 

42.0 

42.0 

42.0 

42.0 

42.0 

42.0 

42.0 

48.0 

o 

oo' 

48.0 

e 

<s 

<s 

cs 

CS 

8440 

11927 

11927 

3606 

3606 

i _ 

3606 

3606 

3606 

3606 

3606 

3606 

3606 

3606 

4753 

4753 

4753 

Tire 

Construction 

0-radial 

l=bias 

- 

- 

- 

I'M 

- 

S 

Vehicle 

Identification 

Timbe'ijack  404 

FSU  SA-8 

FSU  SA-8 

MEXA  10x10 

MEXA  10x10 

MEXA  10x10 
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Table  A5.  General  description  of  the  laboratory  observations  used  for  the  drawbar  modeling  and  the  associated  principal  tire 
characteristics 
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Table  Bl.  Quality  comparison  for  all  VCIj  model  variants 


Model 

Model 

Degrees  of 

Standard 

Adjusted 

No. 

Description 

Freedom 

Error  (psi) 

R-squared 

EXISTING  COMPARISON  MODELS 


Cl* 

Standard  NRMM  Model:  Mobility  Index  with  Deflection 

2.68 

90.3 

Correction  Factor  (MI&DCF) 

C2* 

Rowland’s  Mean  Maximum  Pressure  (MMPr  ) 

73 

2.91 

88.5 

C3* 

Lanninie’s  Mean  Maximum  Pressure  (MMP^ ) 

73 

2.96 

88.1 

C4* 

Maclaurin’s  Mean  Maximum  Pressure  (MMPm  ) 

74 

3.01 

87.7 

PRELIMINARY  NEW  DEVELOPMENT  MODELS 


PI 

CPF 

77 

3.60 

82.4 

P2 

|CPF*DCF 

76 

2.79 

89.4 

P3 

Hard  Surface  Contact  Pressure  Parameter  (CPH) 

77 

2.91 

88.5 

P4 

Soft  Soil  Contact  Pressure  Parameter  (CPS) 

77 

■  3.50 

83.3 

P5 

Modified  CPH  version  1:  Replaced  d  with  3.876*h  wiiich 
originates  from  the  correlation  between  d  and  h 

76 

3.29 

85.3 

P6 

Modified  CPH  version  2:  Assumed  6^  is  negligible  in  the 
contact  length  estimator 

77 

2.99 

87.8 

P7* 

VCI 

Freitag’s  Clay  Numeric  cast  in  a  VCI  form 

16 

2.82 

89.2 

P8 

Tumage’s  Clay  Numeric  cast  in  a  VCI  form 

16 

2.94 

88.2 

P9 

Smith’s  Clay  Numeric  cast  in  a  VCI  form  ( 

75 

2.93 

88.4 

PIO 

Maclaurin’s  Motion  Resistance  Numeric  cast  in  a  VCI  form 

74 

3.77 

80.7 

Pll 

Average  Contact  Pressure  Parameter  No.  1  (CPI ) 

77 

2.66 

90.4 

P12 

Modified  CPI  version  Oa:  Used  d/2  instead  of  i4  LC^j 

77 

2.68 

90.3 

P13 

Modified  CPI  version  Ob:  Used  +  h^)  -  (6  +  instead 

of'ALC^, 

77 

2.66 

90.4 

PRINCIPAL  NEW  DEVELOPMENT  MODELS 


la 

Modified  MI&DCF  Model:  Optimized  Curve  Fit  to  Database 

70, 69, 

68,  or  67 

2.52 

91.4 

lb 

Modified  MI&DCF  Model:  Used  CPF/10  instead  of  WF 

69 

2.80 

89.4 
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Table  Bl.  (continued) 


Model 

No. 

Model 

Description 

Degrees  of 
Freedom 

Standard 
Error  (psi) 

Adjusted 

R-$quared 

Ic 

Modified  MI&DCF  Model:  Used  CPF/10  instead  of  WF  and 
left  out  WLP 

70 

2.90 

88.6 

Id 

Modified  MI&DCF  Model:  Partially  Optimized  to  Database 
(Curve  Fits  and  WF  only) 

2.48 

91.7 

le* 

Modified  MI&DCF  Model:  Fully  Optimized  to  Database  (GF 
and  EF  have  no  effect) 

70,69, 

68,  or  67 

2.44 

91.9 

2 

CPI  and  m 

75 

2.65 

90.4 

3 

Modified  CPI  version  1:  Empirical  adjustments  to  the  3  parts 
of  the  contact  length  using  3  constants 

14 

2.69 

90.2 

4 

Modified  CPI  version  2:  Empirical  adjustment  to  the  hard- 
surface  deflection  using  a  constant 

16 

2.67 

90.3 

j 

5 

CPI  with  modified  TF  (product) 

16 

2.64 

90.5 

6 

CPI  and  w 

75 

2.69 

90.2 

7 

CPI  and  Aspect  Ratio  (h^) 

75 

2.54 

91.2 

8 

Modified  CPI  version  3:  Empirical  adjustment  to  the  contact 
area  using  aspect  ratio 

76 

2.50 

91.5 

9 

CPI  y  with  modified  TF  (product) 

75 

2.46 

91.8 

10 

Modified  CPI  version  4:  Empirical  adjustment  to  the  contact 
area  using  d/b 

16 

2.52 

91.3 

11 

Modified  CPI  version  5;  Empirical  adjustment  to  the  contact 
area  using  h/d 

16 

2.66  ^ 

90.4 

12a 

Modified  CPI  version  6a:  Empirical  adjustment  to  the  contact 
area  using  h/b,  d/b,  and  h/d 

74 

2.52 

91.4 

12b 

- — - ;  j 

Modified  CPI  version  6b;  Empirical  adjustment  to  the  contact 
area  using  h/b  and  d/b 

75 

2.50 

91.5 

13 

Modified  CPI  version  7:  Empirical  adjustment  to  the  contact 
area  using  a  constant  and  h/b 

16 

2.50 

91.5 

14 

Modified  CPI  version  8:  Empirical  adjustments  to  (1)  the 
additional  contact  length  due  to  sinkage  using  a  constant  and 
(2)  the  contact  area  using  aspect  ratio 

75 

2.47 

91.7 

15 

Modified  CPI  version  9:  Empirical  adjustments  to  (1)  V2LC^j 
using  a  constant  and  (2)  using  aspect  ratio 

75 

2.48 

91.7 
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Table  Bl.  (continued) 


Model 

No. 

Model 

Description 

Degrees  of 
Freedom 

Standard 
Error  (psi) 

Adjusted 

R-squared 

16 

Modified  CPI  version  10:  Empirical  adjustments  to  (1)  the 
additional  contact  length  due  to  sinkage  using  a  constant  and 

(2)  5  -  5^  using  aspect  ratio 

75 

2.47 

91.7 

17 

Modified  CPI  version  11:  Empirical  adjustments  to  (1) 

using  a  constant,  (2)  the  additional  contact  length 
due  to  sinkage  using  a  constant,  and  (3)  the  contact  area  using 
aspect  ratio 

74 

2.49 

91.6 

18 

Modified  CPI  version  12:  Removed  1/m  and  made  empirical 
adjustments  \o{\)l^db  -  5^  using  a  constant,  (2)  the 
additional  contact  length  due  to  sinkage  using  a  constant,  and 
(3)  the  contact  area  using  aspect  ratio 

74 

2.50 

91.5 

19 

Modified  CPI  version  13:  Empirical  adjustments  to  (1)  the 
two  parts  of  the  additional  contact  length  due  to  sinkage  using 
constants  and  (2)  the  contact  area  using  aspect  ratio 

74 

2.48 

91,6 

20 

CPi  g  and  Tire  Construction  Factor  (TCF) 

74 

2.40 

92.1 

21 

Modified  CPI  version  14:  Empirical  adjustments  to  (1)  the 
additional  contact  length  due  to  sinkage  using  a  constant  and 
(2)  the  contact  area  using  aspect  ratio  and  TCF 

74 

2.42 

92.1 

22 

CPI  g  with  modified  TF  (product)  and  TCF 

73 

2.42 

92.1 

23 

CPI  g  and  modified  TCF  [Form  1  ] 

73 

2.41 

92.1 

24 

CPI  g  and  modified  TCF  [Form  2] 

73 

2.39 

92.3 

25 

CPi'  andh^ 

73 

2.50 

91.5 

26 

CPH  and  modified  CPS  version  1:  Included  1/m  in  CPS 

76 

2.62 

90.7 

27 

Modified  CPH  version  3  and  modified  CPS  version  2:  Included 
1/m  in  CPS  and  made  empirical  adjustment  to  both  hard 
surface  and  soft  soil  parameter  contact  areas  using  aspect  ratio 

74 

2.47 

91.7 

28 

Modified  CPH  version  3  and  modified  CPS  version  3: 

Empirical  adjustment  to  both  hard  surface  and  soft  soil 
parameter  contact  areas  using  aspect  ratio 

74 

2.51 

91.4 

29 

Numeric  version  1:  MMP  style  numeric  with  constraint  for 
pressure  units 

75 

2.65 

90.5 

30 

Numeric  version  2:  MMP  style  numeric  without  the  constraint 
for  pressure  units 

74 

2.53 

91.3 
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Table  Bl.  (continued) 


Model 

No. 

Model 

Description 

Degrees  of 
Freedom 

Standard 
Error  (psi) 

Adjusted 

R-squared 

31a 

Numeric  version  3a:  Numeric  using  bd  and  4  dimensionless 
ratios  of  tire  variables  [b/d,  d/h,  6/d,  and  h/b]  viiich  allows 
inclusion  of  b,  d,  h,  and  6  with  optimal  exponents  under  the 
constraint  of  pressure  units 

73 

2.53 

91.3 

31b 

Numeric  version  3b:  Numeric  using  bd  and  3  dimensionless 
ratios  of  tire  variables  [b/d,  d/h,  and  6/d]  which  allows 
inclusion  of  b,  d,  h,  and  6  with  optimal  exponents  under  the 
constraint  of  pressure  units 

74 

2.51 

91.4 

32 

Modified  CPI  version  15  and  TCP:  Empirical  adjustments  to 
(1)  Vd  6  -  6^  using  a  constant  and  (2)  the  contact  area  using 
aspect  ratio 

74 

2.39 

92.2 

33 

Modified  CPI  version  16  and  TCP;  Empirical  adjustments  to 
(1)  2^jdb  -  6^  using  a  constant,  (2)  the  6^  part  of  the 

additional  contact  length  due  to  sinkage  using  a  constant,  and 
(3)  the  contact  area  using  aspect  ratio 

73 

2.41 

92.1 

34 

Modified  CPI  version  17  and  TCP:  Empirical  adjustments  to 
j  (1)  6  -  6^  using  a  constant,  (2)  the  additional  contact 

length  due  to  sinkage  using  a  constant,  and  (3)  the  contact  area 
using  aspect  ratio 

73 

2.41 

92.1 

35 

Modified  CPI  version  18  and  TCP:  Empirical  adjustments  to 
(1)  V46  -  6*  using  a  constant,  (2)  '/iLC,,  using  a  constant, 
and  (3)  the  contact  area  using  aspect  ratio 

73 

2.41 

92.1 

36* 

Modified  CPI  version  19  and  TCP:  Empirical  adjustments  to 
(1 )  Vi  using  a  constant  and  (2)  the  contact  area  using 

aspect  ratio 

74 

2.39 

92.2 

37 

Modified  CPI  version  20  and  TCP:  Empirical  adjustments  to 
(1)  the  additional  contact  length  due  to  sinkage  using  a 
constant,  (2)  'A  LC^,  using  a  constant,  and  (3)  the  contact  area 
using  aspect  ratio 

73 

2.40 

92.2 

38 

Modified  CPI  version  21  and  TCP:  Replaced  Vi  LC^j  with 
c^dfl  and  made  an  empirical  adjustment  to  the  cont^t  area 
using  aspect  ratio 

74 

2.39 

92.2 

39 

Modified  CPI  version  22  and  TCP:  Replaced  ViLC^j  with 
ViLCy3  (i.e.,  changed  z  from  to  c^d)  and  made  an 

empirical  adjustment  to  the  contact  area  using  aspect  ratio 

74 

2.40 

92.2 

40 

Modified  CPI  version  23  and  TCP:  Replaced  ViLC^;  with 

Vi (i.e.,  changed  z  fiom  c^h)  and  made  an 

empirical  adjustment  to  the  contact  area  using  aspect  ratio 

74 

2.40 

92.2 
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Table  Bl.  (continued) 


Model 

No. 

Model 

Description 

Degrees  of 
Freedom 

Standard 
Error  (psi) 

Adjusted 

R-squared 

41 

Modified  CPI  version  24  and  TCP:  Replaced  with 

(i.e.,  changed  z  from  to  VaJ)  and  made  empirical 

adjustments  to  (I)  the  additional  contact  length  due  to  sinkage 
using  a  constant  and  (2)  the  contact  area  using  aspect  ratio 

73 

2.42 

92.0 

42 

Modified  CPI  version  25  and  TCP:  Replaced  ViLC,i  with 
(i.e.,  changed  z  fix>m  to  c^d)  and  made  empirical 

adjustments  to  (1)  ^db-b^  using  a  constant  and  (2)  the 
contact  area  using  aspect  ratio 

73 

2.41 

92.1 

43 

Modified  CPI  version  26  and  TCP:  Replaced  '/tLC,;  with 
ViLC^  (i.e.,  changed  z  fiom  to  c,</)  and  made  an 

empirical  adjustment  to  y/dd  -6^  using  aspect  ratio 

74 

2.40 

92.2 

44 

Modified  CPI  version  27  and  TCP:  Replaced  with 

54 LC^  (i.e.,  changed  z  fiom  to  c^d)  and  made  empirical 

adjustments  to  (1)  the  additional  contact  length  due  to  sinkage 
using  a  constant  and  (2)  ^d6-6^  using  aspect  ratio 

73 

2.40 

92.1 

45 

Modified  CPI  version  28  and  TCP;  Removed  1/m,  replaced 

14  with  14LC^5  (i.e.,  changed  z  fiom  to  c,</),  and 

made  empirical  adjustments  to  (1)  the  additional  contact 
length  due  to  sinkage  using  a  constant  and  (2)  )/db-b^  using 
aspect  ratio 

73 

2.40 

92.2 

46 

CRff'andTCP 

75 

2.62 

90.7 

47 

CP5j  and  TCP 

75 

3.15 

86.5 

48* 

1 

1 

Modified  Average  Contact  Pressure  Parameter  No.  2  (CP2) 
version  1  and  TCP:  Used  c^d  for  internal  length  adjustment 
and  made  empirical  adjustment  to  the  contact  area  using 
aspect  ratio 

74 

2.37 

92.4 

49 

Modified  CP2  version  2  and  TCP:  Split  the  contact  length  into 

2  parts  and  simplified  the  adjustment/correction  part  to  Cj  d 

74 

2.39 

92.3 

50 

Recast  CP2  [  and  TCP:  Split  the  contact  length  into  2  parts 

74 

2.37 

92.4 

51* 

A^^andTCP 

73 

2.40 

92.2 

^  -  These  models  were  discussed  in  the  main  text 
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VCIi  Model  Cl 

Standard  NRMM  Model:  Mobility  Index  with  Deflection  Correction  Factor  (MI&DCF) 


CPFWF 

TEFGF 


where  CPF  = 


+  WLF  -  CF  (EF)(TF) 


0.5  ndb 

WF  =  c^,  (w/1000)  + 
where  w<  2,000  i 


TEF  = 


w<  2,000/6/ 

^WFl 

=  0.553 

and 

CWF2  =  0 

2,000  ^w<  13,500/6/ 

^WFl 

=  0.033 

and 

=  1050 

13,500  :£w<  20,000/6/ 

^WFJ 

=  0.142 

and 

=  -0.420 

20,000  ^  w 

^  ^WFl 

=  0.278 

and 

^WF2  -  -3- 115 

10  +  6 

GF  =  1  +  0.05  Cqp  where  c^p  =  1  if  tire  chains  are  used  or  0  if  not 


WLF  = 


K 

CF  =  — 
10 


EF  =  I  +  0.05  Cpp  where  c^p  =  1  if  PWR  <  10  hptton  or  0  if  not 

7F  =  1  +  0.05  Cjp  where  c^p  =  \  if  manual  transmission  or  0  if  automatic 


:.MI  = 


i  w  ^wFi 

[  0.5 /irfft(0.10  + 0.016)  1  +  0.05  c^p 

U/M 


■WF2  ^  W 


2000  10 


^1+0.05  C£pj|l+b.05Cjj,j 

/ 


ycij  =  f  {MI,  DCF) 


where  115  -  /(M, DCF)  =  |ll.48 +0.2Aff--j^^^jz)CF 

MI>115  -  /(A47,DCF)  =  |4.1M'’‘*^jDCF 
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VCIj  Model  C2 

Rowland’s  Mean  Maximum  Pressure  (MMP^) 


MMPg^ 


0.91  KGVW 


where  for  all-drive  vehicles 


m=2  =*  K  =  3.66 
ot=4  -  a:  =  4.io 


m  =  3  ->  a:  =  3.90 
ot  =  5  ■=*  A:  =  4.32 


FC/,  =  3.65  +  O.AllMMPu 


veil  Model  C3  ^ 

Larminie’s  Mean  Maximum  Pressure  (MMPl) 


MMP^  = 


0.91K'GVW 
d'^^  ^fbid 


where  for  all-drive  vehicles  m  -  2  =*  K'  =  1.83 

m=A  =*  K'  =  2.05 


OT  =  3  -»  a:'  =  1.95 
TO  =  5  =»  i:'  =  2.16 


FC/,  =  3.68  +  O.AiOMMP^^ 


veil  Model  CA 

Maclaurin’s  Mean  Maximum  Pressure  (MMPm) 

MMP,- _ . 

FC/,  =  2.53  +  1.35  MMPj^ 

veil  Model  PI 
ePF 


CPF 


'Andb 


FC/,  =  4.61  +  1.45  CPF 
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veil  Model  P2 
CPF*DCF 


CPF*DCF  = 


Yindb  [  6/h  j 


CPF*DCF  = 


1.245  m' 
ndbiblhf^ 


VCI^  =  3.14  +  1.19  {CPF*DCF) 


veil  Model  P3 

Hard  Surface  eontact  Pressure  Parameter  (ePH) 

CPH  = - T= 

Inb^idb-b'^ 

VCl^  =  3.12  +  1.53  CPH 


veil  Model  P4 

Soft  Soil  eontact  Pressure  Parameter  (ePS) 

CPS  =  — 

n  b  LC^j 

where  -  b  ^  LC^j  =  2  ^d(b  +  h^)  -  (6  + 

h  -  b  ^ 

C  2 

yCIj  =  4.83  +  2.75  CPS 


-152- 


VCIi  Model  P5 

Modified  CPH  version  1;  Replaced  d  with  3.876*h  which  originates  fi’om  the  correlation 
between  d  and  h 

_ w _ 

2nbyj3.S76h5- 

3.67  +  1MCPH{ 


CPH[  = 
FC/,  = 


veil  Model  P6 

Modified  CPH  version  2:  Assumed  6^  is  negligible  in  the  contact  length  estimator 

CP/f '  = - ^ - 

Inb^dW 

FC/j  =  3.26  +  1.57  CP/T' 


veil  Model  P7 

Freitag’s  Clay  Numeric  cast  in  a  VCI  form 


N. 


VCI 

c(F) 


W 

nbdyfm 


VCJj  =  2.99  +  1.58 
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VCIi  Model  P8 

Tumage’s  Clay  Numeric  cast  in  a  VCI  form 


FC/,  =  2.44  +  lAON^^ 

veil  Model  P9 

Smith’s  Clay  Numeric  cast  in  a  VCI  form 


VCJj  =  3.09  +  4.09 

veil  Model  PIO 

Maclaurin’s  Motion  Resistance  Numeric  cast  in  a  VCI  form  (Nj^  ) 

J,VCI  _  _ ^ _ 

FC/,  =  2.30  +  1.26 


VCI,  Model  PI  1 

Average  Contact  Pressure  Parameter  No.  1  (CPI) 


FC/,  =  4.19  +  1.70  CP7 
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VCIi  Model  P12 

Modified  CPI  version  Oa:  Used  d/2  instead  of  ViLC^j 

Cpj!  =  _ _ 

nb^l^dd-b^  +  l//«(d/2 -^48  -  6^)] 

KC/,  =  4.13  +  1.73  CP7^ 


veil  Model  P13 

Modified  CPI  version  Ob:  Used  ^ci(b  +  h^)  -  (b  +  hj^  instead  of 

n6|2/(i6-6^  +  l/m\^^d(b  + h^)  -  (5  + -^db-b^ 
VCI^  =  4.18  +  l.lOCPl^ 


veil  Model  la 

Modified  MI&DCF  Model:  Optimized  Curve  Fit  to  Database 


Ml 


w 


^wFi  (w/1000)  +  c. 


WF2 


W 


0.5  w  4ft  (0.10 +  0.016) 


1  +  0.05  c, 


GF 


2000 


I  ^  1+0.05  1+0.05  Cjp) 


DCF'  =  f  % 
^  \  6 


is'' 

bih  j 


0.29 


FC/,  =  fiM,DCF'j 

where  A77  <:  1 15  -»  fiMI ,  DCF'J  =  |  9.95  +  0.232  MI  -  ■  ]  DCF' 

^  MI  +  4.49  ) 

MI>  115  -*  fiMI,DCF'^)  =  ^7.80  MI  DCF' 
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VCI,  Model  Id 

Modified  MI&DCF  Model:  Partially  Optimized  to  Database  (WF  only) 


M'  = 


CPFWF' 


^  +  WLF  -  CF 


TEF  GF 

=  ^WFt  (>^/1000)  + 

where  w<  2,33S  Ibf 

2,338  ^w<  12,564/6/ 
12,564  ^  w 


(EF)  (TF) 


^WFI  ~ 

^WFl  ~  0-0^ 

Cjfpj  =  0.093  and 


w 


c^j{wiim)^c. 


WF2  ^ 


^WF2 

c^=  0.988 
^WF2  ~  0-372 


DCF'-(9^Y 

“  (6/hj 


0.5  nrf6  (0.10 +  0.016) 

0.31 


1  +  0.05  c, 


GF 


2000  10 


|l+ 0.05  C£p.j^l+ 0.05  Cjyj 


VCI,  =  JiMl'^,DCF'a) 

where  MI 'j  ^  ns  -  /(A47',  DCF' )  =  f  10.2  +  0.223  - |DCF' 

I  +  2.62  I 

A47'  >  115  -+  /(Ai7',DCF')  =  |3.21  DCF' 
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VCIi  Model  le 

Modified  MI&DCF  Model:  Fully  Optimized  to  Database  (GF  and  EF  have  no  effect) 


^  CPW  WF^  ^ 

V  =  f-l  +  WLF'  -  CF'  (£F)(7F') 

*  TEF'  GF 

\  ^  / 


WF'^  = 
where 


,  (w/lOOO) 


TEF'  = 


WLF',  = 


,  K 

CF',  = 

*  6.9 


w<  2,338/6/  -» 

^WFl 

=  0.467 

and 

^WF2  “ 

2,338  ^w<  12,564/6/  - 

^WFl 

=  0.044 

and 

^WF2  “ 

12,564  iw  -* 

^WFl 

=  0.093 

and 

^WF2  “ 

8.3  +  1.136 

TF',  =  1  +  0.10  Cjy  where  =  1  manual  transmission  or  0  if  automatic 

_  ( _ ^ ^  _2!:_  _  A]fl+0.05cJfl+0.10cjy] 

0.5  nrfft  (0.083 +0.0113Z»)  I  +  O.OSCq^  2000  6.9  jV  )\  ) 


•  ui' -\ _ - _ ^ +  - l+0.U5c~ 

‘  *  [0.5  nrfft  (0.083 +0.0113Z>)  1  +  0.05Cgp  2000  6.9  jV  ) 

VCf  =  f{Ml'„DCF',) 

where  Ml'  ^115  -•  f(Ml',,DCF',)  =  f  10.2  +  0.222 A/7' - — —  IeCF^ 

I  Ail^  +  2.62  I 

M',  >115  -*  fiM'^DCF',)  =  ^3.216  DCF' 


veil  Model  2 
CPI  and  m 

VCf  =4.19  +  \MCP1  +  0.000166  m® 
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VCI,  Model  9 

CPij  vdth  modified  TF  (product) 
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TF^j,  =  1  +  C|  Cjj-  where  Cjp  =  \  if  manual  transmission  or  0  if  automatic 


Numerical  Optimization  yields  /?,  =  -0.47  and  Cj  =  0.06 


VCIj  Model  10 

Modified  CPI  version  4:  Empirical  adjustment  to  the  contact  area  using  d^ 


Numerical  Optimization  yields  /?,  =  -0.49 

CFi' - . . , - 

n b  \2)/d6  -  6^  +  1/m  ( '/i LC^,  -  ^jd^  -  6^  )  J (b  Idf  ^ 

FC7,  =  5.09  +  0.917  CP/' 
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VCIi  Model  13 

Modified  CPI  version  7:  Empirical  adjustment  to  the  contact  area  using  a  constant  and  h^ 
CPli,  = 

Numerical  Optimization  yields  c,  =  1  and  =  -0.46 
CPli,  = 

VCI^  =5.10  +  1.79  CP7' 
veil  Model  14 

Modified  CPI  version  8:  Empirical  adjustments  to  (1)  the  additional  contact  length  due  to 
sinkage  using  a  constant  and  (2)  the  contact  area  using  aspect  ratio 

CPI',  = 

Numerical  Optimization  yields  c,  =  0.59  and  =  -0.50 
CPI',  = 

FC7,  =  4.76  +  1.73  CP7' 


nc,b  2^db  -  fp-  +  llm[v2LC^,-yjdb  -  (hlb)’’' 
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VCIi  Model  15 

Modified  CPI  version  9:  Empirical  adjustments  to  (1)  V^LC^,  using  a  constant  and 
(2)  yjdb-b'^  using  aspect  ratio 


VCI,  Model  16  ,  .  j 

Modified  CPI  version  10;  Empirical  adjustments  to  (1)  the  additional  contact  length  due 

to  sinkage  using  a  constant  and  (2)  ^db  -  b^  using  aspect  ratio 


Numerical  Optimization  yields  c,  =  0.69  and  /»,  -  -0.62 


~  nb\2>ldb  -^ +  Qmim[viLC^t  - ^jdb  - 
VCI^  =  4.68  +  1.76CP7(o 
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VCI,  Model  17 

Modified  CPI  version  11:  Empirical  adjustments  to  (1)  l^db  -  6^  using  a  constant,  (2) 
the  additional  contact  length  due  to  sinkage  using  a  constant,  and  (3)  the  contact  area 
using  aspect  ratio 

CPl[,  = 

Numerical  Optimization  yields  c,  =  1.2,  Cj  =  0.74,  and  p,  =  -0.50 
CPI[^  = 

FC/,  =  4.79  +  2.08  CP7', 


nb  2c^)/d6-S^  +  c^illm)[v2LC^j  -  ^db  -  6^ )  (b/bf' 


veil  Model  18 

Modified  CPI  version  12:  Removed  1/m  and  made  empirical  adjustments  to  (1) 
2^d6  -  6^  using  a  constant,  (2)  the  additional  contact  length  due  to  sinkage  using  a 
constant,  and  (3)  the  contact  area  using  aspect  ratio 


W 
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VCIi  Model  19 

Modified  CPI  version  13:  Empirical  adjustments  to  (1)  the  two  parts  of  the  additional 
contact  length  due  to  sinkage  using  constants  and  (2)  the  contact  area  using  aspect  ratio 


Numerical  Optimization  yields  c,  =  0.62,  =  0.72,  and  /?,  =  -0.51 


VC1^  =  4.74  +  1.71  CPy  (3 


VCI,  Model  20 

CPI  g  and  Tire  Construction  Factor  (TCF) 


TCF  =  where  Cj.(,p  =  0  if  radial  construction  or  1  if  bias  ply  construction 

Numerical  Optimization  yields  c,  =  0.71  and  />,  =  -0.45 


VCf  =  3.09  +  IXICP1[  +  1.50  rCF 
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VCIi  Model  21 

Modified  CPI  version  14:  Empirical  adjustments  to  (1)  the  additional  contact  length  due 
to  sinkage  using  a  constant  and  (2)  the  contact  area  using  aspect  ratio  and  TCP 


CP1\^  =  — , —  . . . .  p. - 

nb\2^d6-6^  +  c,(l/ifi)('/2£q, -)/d6  -  6^  ) ](*/*)'■  (1  *c^TCF) 
Numerical  Optimization  yields  Cj  =  0.72,  =  -0.47,  and  =  -0.06 


CP1[^  = 


n b  2^1  db  -  6^  +  O.n/miViLC^,  - ^46  -  6^  )  (b Ihf^ (1-0.06 TCP) 


VCI^  =  4.52  +  1.72  CP7 '4 


veil  Model  22 

CPig  wth  modified  TF  (product)  and  TCP 


TFIj22  =  1  +  C2  ^TF  y»here  c^p  =  1  1/  manual  transmission  or  0  if  automatic 
Numerical  Optimization  yields  Cj  =  0.72,  p^  =  -0.46,  and  C2  =  0.01 


VCf  =3.18  +  1.8l(CP7'rFj;^^)  +  \39TCF 


VCI,  Model  23 

CPig  and  modified  TCF  [Form  1] 
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cp;'  = 


n b |27£/6  -  6^  +  Cj(l//w) [^ALC^j  -  -  6^  j  j (hlb)^^ 

ia$  =  [^TCF  ®  construction  or  1  if  bias  ply 


Numerical  Optimization  yields  Cj  =  0.68,  p^  -  -0.50,  and  pj  -  -0.27 


CP7'  = 


*  nb\l^dh  -  62  +  QMIm[v2LC^,-y]db  -  ^)\(b/h)° 


KC/,  =  1.33  +  1.85  CP7'  +  139  TCF'^ 


veil  Model  24 

CPI  g  and  modified  TCF  [Form  2] 


CPI',  = 


rr  _ 

nb\2y]dh  -  6^  +  c,(l/m) -\/<76  -  6^ ) J (A /i) 

TCF'u,j  =  ^Cj.(.p  +  ^TCF  ~  ®  if  radial  construction  or  1  if  bias  ply 


M24  "  I  ^TCF 


Numerical  Optimization  yields  c,  =  0.78,  p■^  =  -0.60,  and  P2  -  ^ 


CPI',  = 


nb\2}Jd5  -  6^  +  0.78/»i(>/2lQ  - ^48  -  6^  )](i/7j)®® 
^CF'h^4  =  +  lj|— j 


KC7,  =  4.85  +  1.88  CP7'  -  4.13  TCF'^ 
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VCIi  Model  25 
CP7'  andh/b 


CP7' 


w 


nb\2y/d6  -  6^  +  Cj(l/m)[y2LC^, -  8^  )](A/A)^> 
Numerical  Optimization  yields  Cj  =  0.60  a/irf  =  -0,45 


•••  CP/'  = 


w 


nb  lyjdb  -  6^  +  0.60/ffi ['ALC^,  -\ldb-b^) 


(b/h) 


0.45 


FC/,  =  4.64  +  1.70 CP/'  +  0.712 1  j 


6.4 


veil  Model  26 

CPH  and  modified  CPS  version  1;  Included  1/m  in  CPS 


CPH - — . - 

Inb^db  -  b^ 

Qpgf  -  _ _ 

’  "  nb(,llm)LC^, 

VCI^  =  3.97  +  1.16  CP^  +  0.190  CPS'( 
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VCI,  Model  27 

Modified  CPH  version  3  and  modified  CPS  version  2:  Included  1/m  in  CPS  and  made 
empirical  adjustment  to  both  hard  surface  and  soft  soil  parameter  contact  areas  using 
aspect  ratio 


veil  Model  28 

Modified  CPH  version  3  and  modified  CPS  version  3;  Empirical  adjustment  to  both  hard 
surface  and  soft  soil  parameter  contact  areas  using  aspect  ratio 


veil  Model  29 

Numeric  version  1:  MMP  style  numeric  with  constrjunt  for  pressure  units 


Numerical  Optimization  yields  =  1.45  and  p^  =  0.22 

•  ^  ^ _ 

'  '  451/0.22  50.33 

FC7,  =  4.06  +  0.717 


VCI,  Model  30 

Numeric  version  2:  MMP  style  numeric  without  the  constraint  for  pressure  units 
N  =  ^ 

Numerical  Optimization  yields  p^  =  1.57,  p^  =  -0.32,  and  p^  =  0.41 
Note:  Pi  +/?2  +P3  =  1*66 

•  N  - 

^  n  b  6®'" 

FC/,  =  5.40  +  O.I28JV2 


veil  Model  31a 

Numeric  version  3a:  Numeric  using  bd  and  4  dimensionless  ratios  of  tire  variables  [b/d, 
d/h,  6/d,  and  h/b]  which  allows  inclusion  of  b,  d,  h,  and  6  with  optimal  exponents  under 
the  constraint  of  pressure  units 


nbd  (b/d)”'  (d/h)”^ (b/d)”^ (h/b)”* 

Numerical  Optimization  yields  /?,  =  0.11,  =  0.28,  p^  =  0.36  and  p^  =  -0.39 
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VCIi  Model  31b 

Numeric  version  3b:  Numeric  using  bd  and  3  dimensionless  ratios  of  tire  variables  [b/d, 
d/h,  and  6/d]  which  allows  inclusion  of  b,  d,  h,  and  6  with  optimal  exponents  under  the 
constraint  of  pressure  units 


w 


nbdibld)’’'(dlhf^6/d/^ 

Numerical  Optimization  yields  /7j  =  0.50,  =  0.67  and  p^  =  0.36 


n  b 


1.50  j  0.81  sO.36 


FCT,  =  4.08  +  \.56N 


3b 


VCI,  Model  32 

Modified  CPI  version  15  and  TCP:  Empirical  adjustments  to  (1)  \/db-5^  using  a 
constant  and  (2)  the  contact  area  using  aspect  ratio 


CPI  is  =  - 1 .  . — - ; -  -  n - 

62  +  (l/ni)('/2£C,;  - c,^/rf6  -  S^jp/bf' 

Numerical  Optimization  yields  c,  =  1.4  and  p^  =  -0.47 

.  rpi'  =  _ - _ _ 

nb[2.&sfd6^  +  {l/m){y2LC^,-lA^/d5-5^)](blh)<>^^ 

yCI,  =  2.84  +  2.47  CP;  (5  +  1.58  rcF 
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VCIi  Model  33 

Modified  CPI  version  16  and  TCP:  Empirical  adjustments  to  (1)  lyjdb  -  6^  using  a 

constant,  (2)  the  yjdb-b^  part  of  the  additional  contact  length  due  to  sinkage  using  a 
constant,  and  (3)  the  contact  area  using  aspect  ratio 

Numerical  Optimization  yields  c,  =  1.2,  Cj  =  1.6  and  p,  =  -0.49 
•••  CP1[,  = 

PC/,  =  2.71  +  2.10  CP/ (g  +  1.71  rCF 


w 


veil  Model  34 

Modified  CPI  version  17  and  TCP:  Empirical  adjustments  to  (1)  ^db-b^  using  a 
constant,  (2)  the  additional  contact  length  due  to  sinkage  using  a  constant,  and  (3)  the 
contact  area  using  aspect  ratio 
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VCIi  Model  35 

Modified  CPI  version  18  and  TCP:  Empirical  adjustments  to  (1)  using  a 

constant,  (2)  ViLC^j  using  a  constant,  and  (3)  the  contact  area  using  aspect  ratio 


(^pj  •  =  _  I — — - 

'*  nb^c^'fdb-  +  \lm[c^V2LC^J-c^^|db- 
Numerical  Optimization  yields  c,  =  1.1,  Cj  =  0.83  and  p,  =  -0.47 

,  w  _ 

'*  „b  2.2^ db  -  62  +  1/OT  0.415 IC^,-l.l^/</6-  82)  (6//i)' 


PC/,  =  2.90  +  1.96  CP/ {g  +  1.59  TCP 


veil  Model  36 

Modified  CPI  version  19  and  TCP;  Empirical  adjustments  to  (1)  YzLC^j  using  a  constant 
and  (2)  the  contact  area  using  aspect  ratio 


(2p]  f  =  -  — _  ^  ^ - 

'  nb\2{d^  +  -yfdb^)\ihlbY' 

Numerical  Optimization  yields  c^  =  0.74  and  p,  =  -0.47 


A  CP/jg  = 


n b  2yjdb  -  b^  +  l/m (o.37/-C,;  -  yjdb-b'^  )  (b Ih) ® 


PC/,  =  2.88  +  l.nCPl{g  +  1.58  TCP 
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VCI,  Model  37 

Modified  CPI  version  20  and  TCP:  Empirical  adjustments  to  (1)  the  additional  contact 
length  due  to  sinkage  using  a  constant,  (2)  ViLC^j  using  a  constant,  and  (3)  the  contact 
area  using  aspect  ratio 


CPI*  = 


w 


«  d  [2  ^  +  c,  (cj  'ALC^,  -^jdb-¥)  Qilbf' 

Numerical  Optimization  yields  c,  =  0.51,  Cj  =  0.75  and  =  -0.50 


•••  CPl!^  = 


W 


nb  lyjdb  -  6^  +  0.51  (o.375iq,  -  ^db  -  6^  {bih) 


0.50 


FC/j  =  2.08  +  1.91  CPI +  1.89  rCF 


VCI,  Model  38 

Modified  CPI  version  21  and  TCP:  Replaced  VtLC^,  with  c^d/2  and  made  an  empirical 
adjustment  to  the  contact  area  using  aspect  ratio 


(^p  Jf  =  ^ _ _ 

nb |2 6  -  6^  +  Hm dtl  -  6  -  6^  j  j (/r lb) 

Numerical  Optimization  yields  Cj  =  0.711  and  =  -0.48 

•*.  CP7oi  =  ^ ^ _ _ _ 

nbllfdbr  6^  +  l/w  (o.356  d  -  }jdb  -  b'^)]ib  Ih) 
KC/,  =  2.75  +  1.79  CP7'^  +  1.63  TCF 
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VCI,  Model  39 

Modified  CPI  version  22  and  TCP:  Replaced  V^LC^,  with  V^LC^  (i.e.,  changed  z  fi-om 
to  c^d)  and  made  an  empirical  adjustment  to  the  contact  area  using  aspect  ratio 


Numerical  Optimization  yields  c,  =  0.16  and  -  -0.46 


FC/,  =  3.02  +  1.83  CP/ ^  +  1.51  TCP 


veil  Model  40 

Modified  CPI  version  23  and  TCP:  Replaced  V^LC^j  \wth  V2LC^^  (i.e.,  changed  z  fi’om 
to  c^h)  and  made  an  empirical  adjustment  to  the  contact  area  using  aspect  ratio 


Numerical  Optimization  yields  Cj  =  0.69  and  pj  =  -0.48 


FC/,  =  3.12  +  1.84  CPI  ^  +  1.49  TCP 
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VCIi  Model  41 

Modified  CPI  version  24  and  TCP:  Replaced  ViLC^j  wth  ViLC^  (i.e.,  changed  z  fi’om 
to  VAd)  and  made  empirical  adjustments  to  (1)  the  additional  contact  length  due  to  sinkage 
using  a  constant  and  (2)  the  contact  area  using  aspect  ratio 


'Ad  i  —-6  -»  LC,  =  d 

2  ^ 


Numerical  Optimization  yields  c,  =  0.81  and  p,  =  -0.46 


FC/,  =  3.03  +  1.84  CW '4  +  1.51  TCF 


VCI,  Model  42 

Modified  CPI  version  25  and  TCP;  Replaced  ViLC^,  vrith  'ALC^  (i.e.,  changed  z  fi’om 

to  Cjrf)  and  made  empirical  adjustments  to  (1)  using  a  constant  and  (2)  the 

contact  area  using  aspect  ratio 


CPl'^  = 


nb  2  c^yjdb  -  6^  +  1/m  ['A  LC^  -  c^^db  -  6^ )  (h  lb) 


Numerical  Optimization  yields  c,  =  0.49,  c^  =  1.4  and  p,  =  -0.48 


FC/,  =  2.75  +  2.51  CPl!^  +  1.63  TCP 
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VCIi  Model  43 

Modified  CPI  version  26  and  TCP:  Replaced  VtLC^,  with  VzLC^  (i.e.,  changed  zfi-om 
to  c^d)  and  made  an  empirical  adjustment  to  yjdb-b'^  using  aspect  ratio 


veil  Model  44 

Modified  CPI  version  27  and  TCP:  Replaced  ViLC^,  with  54 (i.e.,  changed  z  fi-om 
to  c^d)  and  made  empirical  adjustments  to  (1)  the  additional  contact  length  due  to  sinkage 
using  a  constant  and  (2)  ^d  6  -6^  using  aspect  ratio 


Numerical  Optimization  yields  c,  =  0.070,  Cj  =  2.3  and  -  -0.84 


VCI^  =  2.44  +  IMCPl'^  +  1.81  TCP 
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VCIi  Model  45 

Modified  CPI  version  28  and  TCP:  Removed  1/m,  replaced  with  VtLC^  (i.e., 

changed  z  fi’om  to  c^d\  and  made  empirical  adjustments  to  (1)  the  additional  contact 
length  due  to  sinkage  using  a  constant  and  (2)  -6^  using  aspect  ratio 

CP  Jiff  =  ^ ^ 

“  nb\z^d6  -  62  +  C2(y2LC^3-)ld^  “  (A/*)'’')] 

Numerical  Optimization  yields  c,  =  0.11,  Cj  =  0-80  and  p.^  =  -0.92 

CPlU  =  , _ _ tt 

n  A  [2  ^d6  -  b^ib  Ih)  +  0.80  {'ALC^  -  ^db  -  6*  (b  Ih)  ^ 

where  0.11</<— -6  -  LC^  =  2V4(6  +  0.11</)  -  (6  +  0.114)2 

0.114  i  — -6  -  LC^  =  d 

2  ^ 

VCI^  =  1.91  +  2.09  CPy'g  +  1.96  rCF 

veil  Model  46 
CPH^  and  TCP 

CPH^  = - ; — - 

2nb^d6-  6^  (hlbf' 

Numerical  Optimization  yields  />,  =  -0.60 

A  CPH'^  = -  ^  - 

2nbsjdb-b^  (blh)^-^ 


yCJ,  =  3.58  +  1.66  CPff'  +  0.75  TCP 
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VCIi  Model  47 
CP5'  and  TCF 


CPS' 


_ W _ 

nbLC^,(h/bf' 


Numerical  Optimization  yields  />,  =  -0.50 


CPS'  = 


_ w _ 

nbLC^jiblh)^^ 


VCI^  =  2.30  +  3.21  CPS'  +  2.84  TCF 


veil  Model  48 

Modified  Average  Contact  Pressure  Parameter  No.  2  (CP2)  version  1  and  TCF;  Used  CjO 
for  internal  length  adjustment  and  made  empirical  adjustment  to  the  contact  area  using 
aspect  ratio 


2nb  ^4(6 +  c, 4)  -  (5  +  c^df  (h lb)’’' 

Numerical  Optimization  yields  c,  =  0.026  and  /?,  =  -0.52 


2nb  ^d(b  +  0.026  4)  -  (6  +  0.026  df  (b  Ih)  ° “ 
PC/,  =  1.95  +  2.03  CP2'  +  1.95  TCP 


veil  Model  49 

Modified  CP2  version  2  and  TCF:  Split  the  contact  length  into  2  parts  and  simplified  the 

adjustment/correctionpartto  Cjt/ 

CP2i.  =  _ _ 

^  2nb[idb-b^  +  c,4)(A/i)'’‘ 

Numerical  Optimization  yields  c,  =  0.12  and  =  -0.50 

CP2i  =  _ ^ _ 

^  2nb{yjdb-b^  +  0.124)(A//i)°“ 


PC/,  =  1.86  +  2.58  CP2'  +  1.96  TCP 
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VCIi  Model  50 

Recast  CP2  [  and  TCF:  Split  the  contact  length  into  2  parts 

CP2^^  -  _ r _ : _ ^ ^ _ 

'  2n*[^rf5-62  +  (^£/(6+Cii/)-(6  +  c,rf)2  -  V^<f6  -  ^)](A/A)^' 

Numerical  Optimization  yields  =  0.026  and  =  --0.52 

•  QP2^  =  ^ 

’  2nb[^di>-b^  +  {^d(b  +  0.026 rf)  -  (6  +  0.026 df  -  ^db-  62)](fe/A)0“ 

FC/,  =  1.95  +  2.03  CP2'  +  1.95  TCF 


veil  Model  51 
and  TCF 


nbd  (A/rf)'’’  (dthY^ (b/d)”^ 

Numerical  Optimization  yields  p^  -  0.42,  p^  =  0.66  and  p^  =  0.32 

•  N 

KC/,  =  1.56  +  2.03  +  2.03  TCF 


APPENDIX  C 

SUMMARY  OF  ALL  MODEL  VARIANTS  EXPLORED 
DURING  THE  DRAWBAR  MODELING 
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Table  Cl.  Quality  comparison  for  all  drawbar  model  variants  (field  data  only) 


Model 

Model 

Degrees  of 

Standard 

Adjusted 

No. 

Description 

Freedom 

Error  (D/W) 

R-squared 

EXISTING  COMPARISON  MODELS 


Cl* 

Standard  NRMM  Model:  Excess  Rating  Cone  Index  (RCI^) 

81 

0.097 

77.8 

C2* 

Freitag’s  Clay  Numeric  ( ) 

80 

0.095 

78.7 

C3* 

Tumage’s  Clay  Numeric 

80 

0.085 

82.8 

C4® 

Smith’s  Clay  Numeric 

79 

0.082 

84.1 

C5* 

Maclaurin’s  Drawbar  Numeric  (Nj^) 

78 

0.090 

80.6 

PRELIMINARY  NEW  DEVELOPMENT  MODELS 


PI 

Bearing  Ratio  No.  1 :  Ratio  of  Soil  Strength  over  Contact 

Pressure  using  RCI  and  CP2  [ 

79 

0.120 

65.6 

P2 

Bearing  Ratio  No.  2:  Ratio  of  Soil  Strength  over  Contact 

Pressure  using  RCI  and  CPH 

81 

0.098 

77.3 

P3 

Bearing  Ratio  No.  3:  Ratio  of  Soil  Strength  over  Contact 

Pressure  using  RCI  and  CPH!^ 

80 

0.123 

64.3 

P4 

Excess  Bearing  Ratio  No.  1 :  Ratio  of  Soil  Strength  over  Contact 
Pressure  using  RCI^  and  CP2  [ 

79 

0.127 

61.7 

P5 

Excess  Bearing  Ratio  No.  2:  Ratio  of  Soil  Strength  over  Contact 
Pressure  using  RCI^  and  CPH 

81 

0.110 

71.3 

P6 

Excess  Bearing  Ratio  No.  3;  Ratio  of  Soil  Strength  over  Contact 
Pressure  using  RCI^  and  CPH^ 

80 

0.129 

60.7 

P7 

Excess  RCI  Ratio  No.  1 :  Ratio  of  RCI^  over  PC/, 

81 

0.103 

74.7 

P8 

Excess  RCI  Ratio  No.  2:  Ratio  of  RCI^  over  PC/,  with  aspect 
ratio  (product) 

81 

0.087 

82.1 

PRINCIPAL  NEW  DEVELOPMENT  MODELS 


t  A 

la 

Modified  NRMM  Model:  Optimized  Curve  Fits  to  Database 

81 

0.079 

85.3 

lb 

Modified  NRMM  Model:  Removed  CPF  Separation  Cri^a  and 
Optimized  Single  Curve  Fit  to  Database 

81 

0.084 

83.4 

2 

Excess  Bearing  Ratio  No.  4:  Ratio  of  Soil  Strength  over  Contact 
Pressure  using  RCI^  and  CP2^ 

79 

0.084 

83.2 
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Table  Cl.  (continued) 


Model 

No. 

Model 

Description 

Degrees  of 
Freedom 

Standard 
Error  (D/W) 

Adjusted 

R-squared 

3“ 

Drawbar  Numeric  No.  1  (A/p, ):  Numeric  using  bd  and  3 
dimensionless  ratios  of  tire  variables  [b/d,  d/h,  and  6/d]  wiiich 
allows  inclusion  of  b,  d,  h,  and  6  with  optimal  eicponents  under 
the  constraint  of  pressure  units 

78 

0.063 

90.5 

4 

Drawbar  Numeric  No.  2  (JVpj  );  Numeric  using  bd  and  3 

dimensionless  ratios  of  tire  variables  [b/h,  d/b,  and  6/d]  which 
allows  inclusion  of  b,  d,  h,  and  6  with  optimal  exponents  under 
the  constraint  of  pressure  units  (check  model) 

78 

0.063 

90.5 

5 

Drawbar  Numeric  No.  3  ):  Numeric  using  bd  and  2 

Himmsinnless  ratios  of  tire  Variables  [b/h  and  6/d]  which  allows 
inrliioinn  of  b,  d,  h,  and  6  under  the  constraint  of  pressure  units 

79 

90.2 

6“ 

Drawbar  Numeric  No.  4  (A/p^):  Numeric  using  bd  and  2 
dimensionless  ratios  of  tire  variables  [b/d  and  6/d]  which  allows 
inclnsinn  of  b,  d,  and  6  witii  optimal  exponents  under  the 
constraint  of  pressure  units 

79 

0.063 

90.5 

7 

Drawbar  Numeric  No.  5  Numeric  using  bd  and  1 

dimensionless  ratio  of  tire  variables  [b/d]  which  allows  inclusion 
of  b  and  d  with  optimal  exponents  under  the  constraint  of 
pressure  units 

80 

0.068 

88.9 

8 

Drawbar  Numeric  No.  6  (A/p^):  Numeric  using  d  as  the  only  tire 
variable  with  no  constraint  for  pressure  units 

80 

0.066 

89.8 

9 

Drawbar  Numeric  No.  7  (A/p,);  Numeric  using  d  as  the  only  tire 
variable  under  the  constraint  of  pressure  units 

80 

0.069 

88.9 

10® 

Drawbar  Numeric  No.  8  (A/pg);  Numeric  using  d  and  Lj  [based 
on  a  consideration  of  Drawbar  Model  6] 

80 

0.063 

90.6 

11 

Drawbar  Numeric  No.  9  (N^g):  Numeric  with  RCI^  for  soil 
strength  and  using  d  and 

80 

0.086 

82.3 

12 

Drawbar  Numeric  No.  10  (A/p,g):  Numeric  using  Lg  as  the  only 
tire  parameter 

80 

0.067 

89.4 

13 

Drawbar  Numeric  No.  1 1  ):  Numeric  using  d  ^d  a 

modified  version  of  where  an  adjustment/correction  is  made 
using  c^d  in  the  middle  ordinate 

79 

0.063 

90.5 

14 

Drawbar  Numeric  No.  12  (A/p„):  Numeric  using  d  and  Lj  with 
TCP 

79 

0.063 

90.5 

15 

Drawbar  Numeric  No.  13  (A/p,g):  Numeric  using  d  and  Lj  with 
aspect  ratio  for  a  width  consideration 

79 

0.063 

90.5 

16 

Drawbar  Numeric  No.  14  (A/p,^);  Numeric  using  d  and  L  j  with 
b/d  for  a  width  consideration 

79 

0.063 

90.5 
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Table  Cl.  (continued) 


Mode] 

Na 

Model 

Description 

Degrees  of 
Freedom 

Standard 
Error  (D/W) 

Adjusted 

R-$quared 

17 

Drawbar  Numeric  No.  15  (Ar^/5)-  Numeric  using  d  and  with 
^jbld  for  a  width  consideration  [based  on  Freitag’s  laboratory 
data] 

79 

0.076 

86.4 

18* 

Drawbar  Numeric  No.  16  Numeric  using  d  and  with 

{hid  f  [based  on  Freitag’s  laboratory  data) 

79 

0.064 

90.2 

-  These  models  were  discussed  in  the  main  text 


Drawbar  Model  Cl 

Standard  NRMM  Model:  Excess  Rating  Cone  Index  (RCl^) 


RCl^  =  RCI  -  VCI^ 


CPF  i  Apsi 


CPF<  4  psi 


—  ^  =  0.621  - - — - 

RCI^*1.6\ 

=  0.678  - . -Ijg- . 

IV^  RCI^  +  6.80 


Drawbar  Model  C2 
Freitag’s  Clay  Numeric 

„  RCInbd(b\'^ 

O^.0.74-— 
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Drawbar  Model  C3 
Tumage’s  Clay  Numeric 


Drawbar  Model  C4 
Smith’s  Clay  Numeric 


=  0.72 - 2JL_ 

»^20  ^c(S)^7.4 


Drawbar  Model  C5 

Maclaurin’s  Drawbar  Numeric  (N^^) 


N, 


M 


RCInb°^d^^b°* 

w 


DRCI 


W. 


20 


0.73  - 


2.7 

Ar^.3.2 


-3/4 
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Drawbar  Model  P 1 

Bearing  Ratio  No.  1 :  Ratio  of  Soil  Strength  over  Contact  Pressure  using  RCI  and  CP2  [ 


BR^  = 


RCI 

CP2[ 


BR,  = 


RCI 


W 


[inb  }Jd(6  +  0.026d)  -  (6  +  0.026df  (*//r)® 


7  5 

—  =  0.76  - 


BR^  + 10.5 


Drawbar  Model  P2 

Bearing  Ratio  No.  2:  Ratio  of  Soil  Strength  over  Contact  Pressure  using  RCI  and  CPH 


BR^ 


RCI 

CPH 


bR^  =  j 


=  0.74  - — — 

BR^^A.l 


Drawbar  Model  P3 

Bearing  Ratio  No.  3:  Ratio  of  Soil  Strength  over  Contact  Pressure  using  RCI  and  CPH'^ 


BR, 


RCI 

CPH'^ 


bR^  =  2nb^db-b^  (blHf  ^  j 


6  2 

id  =  0.76  -  — — — 
BR^  +  8.9 
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Drawbar  Model  P4 

Excess  Bearing  Ratio  No.  1:  Ratio  of  Soil  Strength  over  Contact  Pressure  using  RCI^  and 
CP2\ 


RCI 


Drawbar  Model  P5 

Excess  Bearing  Ratio  No.  2:  Ratio  of  Soil  Strength  over  Contact  Pressure  using  RCI^  and 
CPH 


BR^ 


RCI^ 

"CPH 


pRC/ 


0.69 


1.79 

B/?^  +  2.60 


Drawbar  Model  P6 

Excess  Bearing  Ratio  No.  3:  Ratio  of  Soil  Strength  over  Contact  Pressure  using  RCl^  and 
CPH', 


^20 


0.69 


2.35 

BR^  +  3.40 
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Drawbar  Model  P7 

Excess  RCI  Ratio  No.  1:  Ratio  of  RCI^  over  VCl^ 


RCIR^, 


VCI^ 


RCIR^j 


RCI  -  VCI^ 

raj 


^"°  =  0.68- _ 

RCIR^J  +  1.20 


Drawbar  Model  P8 

Excess  RCI  Ratio  No.  2:  Ratio  of  RCl^  over  VCl^  with  aspect  ratio  (product) 


—  =  0.66  - - 9:^9 _ 

RCIR^  + 0.60 


Drawbar  Model  la 

Modified  NRMM  Model;  Optimized  Curve  Fits  to  Database 


CPFi  4psi 


CPF<4psi 


=  0.63  - 

IV  20  RCI^ +  S.16 

=  0.81  - 

W20  RCI^  + 22.0 


Drawbar  Model  lb 

Modified  NRMM  Model:  Removed  CPF  Separation  Criteria  and  Optimized  Single  Curve 
Fit  to  Database 

—  =  0.65 - _ 

W20  RCI^  + 12.0 
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Drawbar  Model  2 

Excess  Bearing  Ratio  No.  4:  Ratio  of  Soil  Strength  over  Contact  Pressure  using  RCl^  and 
CP2'^ 


BR.4 


RCJ^ 

CP2^ 


where  CP2^  - 
Numerical  Optimization  yields  c,  =  0  and  p,  =  -2 


BR^^  =  ^—^^2nb  ^jdb-fP-  Qtlb)^ 


=  0.65 - 

0.65 


Drawbar  Model  3 

Drawbar  Numeric  No.  1  {Nj^, ):  Numeric  using  bd  and  3  dimensionless  ratios  of  tire 
variables  [b/d,  d/h,  and  6/d]  which  allows  inclusion  of  b,  d,  h,  and  8  with  optimal 
exponents  imder  the  constraint  of  pressure  units 

b  d  ibldf^  (dlhf^  (6/rf)'»  j 

Numerical  Optimization  yields  =  -0.90,  p^  =  0.82,  and  p^  =  0.45 
.  A/  -  /?C/f 

“  Dl  ^  ^  ^0.82  J 


0.695  - 


8.6 

N^,*7.5 
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Drawbar  Model  4 

Drawbar  Numeric  No.  2  Numeric  using  bd  and  3  dimensionless  ratios  of  tire 

variables  [b/h,  d^,  and  6/d]  which  allows  inclusion  of  b,  d,  h,  and  6  with  optimal 
exponents  under  the  constraint  of  pressure  units  (check  model) 

^  ] 

Numerical  Optimization  yields  =  0.82,  /?2  =  1.72  and  p^  =  0.45 
.  ^  _  RCl{nb^^^d^^l_b^ 

**  “  W  [  /,0*2  ^ 


DRCi 

—  =  0.695 


8.6 


W. 


20 


Drawbar  Model  5 

Drawbar  Numeric  No.  3  Numeric  using  bd  and  2  dimensionless  ratios  of  tire 
variables  [b/h  and  6/d]  which  allows  inclusion  of  b,  d,  h,  and  6  under  the  constraint  of 
pressure  units 


N, 


D3 


RCI 


w 


^nbdib/hf'ibfd)”^^ 


Numerical  Optimization  yields  p^  =  -1.10  and  p^  =  0.60 


••  ^D3  ~ 


RCI 


w 


nh 


1.10  j  0.40  x0.60 


,0.10 


.  0.68S  - 

»'»  Wo,.  0.344 
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Drawbar  Model  6 

Drawbar  Numeric  No.  4  (N^):  Numeric  using  bd  and  2  dimensionless  ratios  of  tire 
variables  [b/d  and  8/d]  which  allows  inclusion  of  b,  d,  and  6  with  optimal  exponents  under 
the  constraint  of  pressure  units 


d  (b/d)”'  (m”^ ) 

Numerical  Optimization  yields  />]  =  -1.0  and  P2  =  0.50 


N, 


D4 


N, 


D4 


DRCI 


W. 


20 


0.69  - 


2.55 

Ar^  +  2.10 


Drawbar  Model  7 

Drawbar  Numeric  No.  5  (N^^y.  Numeric  using  bd  and  1  dimensionless  ratio  of  tire 
variables  [b/d]  which  allows  inclusion  of  b  and  d  with  optimal  exponents  under  the 
constraint  of  pressure  units 

Numerical  Optimization  yields  /»,  =  -1.05 
^0.05  J 


^D5  ~ 


RCI 

w 


^RCi 


0.694 


12.3 

N^s  +  11-4 
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Drawbar  Model  8 

Drawbar  Numeric  No.  6  Numeric  using  d  as  the  only  tire  variable  with  no 
constraint  for  pressure  units 


Ntmerical  Optimization  yields  p,  =  1.64 


DRCI 


W. 


20 


0.687 


2.71 

Noe  ^2.00 


Drawbar  Model  9 

Drawbar  Numeric  No.  7  Numeric  using  d  as  the  only  tire  variable  under  the 
constraint  of  pressure  units 

— =  0.695  - — 

IFjo  ^i,7  +  10.9 


Drawbar  Model  10 

Drawbar  Numeric  No.  8  (N^g):  Numeric  using  d  and  Lg  [based  on  a  consideration  of 
Drawbar  Model  6] 


where  ij  =  2^d5  -  6^ 

nflC/  5  Q 

il  =  0.692  - - — — 

^20  Nos  1 
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Drawbar  Model  11 

Drawbar  Numeric  No.  9  (N^g):  Numeric  with  RCI^  for  soil  strength  and  using  d  and  i  j 

RCL I  \ 

^D9  =  ^  j 

RCI  -  VCL  (  y  \ 

- -<'"«") 

3  41 

—  =  0.672  -  — — — 


Drawbar  Model  12 

Drawbar  Numeric  No.  10  (N^,o):  Numeric  using  as  the  only  tire  parameter 


—  =  0.698  - - - 

W'jo  No,o*2.30 


Drawbar  Model  13 

Drawbar  Numeric  No.  1 1  (N^„  ):  Numeric  using  d  and  a  modified  version  of  where  an 
adjustment/correction  is  made  using  c^dm  the  middle  ordinate 


where  ij  =  2^rf(6  +  c,rf)  -(6  +  Cjrf)^ 
Numerical  Optimization  yields  =  -0.003 


'  N 

( Ind^dib -  0.003<i)  -  (6 - O.OOSrf)"' 

W 

V 

DRCI 

—  =  0.692  - 

4.86 

^20 

A/o„+4.00 

Drawbar  Model  14 

Drawbar  Numeric  No.  12  Numeric  using  d  and  Zj  with  TCF 
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^D,2  =  ^\ndL,il*cJCF)] 
Numerical  Optimization  yields  Cj  =  -0.015 


RCI 


DI2 


W 


2n  (1  -  0.015  TCF) 


nRCi  4  9 

—  =  0.691  - 

W. 


20 


^d/2-^4.0 


Drawbar  Model  15 

Drawbar  Numeric  No.  13  Numeric  using  d  and  Zj  wdth  aspect  ratio  for  a  width 

consideration 


Numerical  Optimization  yields  /?,  =  -0.005 


.  0.692  -  ’•« 

w. 


20 


^d/5-4.1 
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Drawbar  Model  16 

Drawbar  Numeric  No.  14  Numeric  using  d  and  Xj  vdth  b/d  for  a  width 

consideration 


Numerical  Optimization  yields  =  0.0014 


Drawbar  Model  17  _ 

Drawbar  Numeric  No.  15  Numeric  using  d  and  Ij  with  {bid  for  a  width 

consideration  [based  on  Freitag’s  laboratory  data] 


Drawbar  Model  18 

Drawbar  Numeric  No.  16  Numeric  using  d  and  Xj  with  {hldf^  [based  on 

Freitag’s  laboratory  data] 


APPENDIX  D 

TYPICAL  PROCEDURES  FOR  MEASURING 
VCI,  AND  DRAWBAR  PERFORMANCE 
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Multi-Pass  Experiments 

VCIj  is  not  directly  measurable.  Rather  it  is  interpreted  from  the  results  of  multi-pass 
experiments.  In  a  typical  multi-pass  experiment,  a  level,  straight-line,  homogeneous  (as  much 
as  possible)  test  lane  is  marked  off,  and  cone  index  measurements  are  taken  throughout  the 
lane.  Then  the  vehicle  makes  passes  through  the  test  lane  at  a  slow,  steady-state  speed 
(approximately  2  mph)  in  its  lowest  gear.  The  experiment  is  usually  conducted  with  the 
vehicle  first  traversing  forward  through  the  test  lane  for  pass  number  one  and  then  traversing 
backward  (i.e.,  in  reverse  gear)  through  the  test  lane  for  pass  number  two.  The  vehicle  will 
continue  to  make  passes  until  immobilization  occurs.  When  immobilization  is  reached,  the 
immobilization  pass  number  is  recorded,  and  other  supporting  soil  consistency  data  (e.g., 
remold  index,  moisture  content,  density,  etc.)  are  measured  in  a  spot  adjacent  to  the  immobili¬ 
zation  but  out  of  the  zone  of  disturbance.  The  measured  soil  consistency  data  are  intended 
to  represent  the  soil  characteristics  with  the  greatest  influence  on  the  immobilization  for  the 
particular  vehicle  and  terrain  conditions. 

To  establish  the  VCIj  measurement  for  a  particular  vehicle  configuration,  several 
multi-pass  experiments  are  conducted  on  ranging  soil  strengths  in  order  to  acquire  several 
observations  for  passes-made-good  versus  soil  strength  in  the  critical  layer.  The  critical  layer 
is  the  layer  of  soil  (typically  3-9"  or  6-12")  that  contributes  most  significantly  to  the  perform¬ 
ance  of  the  vehicle.  The  soil  strength  value  at  which  the  vehicle  is  capable  of  consistently 
maldng  one  pass  is  then  interpreted  from  a  x-y  plot  of  the  collected  multi-pass  observations, 
and  this  inferred  soil  strength  represents  the  VCIi  performance  measurement.  Figure  D1 
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provides  a  sample  x-y  plot  of  multi-pass  observations  that  were  used  to  determine  the  VCIj 
measurement  for  one  of  the  vehicle  configurations  used  during  this  research. 


Figure  Dl.  Sample  multi-pass  observations  used  for  determining  VCIj 

Drawbar  Experiments 

A  drawbar  experiment  is  conducted  to  determine  the  drawbar  performance  as  a 
function  of  slip  for  a  particular  vehicle  configuration  and  soil  condition  (i.e.,  type,  strength, 
surfece  wetness,  etc.).  In  the  typical  drawbar  experiment,  a  level,  straight-line,  homogeneous 
(as  much  as  possible)  test  lane  is  marked  off,  the  lane  is  prepared  with  simulated  rainfall  (i.e., 
surface  wetness)  or  left  dry  (i.e.,  normal)  depending  on  the  desired  test  condition,  and  then 
soil  consistency  data  are  measured  in  the  lane,  out  of  the  zone  of  passage  for  the  traction 
elements.  The  test  vehicle  is  connected  to  a  load  vehicle  with  a  tension  member  (e.g.,  steel 
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cable,  nylon  strap,  etc.),  and  a  force  load  ceU  is  connected  between  the  vehicles  within  the 
tension  member.  The  test  vehicle  is  also  instrumented  so  that  the  actual  horizontal  distance 
traveled  and  the  ^parent  horizontal  distance  traveled,  as  indicated  by  the  rotational  velocity 

of  the  traction  elements,  can  be  measured. 

The  test  vehicle  enters  the  lane  followed  by  the  load  vehicle  at  a  slow,  steady-state 
speed  (approximately  2  mph).  The  load  vehicle  controls  the  experiment  by  controlling  the 
force  in  the  connection  load  cell  (i.e.,  drawbar  force).  It  applies  braking  force  in  stepwise 
increments  so  that  the  test  vehicle  undergoes  a  stepwse  transition  from  a  zero  pull,  zero  slip 
condition  to  a  high  pull,  high  slip  condition.  The  test  vehicle  simply  attempts  to  maintain 
constant  motion  by  applying  more  horsepower  as  necessary.  The  measured  data  resulting 
from  a  drawbar  experiment  are  principally  drawbar  force,  slip,  and  soil  consistency. 

In  order  to  reduce  the  measured  data  from  a  drawbar  experiment  into  a  single  drawbar 
performance  quantity,  a  x-y  plot  of  drawbar  coefficient  versus  slip  is  generated  for  the 
particular  vehicle  configuration  and  soil  condition.  Then  a  best-fit  curve  (and  possibly  an 
accompanying  equation)  is  established  to  describe  the  trend  using  visual  or  statistical 
methods.  This  best-fit  curve  is  used  to  determine  the  expected  drawbar  coefficient  at  a  single 
slip  magnitude.  The  slip  magnitude  typically  used  is  either  20%  or  the  true  optimum  slip 
which  normally  occurs  between  10%  and  30%  for  fine-grmned  soils.  If  optimum  slip  is  used, 
it  is  determined  by  evaluating  the  work  index  versus  slip  trend  in  addition  to  the  drawbar 
versus  slip  trend,  and  the  optimum  sUp  is  that  occurring  at  the  maximum  work  index.  The 
inferred  drawbar  coefficient  at  20%  or  optimum  slip  represents  the  drawbar  performance 
measurement.  Figure  D2  provides  a  sample  x-y  plot  of  data  that  were  used  to  determine  a 
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drawbar  measurement  for  one  of  the  vehicle  configurations  used  during  this  research.  Note 
that  there  can  be  several  drawbar  measurements  for  a  single  vehicle  configuration  since 
drawbar  experiments  can  (and  often  should)  be  conducted  on  several  different  soil  strengths 
and/or  surface  conditions  depending  on  the  goals  of  the  particular  study. 


Figure  D2.  Sample  drawbar  data 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reportingburdenformis  collection 

Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. - - - - - - 


111.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

August  1999 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final  report 


4.  TITLE  AND  SUBTITLE 

Improving  the  Traction  Prediction  Capabilities  in  the  NATO  Reference 
Mobility  Model  (NRMM) 

6.  AUTHOR(S) 

Jody  D.  Priddy 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Waterways  Experiment  Station 
3909  Halls  Ferry  Road,  Vicksburg,  MS  39180-6199 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

Technical  Report  GL-99-8 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

Available  from  National  Technical  Infomation  Service,  5285  Port  Royal  Road,  Springfield,  VA  22161. 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  fMaxi/nom  200  worefsj 

The  NATO  Reference  Mobility  Model  (NRMM)  has  traction  relationships  that  describe  the  performance  of  vehicles  in 
naturally  occurring  soils  reasonably  well,  but  they  have  weaknesses  stemming  primarily  from  an  inadequate  consideration  of 
the  influence  of  vehicle  characteristics  on  traction.  This  research  sought  to  demonstrate  the  feasibility  of  developing  better 
traction  relationships  for  NRMM  through  exploitation  of  existing  data.  In  order  to  demonstrate  the  feasibility,  a  statistical 
modelin<'  effort  was  conducted  for  developing  new  traction  relationships.  The  research  involved  a  systematic  and  extensive 
exploration  of  new  and  old  traction  modeling  concepts  that  rely  on  relatively  simple  characterizations  of  the  vehicle  and  the 
soil  The  effort  concentrated  on  one-pass  vehicle  cone  index  (VCIl)  and  drawbar  performance  of  pneumatic-tired  vehicles  on 
nonslippery,  highly  plastic  clay  soils  (CH  classification).  The  results  demonstrated  that  better  traction  relationships  can  be 
developed  for  NRMM  and  provided  new  insight  into  the  controlling  mechanisms  in  VCIl  and  drawbar  performance. 


1 15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


14.  SUBJECT  TERMS  r  « 

Drawbar  pull  Mobility  modeling  NRMM  _ 215  _ 

Pneumatic  tire  Traction  Vehicle  Cone  Index  16.  PRICE  CODE 

Wheeled  vehicle 

17.  SECURITY  CLASSIFICATION  1 18.  SECURITY  CLASSIFICATION  1 19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


||  UNCLASSIFIED 

NSN  7540.01-280-5500 


UNCLASSIFIED 


Standard  Form  298  (Rev.  2-i 

Prescribed  by  ANSI  Std,  Z39-18 
298*102 


